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Summary 

The  possibility  of  predicting  the  ultimate  static  bearing 
capacity  of  a  pile  from  its  dynamic  behavior  during  driving  is 
investigated  by  the  wave  equation  method  of  pile  driving  analysis. 
The  method  developed  by  E.  A.  L.  Smith  over  a  period  of  years  is 
followed,  but  its  validity  is  first  shown  by  proving  that  the  math- 
ematical model  used  is  equivalent  to  the  wave  equation.  Computer 
programs  developed  for  purposes  of  this  investigation  are  presented 
in  their  entirety  together  with  detailed  instructions  on  their  use 
so  that  they  may  be  utilized  by  future  investigators  who  may  not  be 
familiar  with  computer  work. 

The  soil  engineering  aspects  of  the  problem  are  explored.  The 
information  available  on  dynamic  soil  properties  from  published  test 
data  is  reviewed  to  obtain  information  applicable  to  pile  driving. 
Using  values  thus  deduced  and  the  computer  programs,  computations  of 
ultimate  resistance  are  made  from  2k   published  pile  driving  records 
and  correlated  with  their  load  test  results. 

Although  not  conclusive  because  of  the  relatively  small  num- 
ber of  correlations  attempted,  the  results  are  very  encouraging.  It 
appears  that  the  wave  equation  method  of  pile  driving  analysis  may 
become  an  accurate  tool  in  predicting  a  pile's  static  bearing  capacity 
from  its  driving  record.  This  appears  to  be  the  case  regardless  of 
the  type  or  size  of  pile  and  driving  equipment,  and  indications  are 
that  predictions  may  be  possible  for  piles  driven  in  cohesive  as  well 
as  granular  soils.  General  limitations  of  the  method  are  summarized. 


XI 


Values  deduced  for  point  and  side  damping  and  ground  quake 
in  various  soils  are  presented.  Corresponding  values  of  friction 
acting  on  the  pile  sides  as  a  percentage  of  the  ultimate  resistance 
are  given  for  the  piles  investigated.  More  research  into  this  problem 
is  indicated,  and  recommendations  are  made  for  such  work. 

It  is  concluded  that  even  with  inexact  knowledge  of  dynamic 
soil  properties  the  wave  equation  method  of  pile  driving  analysis 
appears  to  give  good  results. 
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Chapter  I 
INTRODUCTION 

The  possibility  of  relating  the  behavior  of  a  pile  during 
driving  to  its  static  load  carrying  capacity  by  the  use  of  the  wave 
equation  method  of  pile  driving  analysis  is  encouraging.  Many 
engineers  today  realize  that  the  numerous  dynamic  pile  driving 
formulas  which  have  been  used,  and  are  still  being  used,  have  seri- 
ous limitations  and  cannot  be  depended  upon  to  give  reliable  results, 
but  they  are  still  used  for  lack  of  an  adequate  substitute.  Some 
feel  that  the  dynamic  formula  approach  is  not  valid,  and  have  turned 
instead  to  trying  to  determine  the  bearing  capacity  of  a  pile  from 
purely  static  soil  mechanics  and  structural  considerations  for  the 
cases  where  the  pile-soil  system  is  required  to  carry  only  static 
loads.  Such  an  approach  is  certainly  valid,  and  it  requires  a  thorough 
investigation  of  the  soil  conditions,  together  with  a  careful  evalua- 
tion of  their  engineering  properties,  and  the  use  of  a  great  deal  of 
personal  judgment  in  its  application. 

It  is  believed  that  the  wave  equation  method  of  analysis  of 
the  dynamic  pile  driving  problem  is  an  equally  valid  approach  and  may 
offer  advantages  over  the  static  approach  if  it  can  be  successfully 
applied.  This  is  not  to  say  that  one  approach  should  be  used  to  the 
exclusion  of  another,  but  that  they  are  working  toward  a  common  goal 
and  should  supplement  each  other.  The  wave  equation  method  is  not 
just  another  dynamic  pile  driving  formula,  but,  rather,  it  is  a 
method  of  analysis  which  is  well  founded  mathematically.  It  requires 


a  knowledge  of  the  soil  and  its  properties,  both  static  and  dynamic, 
for  its  successful  application.  Current  knowledge  of  soil  dynamics 
is  incomplete,  but  an  increasing  amount  of  work  is  being  done  in  this 
field. 

Approaching  the  problem  from  the  dynamic  side,  the  over-all 
consideration  is  to  relate  the  dynamic  behavior  of  the  driving  equip- 
ment-pile-soil system  to  the  static  behavior  of  the  structure-pile- 
soil  system  with  due  consideration  to  changes  in  properties  of  the 
soil  during  and  after  driving  of  the  piles.  A  part  of  this  broader 
problem  is  to  relate  the  dynamic  resistance  to  driving  of  a  single 
pile  under  the  last  hammer  blow  to  its  static  bearing  capacity,  with 
time  effects  between  driving  and  static  loading  minimized.  It  is  for 
this  more  limited  problem  that  the  wave  equation  method  of  analysis 
is  useful. 

The  wave  equation  alone  will  not  provide  a  solution  to  the  more 
general  consideration  which  includes  the  group  effect  of  piles,  de- 
velopment of  negative  friction,  consolidation  of  a  clay  layer  beneath 
the  pile  tip,  long-term  changes  in  water  table,  deterioration  of  the 
pile  due  to  insect  attack  or  deleterious  chemicals,  set  up  or  relaxa- 
tion of  the  soil  with  time,  or  development  of  hydrostatic  uplift. 
That  it  does  not  provide  a  solution  to  such  problems  may  be  obvious, 
but  the  statement  is  made  here  in  an  effort  to  clear  up  the  type  of 
misunderstandings  and  objections  already  voiced  in  opposition  to  the 
application  of  the  wave  equation  to  pile  driving  analysis.  In  short, 
it  is  not  a  magical  formula,  but  simply  a  tool  which  is  well  founded 


mathematically  and  is  available  to  assist  in  visualizing  and  evaluating 
the  dynamics  of  the  problem. 

The  wave  equation  can  be  utilized  to  advantage  for  investigation 
of  other  facets  of  pile  driving  such  as  stresses  in  the  pile  and 
selection  of  appropriate  pile  driving  equipment  for  a  particular  set 
of  field  conditions ;   but  the  main  points  of  this  investigation  will 
be  to  consolidate  the  wave  equation  theory  as  applicable  to  pile  driv- 
ing, make  readily  available  computer  programs  to  facilitate  the  work 
of  future  investigators,  explore  the  unsolved  problems  of  interaction 
between  soil  and  pile,  attempt  correlation  between  wave  equation  solu- 
tions and  a  few  pile  driving  records  and  load  tests,  and  to  suggest 
possibilities  for  further  research  work  in  this  general  area. 


Chapter  II 
BACKGROUND 

Although  use  has  been  made  of  piles  for  at  least  2,000  years, 
attempts  to  determine  the  bearing  capacity  of  a  pile  from  its  driving 
record  have  been  generally  unsuccessful  in  spite  of  the  amount  of 
effort  spent  in  this  direction.  An  extensive  treatment  of  the  develop- 
ment and  interrelatedness  of  dynamic  formulas  by  Chellis  (ref.  5)  is 
included,  and  he  shows  that  they  may  be  grouped  as  follows : 

The  first  dynamic  formula  was  proposed  by  Major  Saunders  in 
1851  and  is  perhaps  the  simplest  and  most  direct.  He  equated  the 
weight  of  the  ram  multiplied  by  the  stroke,  with  the  driving  resistance 
multiplied  by  the  set,  and  applied  a  factor  of  safety  of  8. 

Another  group  of  proposed  formulas  contained  a  fixed  coefficient 
which  was  supposed  to  compensate  in  some  degree  for  factors  affecting 
the  results  which  were  not  included  as  terms  in  the  formulas.  The 
Engineering  News,  Wellington,  Vulcan,  and  Bureau  of  Yards  and  Docks 
(not  now  used  by  Navy)  formulas  are  of  this  type. 

Another  group  of  formulas  attempts  to  account  for  the  variables 
by  using  expressions  for  efficiency  of  applied  energy  by  including 
relative  weights  of  the  pile  and  the  hammer.  The  Dutch,  Ritter,  and 
Benabencq  formulas  are  of  this  type. 

Other  similar  formulas  try  to  include  the  effect  of  variables 
by  using  both  fixed  coefficients  and  expressions  for  the  relative 
weights  of  the  pile  and  ram.  The  Eytelwein  and  Navy-McKay  (not  now 
used  by  the  Navy)  are  of  this  type. 


■ 


• 


■ 


The  next  group  contains  either  all  or  part  of  a  series  of 
terms  designed  to  represent  impact  losses  in  the  driving  cap,  soil, 
and  pile  during  driving.  The  Redtenbacher  and  Hiley  formulas  are 
among  those  which  fall  into  this  category. 

Only  a  few  of  the  formulas  proposed  have  been  mentioned. 
Chellis  (ref.  5)  lists  38,  while  it  is  understood  that  the  Editors  of 
Engineering  News  Record  have  on  file  U50  such  formulas.  Some  of  these 
formulas  are  extremely  simple,  while  others  are  quite  complicated. 
Some  try  to  approximate  the  dynamic  situation,  while  others  are  con- 
trived from  a  purely  statistical  approach  without  reference  to  the 
numerous  variables  actually  involved.  For  example,  Marvin  Gates  has 
proposed  (ref.  lU)  that  for  a  pile  driven  with  a  Vulcan  Number  1  Hammer 
the  ultimate  bearing  capacity  may  be  determined  by  multiplying  U8  by 
the  log  of  (10/set  in  inches).  He  claims,  ".  .  .this  relationship 
gives  more  consistently  accurate  results  than  the  most  complex  dynamic 
formula  yet  advanced,"  in  spite  of  the  fact  that  it  contains  but  two 
simple  parameters.  That  this  can  be  true  for  a  formula  which  seems 
to  ignore  the  basic  aspects  of  the  problem  is  not  surprising  when  one 
considers  the  shortcomings  and  omissions  of  -the  other  dynamic  formulas 
proposed. 

Some  of  the  dynamic  formulas  were  based  on  Newtonian  theory 
of  impact  which  has  been  shown  by  A.  E.  Commings  (ref.  10)  as  in- 
applicable to  the  pile  driving  problem.  Other  formulas  subtracted  the 
same  energy  losses  twice,  as  also  shown  by  Cummings  (ref.  10).  The 
Committee  on  the  Bearing  Value  of  Pile  Foundations  of  the  American 
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Society  of  Civil  Engineers  recognized  the  many  limitations  of  dynamic 
pile  driving  formulas  (ref.  8).  After  this  point  of  agreement,  en- 
gineers seem  to  be  divided  into  two  general  schools  of  thought.  One 
group  is  optimistic  and  feels  that  rational  use  can  be  made  of  such 
formulas.  The  other  group  stresses  the  limitations  of  the  formulas, 
and  basically  takes  the  position  that,  since  none  of  them  are  reliable, 
you  may  as  well  save  yourself  work  and  choose  a  simple  formula.   The 
latter  school  of  thought  may  seem  rash  at  first  readingj  but  it  does 
have  some  merit,  for  it  is  simply  recognizing  that  all  of  the  dynamic 
formulas  are  over-simplifications  of  a  very  complex  problem. 

It  is  certain  that  each  of  the  k$0   formulas  available  to 
choose  from  will  give  reliable  results  under  particular  combinations 
of  driving  equipment  and  soil  conditions  for  a  particular  pile,  but 
there  would  seem  to  be  an  almost  infinite  number  of  combinations  of 
the  variables  involved  making  the  choice  of  a  formula  a  difficult  one. 
That  so  many  formulas  have  been  developed  is  indicative  of  the  interest 
in  the  problem  and  the  ingenuity  of  engineers  in  attempting  to  estab- 
lish a  relationship  between  static  and  dynamic  resistance  of  piles. 
Empirical  formulas  certainly  have  their  place  when  their  range  of 
applicability  is  clearly  established  and  when  their  limitations  are 
kept  clearly  in  mind;  however,  this  has  not  always  been  the  case  in 
pile  driving. 

In  spite  of  the  interest  exhibited,  there  has  not  been  a  great 
deal  of  work  done  theoretically  or  experimentally  in  defining  or 
attacking  the  fundamental  problems  involved.  The  complexity  of  the 
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problem,  the  difficulty  in  isolating  the  variables,  and  the  limited 
techniques  and  equipment  available  for  detailed  analysis  have  generally 
limited  advances  in  this  field.  Full  utilization,  however,  of  the 
facilities  available  has  not  been  realized. 

The  basic  problem  is  to  relate  the  dynamic  resistance  to 
driving  of  a  single  pile  to  its  static  bearing  capacity  with  time 
effects  between  driving  and  static  loading  minimized.  Available  to 
assist  in  the  analysis  of  this  problem  is  the  wave  equation  which  was 
developed  over  a  hundred  years  ago  by  De  Saint  Venant  and  Boussinesq 
for  end  impact  on  rods.  It  is  well  founded  mathematically.  D.  V. 
Issaccs,  in  1931,  was  the  first  to  point  out  that  wave  action  occurred 
during  the  driving  of  piles.  In  1938,  E.  N.  Fox  published  a  solution 
to  the  wave  equation  applied  to  pile  driving,  but,  as  no  electronic 
computers  were  available  at  that  time,  he  made  a  number  of  simplifying 
assumptions.  Again  in  19U0  and  19U1,  A.  E.  Cummings  reported  on  the 
work  of  the  foregoing  investigators.  Ten  years  elapsed  before  the  wave 
equation  was  again  proposed  as  being  applicable  to  pile  driving  analysis, 
this  time  by  E.  A.  L.  Smith  (ref.  27).  Then  in  1961  Smith  published  a 
more  comprehensive  article  on  application  of  the  wave  equation  to  pile 
driving  and  described  a  numerical  approach  together  with  a  description 
of  his  efforts  in  developing  a  computer  program  for  the  solution  (ref. 
28).  His  efforts  have  been  untiring  in  developing  this  approach  and 
have  opened  the  door  for  much  needed  work  in  this  field. 

The  wave  equation,  as  will  be  shown,  mimics  the  behavior  of 
a  pile  and  provides  a  good  indication  of  what  will  happen  to  the  pile 


under  conditions  of  intact  for  specific  boundary  conditions.  For  this 
reason,  a  distinction  is  made  between  the  wave  equation  and  dynamic 
formulas  which  appear  to  be  mathematical  but  are  in  reality  empirical. 
The  present  state  of  knowledge  up  through  Smith's  work  can  provide 
a  prediction  of  effects  caused  by  a  ram  hitting  a  pile.  Still  to  be 
investigated  is  the  nature  of  the  effects  caused  by  interaction  be- 
tween the  pile  and  the  soil.  This,  too,  is  a  dynamic  problem,  and 
it  lies  in  the  new  field  of  soil  dynamics.  This  matter  will  be  ex- 
plored, but,  first,  the  wave  equation  theory  and  its  numerical  solution 
wiU  be  developed. 


Chapter  III 
DEVELOPMENT  OF  BASIC  EQUATIONS  AND  THEORY 

A.  Theory 

A.  E.  L.  Smith,  over  a  period  of  several  years,  derived  equa- 
tions suitable  for  numerical  calculations  and  computer  adaptation 
for  pile  driving  analysis  by  the  wave  equation.  His  derivations  are 
based  on  a  mathematical  model  using  weights  and  springs.  This  work 
has  generally  been  accepted  as  correct  and  as   equivalent  to  the  wave 
equation, ,  but  the  full  derivation  establishing  this  equivalence  has 
not  been  published.  In  order  to  establish  this  equivalence  and  the 
validity  of  the  numerical  method,  a  step-by-step  derivation  from  first 
principles  is  presented.  For  completeness,  Smith's  derivations  from 
the  mathematical  model  are  included. 

The  general  plan  for  establishing  this  link  is  as  follows: 

a.  A  partial  differential  equation  is  developed  and  is  shown 
to  be  a  form  of  the  wave  equation. 

b.  This  equation  is  converted  into  a  difference  equation  for 
ease  of  numerical  solution. 

c.  To  simplify  further  numerical  solution,  a  mathematical 
model  consisting  of  weights  and  springs  is  assumed  as  being  equivalent 
to  the  elastic  rod  used  as  a  basis  for  steps  a  and  b  above.  Five 
equations  are  developed  from  this  mathematical  model  which  is  based 

on  the  assumption  that  all  springs  are  perfectly  elastic, and  the  pile 
is  represented  typically  as  shown  by  figure  (3). 
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d.  It  is  then  shown  that  the  equations  of  step  c. can  be  com- 
bined to  form  the  same  equation  developed  in  step  b.  above,  and  it  is 
concluded  that  the  mathematical  model  and  the  five  equations  developed 
from  it  are  equivalent  to  the  elastic  rod,  and  the  wave  equation  is 
applicable  to  it. 

Beginning  with  development  of  the  partial  differential  equation, 
as  shown  by  Timoshenko  and  Goodier  (ref.  32),  consider  the  longitudinal 
waves  in  a  prismatic  bar  of  constant  cross  section  such  as  shown  in 
figure  (l).  The  axis  of  the  bar  is  taken  as  the  x-axis.  It  is  assumed 
that: 

a.  Cross  sections  of  the  bar  remain  plane  during  deformation. 

b.  The  unit  elongation  at  any  cross  section  mn,  due  to  a 

longitudinal  displacement  D,  is  equal  to  -  D  ,  and  the  corresponding 

b   X 

tensile  force  in  the  bar  is  AE  3_-  y   where  A  is  the  cross  sectional 
area  and  E  is  Young's  modulus. 

c.  A  simple  tension  in  the  x  direction  and  the  unit  elongation 

~  is  accompanied  by  lateral  contraction  of  the  amount  -v  -—  where  . 
dX  cxa. 

is  poisson's  ratio.  Inertia  forces  corresponding  to  lateral  particle 
motion  are  neglected  since  the  length  of  the  waves  is  large  in  comparison 
with  the  cross  sectional  dimensions  of  the  bar. 

The  foregoing  leads  to  the  equation  for  force  at  any  section 
in  the  bar  of  uniform  cross  section: 

F  =  EA  iL2  , 

ax 

The  net  force  on  a  section  of  bar  dX,   figure  (1),   is: 
*  V  -  Ife  .  li  -  Al(El||  •  £  (El|f)  dX)  -  I^U  -  Al5|  (El|£)  dX, 


m.,      rn.. 


db< 

■* ►> 


n,       az 


-/xb/r/  a^  />eg//?/7/s/<7  of  /s?A0c?cr 


X 


M/+/ 


~/fh/n. 


/^/00/ze  te) 
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assuming  that  E  may  not  be  constant  for  purposes  to  be  shown  below. 
Considering  the  motion  of  the  particle  dX,   and  applying  Newton's  second 
law  that  a  particle  acted  upon  by  an  unbalanced  force  system  has  an 
acceleration  in  line  with  and  directly  proportional  to  the  resultant 
of  the  force  system, 

A|-(Ei|)dX-AdX/  2»       simplifying 


2>X   *      *!'  J     >>t2 


(1) 


r; 


^t2       ^>X 
and  if  E  is  constant, 

c>  2D  _  JE_     £)2D 

at2        f     >>X2 

2         g 

Introducing  a  constant  c  of  a  value   such  that  c     =  .   this  reduces 

s 

to  the  familiar  form  of  the  wave  equation,  ref.  (lU): 

£>2D  =  c2  ,^2D 
^>t2"    hi2' 

Thus,  it  is  established  that  equation  (l)  is  a  form  of  the  wave  equa- 
tion. It  is  desired  to  modify  it  to  account  for  boundary  conditions 
expected  in  a  pile.  This  can  be  done  in  the  general  case  as  follows: 

with  R  representing  resistance  to  driving.  Equation  (2)  can  be  con- 
verted into  a  difference  equation  for  numerical  solution  by  standard 
methods  (ref.  2k) '• 
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j   r 

.  t2  ^D(x,t+At)  "  2D(x,t)  +  (x,t-At)J 


JjJ^V^H'Wx.t)  -  °(x,t)>  [K(x-^)][Dx^  "  VAx,t)]| 


KT 
where  E  ■  -j-  and  D(x, t)  is  the  displacement  of  point  x  during  time 

interval  t. 

It  is  convenient  at  this  point  to  introduce  a  simplified 
terminology  which  will  facilitate  comparison  with  model  equations  to 
be  developed  later,  and  to  correlate  with  existing  literature  on  this 
subject  (ref.  28). 

Let    D/ 


D(x,t+^t) 

— 

Dm 

D(x,t) 

= 

d 
m 

(x,t-  At) 

= 

d* 
m 

K(*+^) 

= 

K 
m 

K(x-^) 

= 

m-1 

(x+Ax,t) 

= 

d     4.T 

m+1 

D(x-Ax,t)  =  ^-1 
R         =  R. 


m 


Then 
? 


-£-*   [D  -  2d  +  d*l  =  — £-Jk  (d  _  -  d  )  -  K   (d  -  d  .  )  -  R  ( 

^ .  2  l  m    m    mJ   a  ay^  m  m  -*-    m    m*"-L  m    m~-J-    m-i 


Simplifying, 


Dm  "  2d»  "  «S  +  ^#f  dm-l  "  dA-l  "  K  -   WKm  "  Rm] 
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If  the  acceleration  of  gravity,  g,  is  expressed  in  feet  per  second, 
then 

\   -  2dm  "  <   +  *2*$£&«m*  -  dm'Km-l  "  K   -  W^m  "  ^  (3) 

m   u  J 

For  purposes  of  numerical  calculation,  it  is  more  convenient  to  use  a 
mathematical  model  consisting  of  weights  and  springs  such  as  is  shown 
in  figure  (2).  The  following  notation  will  be  used  and  is  consistent 
with  that  used  in  equation.(3)  above  and  in  reference  (28). 

1,  2,  3,  .  .  .  m-1,  m,  .  .  .  are  subscripts  designating  par- 
ticular weights  and  associated  springs, 

1,  2,  3,  .  .  .  n-1,  n,  .  .  .  are  subscripts  designating 
specific  time  intervals.  Zero  designates  the  initial  instant. 

The  following  apply  to  any  time  interval  n; 

Dm  =  displacement  of  weight  m  measured  from  its  initial 
position  in  inches 

C   =  compression  of  spring  m  in  inches 

Fm  »  force  exerted  by  spring  m,  in  pounds 

Z   =  net  force  acting  on  weight  m,  in  pounds 

V",  =  velocity  of  weight  m,  in  feet  Der  second 
m  J  °    ' 

The  following  apply  to  the  preceding  time  interval  n-1: 
cL  =  displacement  of  weight  m  measured  from  its  initial 

position  in  inches 
c   =  compression  of  spring  m,  in  inches 
f   =  force  exerted  by  spring  m,  in  pounds 
v   =  velocity  of  weight  m,  in  feet  per  second 
z   =  net  force  acting  on  weight  m,  in  pounds 


. 


. 


i 


' 


Ik 


The  following  are  usually  constants; 

Wm  =  magnitude  of  weight  m,  in  pounds 

K^  ■  spring  constant  for  spring  m,  in  pounds  per  inch 


R 
m 


=  external  force  or  resistance  acting  on  weight  m, 
in  pounds 


g   ■  acceleration  due  to  gravity  (32.2  feet/second) 

E   «■  Young's  modulus  of  elasticity  in  pounds/square  inch 

A   -  cross  sectional  area  of  pile  in  square  inches 

Q 

*"      -    unit  length  in  inches 

A  t  =  time  interval  used  for  numerical  calculations, 
in  seconds 

Tm  «•  critical  time  interval  for  spring  m,  in  seconds 

T  .  =  minimum  value  of  T  ,  in  seconds 
min  m' 

0   -  oror-all  "factor  of  safety"  ■  Tmin/A  t 

0m±n~     individual  factor  of  safety  for  spring  m  a  Tm/A.t 

The  numerical  calculations  will  be  a  step-by- step  process  in 

which  the  five  variables,  D,,  C  ,  F  ,  Z  ,  and  V  ,  will  be  calculated 

»  m'  m*  m5  m'     m' 

for  each  weight  or  spring  in  each  successive  time  interval.  Formulas 
for  these  five  variables  are  developed  by  Smith  (ref .  28)  as  follows: 

Dm  ■  dm  +  vm(12At)  (U) 

The  coefficient  12  is  included  since  D„,  and  d  are  expressed  in  inches 

mm 

and  v  in  feet  per  second. 

m        r 

Considering  C  next,  and  referring  to  figure  (3),  let  the  dashed 
squares  represent  the  initial  positions  of  weights  m  and  m  +  1,  and  let 
the  solid  squares  represent  their  positions  at  the  end  of  time  interval 
n.  The  initial  length  of  spring  m  is  il-,   and  at  the  end  of  time  interval  i 
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it  isX  .     D     and  Dm+-j_  are  the  displacements  of  weights  m  and  m  +  1 
respectively,   during  time  interval  n.     Then, 

but  A    *  Dm+1  mjf  -  Dm 

therefore,  Cm  =  Dm  -  Dm+1  (5) 

The  force  Fm  is  directly  proportional  to  Cm  and  to  the  spring 

constant  K  ,  so  that 

Fm  =  C  K  (6) 

m     m  m 

From  figure   (2),   weight  m  is  acted  upon  by  springs  m  -  1  and 

m  and  by  the  external  force  R     .     The  accelerating  force  is  given  by 

Z     «=  F     n   -  F     -  R  (7) 

m         m-1         mm  v '  ' 

From  Newton's  second  law  that  change  in  velocity  is  equal  to 

force  times  time  divided  by  mass, 

V  =  v  +  Z  (4£S)  (8) 

m    m    m  Wm  v  ' 

It  remains  to  be  shown  that  the  model  of  figure  (2)  and  equa- 
tions (h)   through  (8)  are  equivalent  to  the  wave  equation.  This  is 
done  by  using  equations  (l)  through  (8)  in  the  following  manner: 

Rewriting  equation  (8)  for  the  previous  time  interval  gives 

vm  "  *£  +  Z™(^f  >  (8a) 

Similarly,  equation  (7)  for  the  previous  time  interval  is 

*m  =  Vl  "  fm  "  *m  (7a) 


■ 
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Substituting  equation     (7a)  into  equation  (8a)  gives 

v™  "  vm"  +  (fm  l   -  *»  -  RJ  ^r^  (9) 

m        m       v  m-1        m        m'     w_  K7 ' 

Substituting  equation  (?)  into  equation  (U)  gives 

^  =  dm  +  [4  +  (Vl  "  fm  "  Bm)  ^]  12  At 

m 

But  equation  (U)  for  the  previous  time  interval  is 

^.dg^(12M) 
d-d* 

Substituting  (Ub)  into  (iia)  yields 

m 
Substituting  equations  (£.)  and  (6)  for  the  previous  time  interval 

into  the  above  equation  and  simplifying  gives 

Dm  "  2dm  "  <%  +  iTiJ  [(dm-1  "  ^^1  "  (*m  -  W^  -  Rj 

m 

which  is  equation  (3).  Thus  equations  (U)  through  (8)  are  equivalent 
to  the  wave  equation,  and  the  mathematical  model  of  figure  (2)  has 
been  shown  to  be  a  valid  one. 

Of  course,  solving  the  partial  differential  equation  by  a 
numerical  integration  process  is  an  approximate  one,  but  the  results 
should  be  accurate  within  five  per  cent  (ref .  28)  provided  that  the 
time  interval  chosen  is  small  enough  and  is  coordinated  with  the  length 
of  pile  selected  to  be  represented  as  a  block  and  spring  in  the  mathe- 
matical model  so  that  a  stable  set  of  calculations  is  produced.  This 
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matter  will  be  discussed  in  more  detail  later,  but  experience  indicates 
that  5- to  10-foot  sections  of  a  pile  may  be  so  represented  with  a  time 
interval  of  from  0.00025  to  0.00033  seconds,  depending  upon  the  type 
of  pile  material.  As  to  the  numerical  method,  it  may  be  considered 
sufficiently  accurate  in  light  of  present  knowledge  upon  which  the 
many  factors  influencing  any  pile  driving  operation  are  evaluated. 

B.  Mathematical  Model  Representation 

That  a  pile  may  be  treated  mathematically  as  a  series  of 
weights  and  springs  not  only  facilitates  the  computations  but  also 
assists  in  visualizing  the  problem.  For  that  reason,  the  remainder  of 
this  paper  will  deal  principally  in  terms  of  the  mathematical  model  of 
a  pile.  It  is,  therefore,  important  that  this  method  of  pile  repre- 
sentation be  thoroughly  understood  at  the  outset.  The  following  sec- 
tions will  adhere  closely  to  reference  (28). 

The  Hammer-Ram  is  usually  a  short,  heavy,  rigid  object  which 
may  be  represented  by  an  individual  weight  without  elasticity.  In 
figure  (U)  the  first  weight  W^  represents  the  ram. 

The  velocity  of  the  ram  at  the  instant  of  impact  is  needed  to 

start  the  numerical  calculations.  Ordinarily  the  rated  foot  pounds  of 

energy  of  the  hammer  is  given  by  the  manufacturer.  The  efficiency  is 

sometimes  given,  and  sometimes  must  be  assumed.  From  these  data  the 

velocity  may  be  calculated  at  impact  by  the  following  formula: 

velocity  at  impact    plated  energy  (ft-lb)  x  efficiency  x  6U7u    /-ns 
(ft/sec)     "  V  weight  (lb)  uu; 
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The  Capblock  is  a  short,  springy  object  of  wood,  plastic,  or 
similar  material  which  is  comparatively  light,  and  which  may  there- 
fore be  represented  by  a  spring.  In  figure  (U)  this  spring  is 
identified  as  K, .  The  form  of  the  stress-strain  diagram  (or  the 
hysteresis  loop)  that  is  produced  as  the  capblock  is  suddenly  com- 
pressed and  then  allowed  to  re-expand  is  assumed  to  be  as  shown  by 
figure  (5).  Compression  occurs  along  the  line  AB  whose  slope  is 
determined  by  the  elastic  constant  K,  of  the  capblock   Restitution 
occurs  first  along  the  line  BD,  and  then,  because  the  capblock  cannot 
transmit  tension,  it  is  completed  along  line  DA.  ABDA  is  the  hyster- 
esis loop.  The  computations  are  made  so  that 

Area  BCD     /   v2     Energy  output 
Area  ABC     ^el'     Energy  input 

where  e,  is  the  coefficient  or  restitution  of  the  capblock.  This 

relationship  may  be  derived  in  the  following  manner: 

Let  W    be  the  weight  of  the  ram 
ram  ° 

Wcap  be  the  weight  of  the  capblock 

vr    be  the  velocity  of  the  ram  before  impact 

vcap  be  the  velocity  of  the  caP  before  impact 

v1    be  the  velocity  of  the  ram  after  impact 
ram 

v1    be  the  velocity  of  the  cap  after  impact 
cap 

Then,  from  Newton's  laws, 

Wram  vram  +  Wcap  vcap  =  Wram  vram  +  Wcap  vcap 


and  by  definition, 

e  =   cap  ram 
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If     v     ■  v1    =  0 
Xi      cap     ycap      > 

then    6l   =  iSE 
J-    v 

ram 

2  ? 

energy  output   ( 1/2 )  (Wram/g )  ( vram)    (el  vram)     2 

energy  input  -  ^^^^2   "  (W*   "  * 

From  figure  (5), 

area  BCD    2  =  energy  output  ,   » 

area  ABC    1   "  energy  input  ^ 

During  compression  the  spring  force  may  be  computed  from  equation 

(6), 

Fj  -  CA  (12) 

Inelasticity  of  the  capblock  must  be  considered  during  restitution, 
however.     From  equation  (11)  and  figure  (£), 

e2  «  (1/2)  BCx  CD  .  CD 
1       (1/2)  BC  x  AC       AC 


but,  AC     =     C 


1  max 


and  CD     =     e2  C 

1     1  max 

by  similar  triangles  ADE  and  BCD, 

AE  AD 

BC  CD 


C-,  -  CD 

AE     .       lmaX •   BC 

CD 

_     _  1-4 

and  AE     -     -BC  (     m     ) 
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Writing  the  equation  for  line  DB  gives 

F   -  H  Cn  -  BC  (-^) 

x    CD  ±  ej 


and  substituting  for  CD, 


80     •c1-ar<i^) 


1   e2  C       X       e2 

1  1  max  l 


Sinoe  X  '  °1  max  Kl  • 


cl  max  ^1  _        M) 
1    e2c        imax  V-Z-> 

1  1  max 
which  reduces  to 

Fl  =  -±±  -  C.     £,("?-  1)  (13) 

12    1  max  i  e2 

Calculations  involving  the  capblock  will,  therefore,  utilize 
alternately  equation  (12)  during  compression  and  equation  (13)  during 
restitution.  Cn  _._„  should  not  be  thought  of  as  a  constant,  but  as 
the  maximum  compression  occurring  in  spring  K^  before  restitution  be- 
gins. It  often  happens  that  recompression  occurs  and  a  new  value  of 
C,  max  is  reached  before  restitution  begins  again.  This  is  illus- 
trated by  figure  ($A)  where  the  stress-strain  diagram  begins  at  A, 
goes  to  B  during  compression  following  equation  (12),  goes  to  C 
during  restitution  as  governed  by  equation  (13),  at  C  recompression 
begins,  goes  to  B  following  equation  (13),  continues  to  D  by  equa- 
tion (12),  and  thence  to  E  by  equation  (13).  As  previously  described 
the  return  to  A  is  along  EA  because  no  tension  can  be  transmitted. 
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That  recompression  can  occur  may  not  be  apparent  readily,  but  a 
feeling  for  it  may  be  acquired  in  solving  a  problem  or  two  manually. 

As  will  be  noted  later,  equations  (12)  and  (13)  are  used  in 
the  case  of  a  cushion  block,  but  with  the  subscripts  changed  to  2. 

Few  tests  have  been  made  to  determine  the  elastic  character- 
istics of  capblocks  under  impact  conditions.  Raymond  International, 
Inc.,  in  conducting  a  small  number  of  tests  found  that  the  character- 
istics of  a  wooden  capblock  vary  during  driving,  and  concluded  that 
in  order  to  be  on  the  conservative  side  in  computing  pile  penetration 
per  blow,  a  hardwood  capblock  with  the  grain  vertical  and  of  an 
original  six-inch  height,  may  be  assumed  to  have  the  following 
char ac  teristics: 

Spring  constant  K  =  20,000  A  (lb/inch  of  compression) 

where  A  is  the  horizontal  cross  sectional  area  in  square  inches.  The 
coefficient  of  restitution  was  0.5.  The  tests  also  showed  that  a 
Micarta  capblock  (Nema  Grade  "C" )  of  12  inches ■  height  had  a  spring 
constant  given  by  the  following  relation: 

K  =  US, 000  A 

and  the  coefficient  of  restitution  was  0.8. 

Pile  cap,  or  Follower,  or  Helmet— like  the  ram  of  the  hammer, 
the  pile  cap  is  ordinarily  a  short,  heavy,  rigid  object  that  can  be 
represented  by  a  single  weight  without  elasticity,  such  as  Wg  in 
figure  (U).  If  the  pile  cap  is  long  and  slender,  as  is  the  case  when 
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used  as  a  follower  to  drive  piles  below  water  or  the  ground,  it  is 
better  to  represent  it  by  a  number  of  weights  and  springs  as  shown 
in  figure  (6e).  In  this  case  elastic  constants  must  be  computed  by 
the  formula 

i  -  f  m 

where  A  is  the  cross  sectional  area  in  square     inches,  E  is  the 
modulus  of  elasticity  in  pounds  per  square  inch,  and  A-  is  the  unit 
length  in  inches  represented  by  one  spring. 

Cushion  Blocks — (called  "head  packing  in  the  U.K.) — in  figure 
(U)  springs  Y^  to  Kq  inclusive  represent  the  elasticity  of  the  pile 
itself.     However,   if  a  precast  concrete  pile  is  being  driven,   a  cushion 
block  must  be  used  under  the  pile  cap  V\L  to  protect  the  concrete  from 
shattering.     In  this  case,  spring  K2  would  represent  the  cushion  in 
combination  with  the  first  spring  of  the  pile.     Figure  (6c)  applies 
to  the  cushion  block  as  well  as  the  capblock,  and  they  are  treated 
mathematically  in  the  same  manner.     Dynamic  tests  similar  to  those 
described  for  capblocks  show  that  U-inch-thick  pine  boards  with  the 
grain  horizontal  as  used  on  top  of  a  precast  pile  to  distribute  the 
blow  evenly  may  be  assumed  to  have  a  spring  constant  equal  to  3,li80  A 
and  a  coefficient  of  restitution  of  0.5.     Of  course,  the  character- 
istics of  the  wood    vary    during  driving,  but  use  of  the  above 
relationship  is  on  the  conservative  side  when  used  for  calculating 
penetration  of  the  pile  hammer  blow. 

The  Pile  will  be  either  of  wood,  concrete,   or  steel.     All  of 
these  materials  have  traditionally  been  treated  structurally  as  behaving 
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elastically  within  the  working  load  ranges.  Because  the  pile  is 
compressible  due  to  its  length,  it  is  subject  to  wave  action  under  the 
blow  of  a  hammer.  This  wave  action  may  be  analyzed  by  dividing  the 
pile  into  short  unit  lengths  of  5  to  10  feet.  The  weight  of  each 
length  is  noted  as  W  through  W1Q,  and  the  elasticity  of  each  length 
by  an  individual  spring  as  shown  by  K2  through  K0  in  figure  (U).  If 
the  pile  is  of  uniform  section,  the  weights  and  springs  representing 
it  are  identical.  If  the  pile  is  tapered  or  is  a  composite  pile,  or 
is  in  any  way  of  irregular  cross  section,  then  the  weights  and  springs 
are  modified  to  represent  the  pile  under  investigation.  That  com- 
paratively large  unit  lengths,  such  as  10  feet,  give  desired  accuracy 
in  computing  the  action  of  impact  waves  may  not  be  clear.  Although 
water  waves  are  transverse  waves  and  the  waves  in  a  pile  are  longitudinal 
waves,  water  waves  nevertheless  can  be  used  as  an  analogy  to  illustrate 
the  principles  involved. 

If  a  long  strip  of  flexible  material  were  allowed  to  float  on 
the  surface  of  a  body  of  water,  it  would  follow  the  wave  action.  If, 
for  purposes  of  mathematical  analysis,  this  flexible  strip  were  rep- 
resented by  a  number  of  short  floats  of  rigid  construction  but  con- 
nected by  flexible  links,  the  mathematical  model  would  appear  as  in 
figure  (7).  This  representation  would  be  an  approximation,  but  it 
would  involve  negligible  error  because  the  small  floats  would  ride 
the  waves  almost  exactly  like  the  flexible  strip.  If,  however,  the 
rigid  floats  are  made  comparatively  long  and  they  approach  or  exceed 
the  length  of  the  water  wave,  then  the  model  would  be  as  shown  by 
figure  (8).  Obviously  these  long  floats  cannot  follow  the  wave  form 
closely. 
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Applying  this  analogy  to  a  pile  being  hit  by  a  ram,  the 
importance  of  dividing  the  pile  into  unit  lengths  considerably  shorter 
than  the  stress  wave  produced  by  a  hammer  blow  can  be  appreciated. 
Fortunately,  a  pile  driving  impact  usually  produces  a  fairly  long 
wave  form,  and  unit  lengths  of  the  order  of  5  to  10  feet  will  produce 
acceptable  accuracy.  In  the  special  case  where  the  exact  form  of 
the  impact  wave  is  being  investigated,  a  smaller  unit  length  of  1  or 
2  feet  may  be  advisable.  Under  these  conditions,  it  may  also  be  ad- 
visable to  divide  the  ram  into  a  series  of  weights  and  springs  as  is 
done  with  the  pile.  For  normal  pile  driving  analysis,  such  refinements 
are  not  considered  justified* 

There  is  a  very  important  relationship  between  the  unit  length 
of  pile  selected  and  the  time  interval  chosen  for  the  calculations, 
the  time  interval  to  be  used  in  the  calculations  is  a  function  of  the 
weight  of  the  unit  length  selected  and  of  the  elasticity  and  density 
of  the  pile  material.  This  relationship  is  derived  and  presented 
later,  and  governs  the  upper  bound  of  the  time  interval  to  be  used. 
On  the  other  hand,  choice  of  too  small  a  time  interval  results  in 
excessive  computation  for  the  accuracy  required. 

The  spring  constants  are  determined  by  the  formula  given  for 
K  under  "Capblock"  above. 

The  Ground — as  indicated  by  figure  (k)   provision  is  made  for 
ground  resistance  to  be  applied  on  each  section  of  the  pile.  In 
order  to  make  a  calculation  using  the  wave  equation,  an  assumption 
must  be  made  as  to  the  amount  of  resistance  which  will  be  offered  by 
the  ground  at  each  of  the  pile  sections.  Furthermore,  a  stress-strain 
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relationship  of  the  soil  and  viscous  damping  must  be  considered.  For 
the  purpose  of  this  section,  it  will  be  assumed  that  these  factors 
can  be  evaluated  so  that  the  method  of  their  inclusion  in  the  analysis 
may  be  developed. 

Resistance  at  the  pile  point— according  to  Chellis  (ref.  5) 
the  ground  compresses  elastically  for  a  certain  distance  and  then 
fails  plastically  with  constant  ultimate  resistance  R^.     This  concept 
is  illustrated  by  figure  (9).  The  elastic  deformation  will  be  called 
"quake,"  and  is  denoted  by  the  letter  'Q.'  The  plastic  deformation  will 
be  referred  to  as  "permanent  set,"  and  is  identified  as  's'$  in  figure 
(9)  it  is  the  distance  AB  or  OC.  There  is  a  factor  of  importance 
which  has  not  yet  been  included,  and  this  is  the  rate  of  stress 
application.  It  is  clear  that  the  ground  will  offer  more  instanta- 
neous resistance  to  a  rapid  advancement  of  the  pile  tip  than  to  a  slow 
advancement.  Viscous  damping,  therefore,  will  be  included  in  the 
analysis  to  provide  for  the  effect  of  rate  of  penetration. 

As  has  been  shown,  the  wave  equation  calculation  will  give  the 
instantaneous  velocity  of  the  point  of  the  pile  in  any  time  interval. 
If  J  is  designated  as  a  damping  constant  and  v  the  velocity  of 
the  pile  point,  then  their  product,  Jvp,  can  be  used  to  increase  or 
decrease  ground  resistance  at  the  point  of  the  pile  so  as  to  produce 
damping.  The  instantaneous  damping  resistance  at  point  x,  for  example, 
in  figure  (9)  is  Jv X  .  This  value  would  be  used  for  %o  for  a  P31*" 
ticular  time  interval  applied  as  shown  by  figure  (U).  This  damping 
resistance  is  a  resistance  to  driving  the  pile  and  does  not,  of  course, 
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contribute  to  the  bearing  capacity  of  the  pile  in  supporting  static 
loads . 

Resistance  Along  the  Sides  of  the  Pile — side  resistance  is 
calculated  in  the  same  manner  as  at  the  pile  point,  except  that  a 
different  damping  constant  J1  is  used,  and  it  applies  to  R^  through 
R  in  figure  (U).  As  the  pile  is  driven,  the  soil  under  the  pile 
point  is  displaced  and  caused  to  flow  aside  very  rapidly.  On  the 
other  hand,  the  soil  at  the  sides  of  the  pile,  as  it  is  being  dirven, 
is  not  displaced  nearly  as  much.  It  may  be  assumed,  then,  that  J1 
should  be  less  than  J. 

Distribution  of  Ground  Resistance  Between  the  Pile  Sides  and 
Point — the  wave  equation  computation  will  utilize  any  distribution 
chosen.  The  problem  of  what  distribution  to  choose  will  be  discussed 
at  length  later. 

Figure  (10)  shows  the  model  representation  of  the  ground 
resistance  along  the  sides  of  the  pile.  It  must  be  remembered  that 
ground  quake  is  an  elastic,,  or  recoverable,  deformation  and  that  the 
set  associated  with  it  is  of  a  temporary  nature.  The  deflection  of 
the  spring  K^  must,  therefore,  be  limited  to  a  maximum  value  less  than 
or  equal  to  Q,  the  ground  quake.  As  the  maximum  allowed  deflection 
of  spring  K*  is  reached,  provision  is  made  for  plastic  deformation 
of  the  ground,  D*  in  the  normal  downward  direction,  or  in  an  upward 
direction  in  the  case  of  pile  rebound.  Figure  (11 )  illustrates  the 
way  in  which  D^  is  handled  in  the  calculations.  It  may  be  noted  that 
D^  does  not  take  on  a  value  until  1^  exceeds  Q.  From  that  point  on, 
Dm  cannot  be  less  than  Dm  -  Q  or  more  than  Dm  +  Q. 
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The  value  of  the  ground  resistance  on  a  particular  pile  sec- 
tion is  computed  by  the  equation 

The  point  of  the  pile  is  dealt  with  in  a  similar,  but  slightly- 
different,  manner.  Only  downward  plastic  ground  movement  Dl  is  con- 
sidered. The  maximum  value  of  D*  is  the  permanent  set  of  the  pile 
which  is  designated  as  s  .  Figure  (12)  illustrates  the  way  in 
which  the  calculations  are  handled.  The  value  of  the  ground  re- 
sistance at  the  point  is  computed  by  the  equation 

Rp  -  (Dp  -  D£)  K£  (1  +  Jvp)  (16) 

C.  Illustrative  Problem 

In  order  to  bring  the  theory  more  sharply  into  focus,  a 
hypothetical  problem  will  be  set  up  and  solved  manually  using  the 
general  format  described  by  E.  A.  L.  Smith  (ref.  27). 

Let  figure  (13)  represent  a  steel  pile  30  feet  long.  The 
weight  of  the  ram  is  £,000  pounds,  and  its  velocity  at  impact  as 
determined  by  equation  (10)  is  found  to  be  10  feet  per  second.  The 
capblock  is  hardwood,  has  a  coefficient  of  restitution  of  0.5,  and 
its  spring  constant  as  determined  by  equation  (5)  is  2,000,000  pounds 
per  inch.  The  pile  cap  weighs  500  pounds,  aid  is  not  able  to  transmit 
tension.  No  cushion  block  or  head  packing  will  be  used,  so  K?  will 
be  the  elasticity  of  the  first  pile  section  whose  coefficient  of 
restitution  will  be  taken  as  1.00.  The  pile  will  be  divided  into 
10-foot  unit  lengths,  each  of  which  weighs  500  pounds.  The  spring 
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constants  Kg,  Kj,  and  K^  for  the  pile  as  computed  by  equation  (6) 
are  each  taken  to  be  k, 000, 000  pounds  per  inch.  Ground  quake  will  be 
0.10  inch.  The  ultimate  ground  resistance  has  been  determined  to  be 
220,000  pounds.  It  is  distributed  as  follows:  10,000  pounds  on  each 
of  the  two  top  sections  of  the  pile  and  200,000  pounds  on  the  bottom 
10- foot  section  of  pile.  Dividing  the  ultimate  ground  resistance  by- 
ground  quake  will  give  the  soil  spring  constants  as  follows: 
K'  ■  k£  =  100,000  pounds  per  inch 

and  Kl     «  2,000,000  pounds  per  inch 

The  time  interval  for  the  calculations  will  be  l/IiOOO  seconds.  The 
damping  constants  J  and  J'  will  both  be  taken  as  equal  to  0.1  for 
purposes  of  this  example,  although,  as  noted  before,  J'  may  be  smaller 
than  J. 

The  calculations  may  then  be  performed  by  slide  rule  and 
recorded  in  tabular  form  as  shown  on  the  following  pages.  Although 
space  does  not  permit  here,  it  is  recommended  that  when  making  such  a 
table  a  single  large  sheet  be  used.  In  referring  to  Table  A^  on  each 
line  values  are  first  computed  for  D.. ,  Dp,  D»,  D,  ,  D^j  then  C   C_, 
C^,  C]^  Ccj;  etc.,  until  all  computations  have  been  completed.  The 
next  line  is  approached  in  the  same  order,  and  so  on. 

It  may  be  of  interest  to  note  from  Table  A  that  C^  reaches  a 
first  maximum  at  time  interval  8j  restitution  takes  place  from  time 
interval  9  through  22.  Recompression  occurs  from  time  interval  23 
through  27  at  which  time  a  second  C,     is  reached.  Restitution 
continues  until  the  end  of  the  calculations. 
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TABLE  A 

Fl- 

hm 

Cn      ■ 

(2xl06)(C1), 

Zl  ■ 

z 

Time 

^+0.003  n 

0 

DX  -D2 
0 

except  during 
restitution 

0 

-Fl 

1 

Interval 

vl+  621,000 

0 

0 

10.00 

1. 

0.03 

0.03 

60, 000 

-  60,000 

9.9035 

2. 

0.06 

0.057 

111*,  000 

-111*,  000 

9.7202 

3. 

0.0892 

0.0781* 

156,800 

-156, 800 

9.U68 

u. 

0.1176 

0.0931* 

186,800 

-186,800 

9.167 

5. 

0.1U51 

0.1025 

205,000 

-205,000 

8.837 

6. 

0.1716 

0.1073 

211*,  600 

-211*,  600 

8.1*91 

7. 

0.1971 

0.1092 

218,1*00 

-218,1*00 

8.139 

8. 

0.2215 

0.1096 

219,200 

-219,200 

7.786 

9. 

0.21*1*9 

0.1091 

215,200 

-215,200 

7.1*1*0 

10. 

0.2672 

0.1081 

207,200 

-207,200 

7.106 

11. 

0.2885 

0.1069 

197,600 

-197,600 

6.788 

12. 

0.3088 

0.1058 

188,800 

-188,800 

6.1*81* 

13. 

0.3283 

0.1053 

181*,  800 

-181*,  800 

6.186 

lit. 

0.31*69 

0.1051 

183,200 

-183,200 

5.891 

15. 

0.361*6 

0.101*9 

181,600 

-181,600 

5.599 

16. 

0.3811* 

0.101*3 

176,800 

-176, 800 

5.3H* 

17. 

0.3973 

0.1033U 

169,120 

-169,120 

5.01*2 

18. 

0.1*121* 

0.10221* 

160,320 

-160,320 

U.781* 

19. 

0.1*267 

0. 10161* 

155,520 

-155,520 

U.531* 

20. 

0.1*1*03 

0.10238 

161,1*1*0 

-161,1*1*0 

U.271* 

21. 

0.1*5312 

0.101*79 

180,720 

-180.720 

3.981* 

22. 

0.1*6507 

0.10855 

210, 800 

-210,800 

3.61*1* 

23. 

0.1*7600 

0.11266 

225,320 

-225,320 

3.281 

21*. 

0.1*8581* 

0.11687 

233,71*0 

-233,71*0 

2.905 

25. 

0.1*91*56 

0.12095 

21*1, 900 

-21*1,900 

2.516 

26. 

0.50211 

0.121*3 

21*8,600 

-21*8,600 

2.116 

27. 

0.508U6 

0.12612 

252,21*0 

-252,21*0 

1.710 

28. 

0.51376 

0.12580 

21*9,680 

-21*9,680 

1.309 

29. 

0.51769 

0.12291* 

226, 800 

-226,800 

0.91*1* 

30. 

0.52052 

0.11925 

197,280 

-197,280 

0.626 

31. 

0.5221*0 

0.11681* 

178,000 

-178,000 

0.31*0 

32. 

0.5231*2 

0.11697 

179,01*0 

-179,01*0 

0.052 

33. 

0.52357 

0.11913 

196,320 

-196,320 

-0.261* 

3k. 

0.52278 

0.1211*2 

211*,  61*0 

-211*,  61*0 

-0.610 

35. 

0.52095 

0.12185 

218,080 

-218,080 

-0.962 

36. 

0.51806 

0.11971 

200, 960 

-200,960 

-1.28U 

37. 

0.5ll*2l 

0.11586 

170,160 

-170,160 

-1.558 

38. 

0.50951* 

0.1121*7 

11*3,01*0 

-11*3,01*0 

-1.789 

39. 

0.501*17 

0.11129 

133,600 

-133,600 

-2.001* 

TABLE  A  (cont'd.) 
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D2  " 

C2  = 

F2  = 

Z2  = 

v2  - 

Time 

d2+  .003  v2 

D2  -D3 

too6c2 

»l-»2 

2 

Interval 

V  62,100 

1 

0 

0 

0 

60',  000 

0.965 

2 

0.003 

0.003 

12,000 

102,000 

2.605 

3 

0.01081 

0.0102 

l*o, 800 

116,000 

1*.1*75 

h 

0.021*23 

0.02123 

81*,  920 

101,880 

6.115 

S 

.01*257 

0.03357 

131*,  280 

70,720 

7.255 

6 

.061*31* 

o.ol*l*5 

178,000 

36,600 

7.81*5 

7 

.08788 

0.05208 

208,320 

10,080 

8.007 

8 

.11190 

0.0557 

222,800 

-3,600 

7.9U9 

9 

.13575 

0.05655 

226,200 

-11,000 

7.772 

10 

.15905 

0.05595 

223,800 

-16,600 

7.505 

11 

.18157 

0.05517 

220,680 

-23,080 

7.13U 

12 

.20297 

o.o5la7 

216,680 

-27,880 

6.681* 

13 

0.22302 

0.05262 

210, 1*80 

-25,680 

6.270 

ll* 

.21*183 

0.05053 

202,120 

-18,920 

5.965 

15 

0.25972 

0.01*825 

193,000 

-11,1*00 

5.781 

16 

0.27706 

0.01*65 

186,000 

-  9,200 

5.633 

17 

0.29396 

0.01*59 

183,600 

-ll*,l*80 

5.U00 

18 

0.31016 

0.01*68 

187,200 

-26, 880 

U.967 

19 

0.32506 

0.01*91*8 

197,920 

4*2,1*00 

U.285 

20 

0.33792 

0.05303 

212,120 

-50,680 

3.1*69 

21 

0.31*833 

0.05667 

226,680 

-1*5,960 

2.729 

22 

0.35652 

0.05980 

239,200 

-28,1*00 

2.272 

23 

0.3633U 

0.0621*1* 

21*9, 760 

-2U,UU0 

1.878 

21* 

0.36897 

0.06359 

251*,  360 

-20,620 

1.51*6 

25 

0.37361 

0.06275 

251,000 

-  9,100 

1.1*00 

26 

0.37781 

0.060l*6 

21*1,81*0 

+  6,760 

1.509 

27 

0.3823U 

0.057U1 

229,61*0 

+22,600 

1.873 

28 

0.38796 

0.05635 

225, 1*00 

+21;,  280 

2.261* 

29 

0.391*75 

0.05808 

232,320 

-  5,520 

2.175 

30 

0.1*0127 

0.06092 

21*3,680 

-1*6,1*00 

1.1*29 

31 

0.1*0556 

0.06208 

21*8,320 

-70,320 

0.298 

32 

0.1*061*5 

0.05977 

239,080 

-60, 0U0 

-0.669 

33 

0.1*01*1*1* 

0.051*66 

218,61*0 

-22,320 

-1.028 

3U 

0.1*0136 

0.01*91*1 

197, 61+0 

+17,000 

-0.75U 

35 

0.39910 

0.01*671 

186,81*0 

+31,21*0 

-0.251 

36 

0.39835 

0.01*636 

185, 1*1*0 

+15,520 

-0.001 

37 

0.39835 

0.01*917 

196,680 

-26,520 

-0.1*28 

38 

0.39707 

0.0508 

203,200 

-60,160 

-1.398 

39 

0.39288 

0.01*769 

190, 760 

-57,160 

-2.320 

TABLE  A  (cont'd.) 


29b 


D. 


2+0.003  v3 


1 

0 

2 

0 

3 

0.0006 

:tU 

0.0030 

5 

0.0090; 

6 

0.0198 

7 

0.0358 

8 

0.0562 

9 

0.0792 

10 

0.1031 

11 

0.126U 

12 

O.H188 

13 

0.170U 

1U   - 

0.1913 

15 

0.2121*7 

16 

0.23056 

17 

0.2U800 

18 

0.26332 

19 

0.27558 

20 

0.28U89 

21 

0.29166 

22 

0.29672 

23 

0.30090 

2U 

0.30538 

25 

0.31086 

26 

0.31735 

27 

0.32U93 

28 

0.33161 

29 

0.33667 

30 

0.3U035 

31 

0.3U3U8 

32 

0.3U668 

33 

0.3U978 

3U 

0.35195 

35 

0.35239 

36 

0.35200 

37 

0.3U918 

38 

0.3U627 

39 

0.3U519 

C3  = 


0 
0 

0.0006 

0.0029 

0.00823 

0.01693 

0.02809 

0.0396 

O.OU883 

0.05U2U 
0.055U7 
0.05U02 

0.05175 

0.05003 

O.OU965 

0.05095 

0.0539 

0.05791 

0.06122 

0.06289 
0.062U9 
0.06151 
0.05803 
0.05557 
0.05U55 
0.05181 

0.05893 
0.06166 
0.06228 
0.06096 
0.05898 

0.05751 
0.05670 
0.05571 
0.05363 
0.05139 
0.0U7UU 

0.QUU19 
0.0UU66 


1ix106C. 


0 

0 

2,U00 

11,600 

32,920 

67,720 

112,360 

158, Uoo 

195,320 
216,960 

221,880 
216,080 
207,000 
200,120 
198,600 
203,800 
215,600 
231, 6U0 
2UU, 880 
251,560 
2U9,960 
2U6, OUO 
232,120 
222,280 
218,200 
207,2UO 
235,720 
2U6,6U0 

2U9,120 
2U3,8UO 

235,920 

230, OUO 

226,800 

222, 8U0 

2lU,520 

205,560 

189,760 

156, 760 

178, 6U0 


52.1 


Q 


0.0031 

0.026U 

O.OU88 

0.070U 

0.0913 

0.1111*7 

0.13056 

0.1U800 

0.16332 

0.17558 

0.18U89 

0.19166 

O.19672 

0.20090 

0.20538 

0.21086 

0.21735 

0. 221*93 

0.23161 

0.23667 

0.21*035 

0.2U3U8 

0.21*668 

0.21*978 

0.25195 

0.25239 

0.25200 

0.2U918 

0.21*627 

0.21*519 
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inter- 
val 

D3 

D3"D3 

K^(D3-Dp 

1  +  J'v3 

1 
2 

0 

0 

0 

0 

3 

0 

0.0006 

60 

1.0193 

1* 

0 

.0030 

300 

1.0809 

5 

0 

0.0090 

900 

1.1985 

6 

0 

0.01980 

1,980 

1.3601 

7 

0 

0.03580 

3,580 

1.5336 

8 

0 

0.0562 

5,620 

1.6795 

9 

0 

0.0792 

7,920 

1.7680 

10 

0.0031 

0.10 

10,000 

1.7952 

11 

0.0261* 

0.10 

10,000 

1.773 

12 

0.01*88 

0.10 

10,000 

1.71*58 

13 

0.070U 

0.10 

10,000 

1.7186 

111 

0.0913 

0.10 

10,000 

1.6965 

15 

0.1111*7 

0.10 

10,000 

1.6721* 

16 

0.13056 

0.10 

10,000 

I.636U 

17 

0.11*800 

0.10 

10, 000 

1.5811* 

18 

0.16332 

0.10 

10,000 

1.50UU 

19 

0.17558 

0.10 

10,000 

1.1*087 

20 

0.181*89 

0.10 

10,000 

1.3103 

21 

0.19166 

0.10 

10,000 

1.2257 

22 

0.19672 

0.10 

10,000 

1.1686 

23 

0.20090 

0.10 

10,000 

1.1391* 

2k 

0.20538 

0.10 

10,000 

l.ll*9l* 

25 

0.21086 

0.10 

10, 000 

1.1826 

26 

0.21735 

0.10 

10,000 

1.2163 

27 

0.221*93 

0.10 

10,000 

1.2525 

28 

0.23161 

0.10 

10, 000 

1.2225 

29 

0.23667 

0.10 

10, 000 

1.1686 

30 

0.21*035 

0.10 

10,000 

1.1227 

31 

0.21*31*8 

0.10 

10,000 

1.101*1* 

32 

0.21*668 

0.10 

10,000 

1.1066 

33 

0.21*978 

0.10 

10,000 

1.1033 

3U 

0.25195 

0.10 

10,  000 

1.0721* 

35 

0.25239 

0.10 

10,000 

1.011*6 

36 

0.25239 

0.09961 

9,961 

0.9536 

37 

0.25239 

0.09679 

9,679 

0.9059 

38 

0.25239 

0.09388 

9,388 

0.9029 

39 

0.25239 

0.09280 

9,280 

0.961*1 

Ki(D3^)-    Z3 


0 

0 

60 

321* 

1,079 

2,690 

5,1*90 

9,1*30 

11*,  000 

17,952 

17,730 

17,1*58 

17,186 

16,965 
16,721* 
16,361* 
15,311* 
15,01*1* 
11*,  087 
13,103 
12,257 
11,686 
11,391* 
11,1*9U 
11,826 
12,163 
12,525 
12,225 
11,686 
11,227 
11,01*1* 
11,066 
11,033 
10,721* 
10,11*6 
9,1*90 
8,760 
8,1*1*0 
8,91*2 


rp-r  ^""tto 

3    ^2.100 

0 

0 

12,000 

0.1931 

38,31*0 

0.8091 

72,996 

1.985 

100, 281 

3.601 

107,590 

5.336 

90,1*70 

6.795 

51*,  970 

7.680 

16,880 

7.952 

-11,112 

7.773 

-18,930 

7.1*58 

-16,858 

7.186 

-13,706 

6.965 

-11*,  965 

6.721* 

-22,321* 

6.361* 

-31*,  161* 

5.811* 

-1*7, 8ll* 

5.01*1* 

-59,1*81* 

U.087 

-61,01*7 

3.103 

-52,51*3 

2.257 

-35,537 

1.686 

-18,526 

1.391* 

+  6,21*6 

1.1*91* 

+20,586 

1.826 

+20,971* 

2.163 

+22,1*37 

2.525 

-18,605 

2.225 

-33,1*65 

1.686 

-28,1*86 

1.227 

-11,387 

1.01*1* 

+  1,356 

1.066 

-  2,026 

1.033 

-19,193 

0.721* 

-35,921* 

0.11*6 

-37,826 

-0.1*61* 

-29,610 

-0.9U1 

-  1,81*0 

-0.971 

+38,000 

-0.359 

+  3,178 

-0.308 

Di 
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Time 
Interval 

di +0.003  v. 

Dl*-D5 

l*xl06C, 

DU-Q 

1 
2 
3 

0 

0 

0 

k 

0.00012 

0.00012 

1*80 

«• 

5 

0.00077 

0.00075 

3,000 

- 

6 

0.00287 

0.00269 

10, 760 

. 

7 

0.00771 

0.0069 

27,600 

- 

8 

.01660 

0.01382 

55,280 

- 

9 

.03037 

0.02329 

93,160 

. 

10 

0.01*886 

0.03376 

135,01*0 

- 

11 

0.07093 

0.01*313 

172,520 

- 

12 

0.091*79 

0.01*978 

199,120 

- 

13 

0.11865 

0.05365 

211*,  600 

0.01865 

11* 

0.11*127 

0.05636 

225,1*1*0 

0.01*127 

15 

0.16182 

0.05975 

239, 000 

0.06182 

16 

0.17961 

0.061*02 

256,080 

0.07961 

17 

0.191*10 

0.06663 

266,520 

0.091*10 

18 

0.205U1 

0.06668 

266,720 

0.105U1 

19 

0.211*36 

0.061*77 

259,080 

0.111*36 

20 

.22200 

0.06219 

21*8, 760 

0.12200 

21 

0.22917 

O.O6OO8 

21*0,320 

0.12917 

22 

0.23521 

0.05518 

220,720 

0.13521 

23 

0.21*287 

0.05865 

23U, 600 

0.11*287 

21* 

0. 21*981 

0.051*75 

219,000 

0.11*981 

25 

0.25631 

0.05603 

221*,  120 

0.15631 

26 

0.26191* 

0.05696 

227,81*0 

0.16191* 

27 

0.26600 

0.05670 

226, 800 

0.16600 

28 

0.26995 

0.05627 

225,080 

0.16995 

29 

0.271*39 

0.05651* 

226,160 

0.171*39 

30 

0.27939 

0.0571*5 

229,800 

0.17939 

31 

0.281*50 

0.05831* 

233,360 

0.181*50 

32 

0.28919 

0.0585U 

23U, 160 

0.18917 

33 

0.29308 

0.05777 

231,080 

0.19308 

3U 

0.29621* 

0.05626 

225,  OUO 

0.19621* 

35 

0.29876 

0.051*1*0 

217,600 

0.19876 

36 

0.30061 

0.0521*3 

209,720 

0.20061 

37 

0.30171* 

0.01*898 

195,920 

0.20171* 

38 

0.30208 

0.01*61*1 

185,61*0 

0.20208 

39 

0.30053 

0.01*358 

17U,320 

0.20053 

TABLE  A  (cont'd.) 
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Ri 


V, 


%MJ-      z 


i.  V!k  &V&!  x  +  J'vu  (1  +  J  V    F3-VRi*  Vsriro 


U  z^ 


1 

2 

3 
1* 

0 

0.00012 

12 

5 

0 

0.00077 

77 

6 

0 

0.00287 

287 

7 

0 

0.00771 

771 

8 

0 

0.01660 

1,660 

9 

0 

0.03037 

3,037 

10 

0 

0.01*886 

U,  886 

11 

0 

0.07093 

7,093 

12 

0 

0.09U79 

9,1*79 

13 

0.01865 

0.10 

10, 000 

11* 

0.01*127 

0.10 

10, 000 

15 

0.06182 

0.10 

10,000 

16 

0.07961 

0.10 

10,000 

17 

0.091*10 

0.10 

10,000 

18 

0.105U1 

0.10 

10,000 

19 

0.111*36 

0.10 

10,000 

20 

0.12200 

0.10 

10, 000 

21 

0.12917 

0.10 

10, 000 

22 

0.13521 

0.10 

10, 000 

23 

0.11*287 

0.10 

10,000 

2k 

0.1U981 

0.10 

10,000 

25 

0.15631 

0.10 

10,000 

26 

0.16191* 

0.10 

10,000 

27 

0.16600 

0.10 

10,000 

28 

0.16995 

0.10 

10, 000 

29 

0.171*39 

0.10 

10, 000 

30 

0.17939 

0.10 

10,000 

31 

0.181*5  0 

0.10 

10, 000 

32 

0.18917 

0.10 

10,000 

33 

0.19308 

0.10 

10, 000 

31* 

O.1962I* 

0.10 

10, 000 

35 

O.19876 

0.10 

10,000 

36 

0.20061 

0.10 

10, 000 

37 

0.20171* 

0.10 

10,000 

38 

0.20208 

0.10 

10, 000 

39 

0.20208 

0.0981*5 

9,81*5 

1.003865 

1.0218 

1.0699 

1.1613 

1.2963 

1.U59 

1.6163 

1.7356 

1.7953 

1.7952 

1.75U1 

1.6851 

1.5929 

1.1*830 

1.3770 

1.298U 

1.25U7 

1.2390 

1.23U6 

1.2555 

1.2313 

1.21678 

1.18768 

1.1351*8 

1.13156 

1.1U806 

1.16656 

1.17038 

1.15568 

1.1301*8 

1.10538 

1.081*03 

1.0616 

1.0378 

1.0112 

0.91*85 


12.05 
78.75 
307 
895 

2,11*5 

1*,1*20 

7,900 
12,300 
17,000 
17,952 
17,5U1 
16,851 
15,929 
H*, 830 
13,770 
12,981* 
12,51*7 
12,390 
12,31*6 
12,555 
12,313 
12,168 
11,877 
11,355 
11,316 
11,1*81 
11,666 
11,701* 
11,557 
11,305 
11,051* 
10, 81*0 
10,616 
10,378 
10,112 

9,31*0 


2,1*00 

11,108 

29,81*1 

56,653 

83,865 

100,975 

97,71*0 

71*,  020 

37,060 

-1*0 

-2Sy$S0 

-1*2,861 

-57,251 

-68, 209 

-65,750 

-1*8,850 

-27,181* 

-  9,71*7 

-  2,750 
+12,971* 
-15,035 

-  9,033 
-18,088 
-32,1*77 

-  2,1*35 
+10,21*1* 
+11,1*79 
+  2,371* 

-  9,11*1* 
-15,677 
-15,585 
-13,251* 
-13,920 
-11*,  776 
-16,538 
-38,992 

-  5,020 


0.03865 

0.21765 

0.699 

1.613 

2.963 

U.590 

6.163 

7.356 

7.953 

7.952 

7.51a 

6.851 

5.929 

U.830 

3.770 

2.981* 

2.5U7 

2.390 

2.31*6 

2.555 

2.313 

2.1678 

1.8768 

1.351*8 

1.3156 

1.1*806 

1.6656 

1.7038 

1.5568 

1.301*8 

1.0538 

0.81*03 

O.6163 

0.378 

0.112 

-  .515 

-0.596 
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D5  = 
0.003  v^ 


D*-Q 


D' 


je. 


TABLE  A  (cont'd.) 


DrDp 


(k£)         1+J  v^     (K£)(D5-Dp) 
(D^-D')  (1+Jv5) 


Vr 


V5  + 

z5  a     z5 


0 
0.00002 

0 

0    0.00002 

Uo 

0.00019 

0    0.00019 

380 

.00086 

0    O.OOO86 

1,720 

.00278 

0   0.00278 

5,560 

.00708 

0   0.00708 

lh,l60 

.01510 

0   0.01510 

30,200 

.02780 

0   0.02780 

55,600 

.0U501 

0   o.oU5oi 

90,020 

.06500 

0   0.06500  130,000 

.08U91 

0    0.08U91  169,820 

0.10207 

0.0207  0.0207  0.10 

200,000 

0.11559 

0.01559  0.01559  0.10 

200, 000 

0.127U7 

0.027U7  0.027U7  0.10 

200, 000 

0.13873 

0.03873  0.03873  0.10 

200,000 

.1U959 

O.OU959  O.OU959  0.10 

200, 000 

.15981 

0.05981  0.05981  0.10 

200, 000 

0.16909 

0.06909  0.06909  0.10 

200,000 

0.17733 

0.07733  0.07733  0.10 

200,000 

.18U22 

0.081422  0.08U22  0.10 

200,000 

0.19506 

0.09506  0.09506  0.10 

200, 000 

0.20028 

0.10028  0.10028  0.10 

200,000 

0.20U98 

0.10U98  0. 101+98  0.10 

200, 000 

0.20931 

0.10931  0.10931  0.10 

200,000 

0.21368 

0.11368  0.11368  0.10 

200, 000 

0.21785 

0.11785  0.11785  0.10 

200,000 

0.2219U 

0.1219U  0.1219U  0.10 

200, 000 

0.22616 

0.12616  0.12616  0.10 

200,000 

0.23063 

0.13063  0.13063  0.10 

200,000 

0.23531 

0.13531  0,13531  0.10 

200,000 

0.23998 

0.13998  0.13998  0.10 

200, 000 

0.2UU36 

0.1UU36  0.11*136  0.10 

200,000 

0.2U818 

0.1U818  0.1U818  0.10 

200,000 

0.25276 

0.15276  0.15276  0.10 

200, 000 

0.25567 

0.15567  0.15567  0.10 

200,000 

0.25695 

0.15695  0.15695  0.10 

200,000 

1.0008 

1.005UU 

1.022a 

1.0639 

1.1U3U 

1.267U 

1.U23U 

1.5737 

1.6662 

1.6636 

1.5721 

1.1*5  08 

1.3959 

1.3755 

1.3620 

1.3U06 

1.309U 

1.27U7 

1.2296 

1.2111* 

1.1739 

1.1567 

1.1511 

1. 11*56 

1.1391 

1.1361* 

1.11*05 

1.11*90 

1.1560 

1.1558 

1.11*60 

1.1273 

1.1527 

1.0969 

1.01*25 


0 
1*0.03 
382 

1,693 
5,920 
16,190 
38,250 
79,100 
11*1,700 
216,200 
282,300 
31U, U20 
290,160 
279,180 
275,100 
272, UOO 
268,120 
261,880 
251*,  91*0 
2U5,920 
21*2,280 
23U, 780 
231,31*0 
230,220 
229,120 
227,820 
227,280 
228,100 
229,800 
231,200 
231,160 
229,200 
225,1*60 
230,51*0 
219,380 
208,500 


U80 
2,960 
10,378 
25,907 
U9,360 
76,970 
96,790 
93,U20 
57,1*20 

-  1,600 
-56,860 
-75,U20 
-31*,  080 
*12,660 

-  8,380 
-13,320 
-19,360 
->21,560 
-31*,  220 
-11,320 
-23,280 
-10,660 

-  3,500 

-  3,1*20 

-  1*,0U0 

-  1,660 
+  2,520 
+  5,260 
+  U, 360 

-120 

-  6,120 
-11,600 

-i5,7ao 

-3a, 620 
-33,7ao 
-3a, 180 


0.007V 

0.055/ 

.222: 

.639! 
1.1*31* 
2.67a 

a.  23a 

5.737 
6.662 
6.636 
5.721 

a. 508 

3.959 
3.755 
3.620 

3.ao6 
3.09a 

2.7a7 

2.296 
2.11a 

1.739 
1.567 
1.511 
1.U56 
1.391 
1.36a 
1.1*05 
1.1*90 
1.560 
1.558 
1.1*60 

1.273 
1.527 
0.969 

o.a25 
-OJ.25 
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It  may  also  be  noted  that  permanent  ground  displacement  was 
not  recorded  until  it  exceeded  the  value  assigned  ground  quake  which 
was  0.10  in  this  problem.  The  time  intervals  at  which  these  occurred 

were 

i 
for  D  ,  time  interval  10 


D,  ,  time  interval  13 
Bl   ,  time  interval  1$ 


The  calculations  were  terminated  when  all  the  velocities  became  nega- 
tive j  effects  after  this  are  considered  to  be  secondary  and  of  little 
importance  for  the  purposes  of  this  study. 

A  graphical  presentation  of  the  results  is  shown  in  figures 
(1U)  through  (16)  inclusive. 

The  reason  that  a  time  interval  of  0.00025  seconds  was  chosen 
may  not  be  clear.  Generally,  it  may  be  said  that  the  best  time  inter- 
val to  use  is  one  that  produces  a  completely  stable  calculation.  This 
means  that  when  the  results  are  plotted  there  should  be  no  sharp 
wiggles  or  peaks.  If  such  irregularities  are  present,  either  an  error 
has  been  made  in  the  calculations  or  else  the  time  interval  chosen  is 
too  large.  Once  the  unit  lengths  of  the  pile  have  been  chosen,  the  time 
interval  must  be  selected  to  correspond.  Choosing  too  small  a  time 
interval  would  result  in  an  excessive  amount  of  calculation  with  only 
a  very  limited  increase  in  accuracy. 

Each  spring  in  the  figure  (13)  has  a  critical  time  interval 
which  is  the  time  that  is  required  far  a  stress  wave  to  traverse  the 
spring  and  its  associated  weight.  Stress  waves  travel  in  both  directions 
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at  a  speed  of  yE/p     where    P   is  the  mass  per  unit  of  volume. 


The 


following  equations  may  be  derived  for  the  critical  time  interval  for 
spring  K^  which  is  called  Tm. 

For  wave  motion  downward: 


.    .fU         i  -U/Vi 


For  wave  motion  upward: 


The  smallest  value  of  Tm  obtained  by  using  these  formulas  is  the 
critical  one.  About  one  half  of  the  critical  value  should  be  selected 
as  the  time  interval  to  be  used  in  any  particular  problem  so  as  to 
prevent  instability  from  arising  due  to  other  factors  not  included  in 
the  above  formulas  such  as  quake,  damping,  and  coefficients  of  restitu- 
tion. When  dividing  the  pile  into  ten-foot  sections  for  purposes  of 
mathematical  representation,  the  recommendations  of  Smith  (ref.  28) 
have  been  found  to  be  generally  satisfactory: 

For  steel  piles    —  0.00025  seconds 

For  wood  piles     —  0.00025  seconds 

For  concrete  piles  —  0.00033  seconds 
During  some  of  the  test  cases  run,  instability  resulted  in  using  this 
time  interval  for  wood  piles,  and  recomputations  were  made  using  smaller 
time  intervals.  Instability  can  generally  be  detected  by  plotting 
the  results  as  noted  previously  or  by  comparing  the  velocity  of  the 
pile  cap  and  the  pile  tip  with  the  velocity  of  the  ram  at  impact.  If 
either  of  these  velocities  exceeds  twice  the  ram  velocity  at  impact,  in- 
stability is  probably  present.  The  velocity  check  should  not  be  relied 
upon  alone,  however,  to  establish  the  stability  of  the  calculations. 
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IV.  COMPUTER  ADAPTATION 

A.  General 

As  shown  by  the  illustrative  example,   pile  driving  may  be 
analyzed  by  the  wave  equation  using  a  numerical  solution  method  and  a 
slide  rule.     This  approach  is  a  valid  one  and  will  produce  satis- 
factory results.     The  main  problem,   of  course,    is  the  amount  of  time 
and  work  required  to  make  the  calculations.     The  example  problem  was 
selected  as  being  about  as  simple  as  possible,  with  convenient  figures 
to  deal  with  numerically,   and  still  illustrate  all  aspects  of  the  cal- 
culations.    Nevertheless,   it  required  about  twelve  hours  of  manual 
computation  to  get  the  final  set.     It  should  also  be  noted  that  an 
error  made  at  the  beginning  of  a  solution  will  carry  throughout  the 
calculations.     The  usual  problem  encountered  in  practice  would  probably 
take  several  man-days  to  complete. 

Fortunately,   there  is  a  quicker  and  better  way  of  making  these 
calculations.     Use  may  be  made  of  high-speed  digital  computers  which 
can  solve  up  to  800  such  practical  problems  in  about  5  minutes,   de- 
pending upon  the  computer  used.     As  a  consequence,   use  of  the  wave 
equation  for  pile  driving  analysis  has  become  a  practical  matter. 

Although  it  is  not  the  purpose  here  to  go  into  the  details  of 
programming,   it  may  be  well  to  mention  some  of  the  broader  aspects  for 
those  totally  unfamiliar  with  computers.     For  a  computer  to  be  effective 
in  solving  engineering-type  problems^    the  problem-solving  procedures 
must  be  presented  to  the  computer  in  a  language  which  it  can  "understand." 
The  computer's  basic  language  consists  of  elementary  instructions  such 
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as  add  or  subtract.  The  problem-solving  procedure  must  be  translated, 
then,  into  instructions  which  the  computer  can  obey.  This  translation 
is  called  programming  and  can  be  carried  out  entirely  by  a  person,  or 
the  computer  can  assist  in  the  process  by  use  of  a  compiler. 

A  compiler  is  a  large  set  of  computer  instructions  which  can 
accept  a  problem-solving  procedure,  written  in  a  form  resembling  the 
procedure,  and  produce  from  it  the  proper  elementary  machine  instruc- 
tions that  will  solve  the  problem.  The  algebraic  compiler  system  used 
for  the  programs  here  is  called  Fortran,  which  stands  for  formula 
translation.  A  procedure  to  be  followed  in  solving  a  problem  with 
Fortran  is  specified  by  a  series  of  statements,  of  which  there  are 
several  types.  One  type  specifies  arithmetic  operations:  another  calls 
for  reading  a  data  card  or  printing  results ;   a  third  specifies  the 
sequence  in  which  the  statements  are  to  be  executed;  and  a  fourth  type 
is  made  up  of  statements  which  provide  information  about  the  procedure 
without  themselves  causing  any  action.  All  of  these  statements  taken 
together  form  a  source  program.  The  programs  included  in  the  appendices 
of  this  writing  are  source  programs.  When  punched  on  cards  and  put 
through  the  machine,  it  is  translated  by  the  Fortran  system  into  simple 
sets  of  machine  instructions  which  is  called  an  object  program.  For 
those  who  wish  to  pursue  the  subject  further,  reference  (21)  provides 
in  the  space  of  83  pages  all  the  information  necessary  for  a  person  to 

learn  to  develop  programs  of  his  own. 
B.  Description  of  Programs  Developed 

1.  Basic  Programs 

A  basic  program  in  Fortran  language  for  the  IBM  7090  computer 

was  prepared,  and  it  follows  the  method  of  calculation  illustrated  in 
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the  example  problem.  This  program  requires  that  input  data  be  supplied 
to  indicate  the  number  and  size  of  the  weights  in  the  mathematical 
model,  the  time  increment  to  be  used,  the  spring  constants  of  the  pile, 
the  spring  coefficients  of  the  soil,  the  coefficients  of  restitution 
of  the  capblock  and  cushion  block  (if  used),  the  velocity  of  the  ram  at 
the  instant  of  impact,  the  value  of  ground  quake,  values  of  point  and 
side  damping,  whether  or  not  the  pile  cap  can  transmit  tension,  and 
whether  or  not  side  resistance  is  present.  Output  of  the  basic  program 
includes  the  input  data  for  purposes  of  identification  and  checking, 
and  then  for  each  time  increment  the  force  and  resistance  applied  to 
each  block  of  the  model  together  with  displacement  and  velocity  of  each 
block  at  that  particular  instant.  The  permanent  set  of  the  tip  of  the 
pile  is  also  recorded.  The  program  is  arranged  to  stop  the  computer 
after  300  time- intervals  or  before  if 

(a)  the  velocity  of  the  pile  cap  exceeds  twice  the  velocity  the 
ram  had  at  the  instant  of  impact; 

(b)  the  velocity  of  the  pile  tip  exceeds  twice  the  velocity  the 
ram  had  at  the  instant  of  impact; 

(c)  the  velocities  of  all  the  blocks  are  simultaneously  negative 
or  equal  to  zero.  Velocities  toward  the  top  of  the  pile  are  taken  to 
be  negative. 

The  basic  program  can  be  used  to  study  stresses,  velocities, 
accelerations,  and  permanent  set  of  the  pile  as  they  vary  with  time. 
The  program  may  be  run  with  one  or  more  sets  of  data  at  a  time.  This 
program  in  its  entirety,  together  with  instructions  on  how  to  use  it, 
is  included  as  appendix  A. 
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Having  developed  this  program,  the  question  was  how  best  to 
use  it  in  trying  to  determine  if  correlations  could  be  achieved  with 
available  field  data  on  driving  records,  load  tests,  and  borings.  Of 
the  variables  that  must  be  provided  as  input  data,  the  greatest  un- 
certainty was  with  the  values  to  assign  ground  quake,  point  and  side 
damping,  and  distribution  of  resistance  along  the  sides  of  the  pile. 

2.  Vary  RU  Program. 

Correlation  was  approached  in  two  different  ways.  The  first 
method  was  particularly  suited  to  the  very  few  cases  in  which  careful 
and  complete  field  tests  were  reported  and  analyzed  comprehensively. 
In  these  cases,  reasonably  accurate  estimates  of  the  variables  except 
ground  quake  and  damping  factors  could  be  made.  Another  program  was 
developed  for  this  approach  which  was  basically  the  original  program, 
but  automatically  for  the  one  set  of  input  data  varies  ground  quake, 
ultimate  ground  resistance,  and  damping  factors  over  a  wide  range  of 
values.  The  output  includes  for  each  variation  of  basic  information 
indicated  the  ultimate  ground  resistance,  the  blows  per  inch  at  final 
set,  maximum  tension  and  compression  at  the  head,  mid-length,  and  tip 
of  the  pile,  the  time  interval  the  calculations  stopped  and  why,  and 
the  original  input  data.   Aiis  program  in  its  entirety  is  included 
with  instructions  on  how  to  use  it  as  appendix  B. 

3.  Researcher  Program. 

The  second  method  of  approaching  correlation  was  more  of  a 
statistical  approach,  and  it  xsras  designed  to  try  and  glean  as  much 
information  as  possible  from  published  field  data  which  were  not  as 
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comprehensively  reported  as  data  used  in  the  first  method  of  approach. 
References  (15)  and  (16)  were  used  extensively  in  this  attempt.  Certain 
pile  sections,  lengths,  and  driving  hammers  were  chosen  because  more 
reasonably  good  data  were  available  for  them  than  for  other  types. 
Another  computer  program  was  developed  around  the  original  program  and 
designed  to  systematically  produce  in  quantity  data  from  which  curves 
of  ultimate  ground  resistance  vs.  number  of  blows  per  inch  could  be 
plotted  for  various  values  of  damping  constants  and  ground  quake  and 
side  resistance.  This  program,  together  with  instructions  for  its  use, 
is  included  as  appendix  C.  It  is  designed  to  be  run  with  one  or  more 
sets  of  data.  Its  output  includes  for  each  of  the  above  automatic 
variations  of  input,  the  ultimate  ground  resistance,  values  of  soil 
spring  coefficients,  the  permanent  set,  the  blows  per  inch  at  the  final 
set,  maximum  tension  and  compression  at  the  head,  mid-length,  and  tip 
of  the  pile,  the  time  interval  the  calculations  stopped  and  why,  and 
the  original  input  data.  Thus  about  eighty  complete  problem  solutions 
are  obtained  for  each  set  of  data  run.  Ordinarily  three  sets  of  data 
can  be  run  at  a  time  in  this  manner  and  take  no  more  than  five  minutes' 
IBM  7090  computer  time.  Manual  solution  would  require  at  least  one  man- 
year  for  just  one  set  of  data.  Running  time  on  a  709  computer  would 
be  20  to  25  minutes,  but  only  about  a  minute  if  the  709U  is  used. 

U.  Validation  of  Computer  Programs 

Before  using  any  computer  program,  great  care  must  be  taken  to 
prove  that  it  is  correct  and  that  all  of  its  features  behave  as  they 
are  intended  to.  Making  a  manual  solution  is  the  only  way  to  be  cer- 
tain that  the  program  is  correctly  written.  To  prove  the  original,  or 
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basic,  program,  the  data  of  the  illustrative  example  presented 
herein  was  used  as  input,  and  the  solution  compared  with  the  manual 
solution.  A  tabulation  of  the  results  of  the  two  solutions  for 
displacements  of  the  ram  in  each  time  interval  are  found  as  Table  B. 
It  can  be  readily  seen  that  the  results  compare  as  closely  as  one 
would  expect  considering  that  a  slide  rule  was  used  in  the  manual 
computation.  As  a  further  check  of  the  program,  the  illustrative 
problem  presented  by  E.  A.  L.  Smith  (ref.  28)  was  computed  using  the 
basic  program,  and  the  results  plotted.  They  compared  in  every 
detail  with  Smith's  solution. 

As  additional  programs  were  developed  around  this  basic 
program,  they  were  tested  with  both  of  these  illustrative  examples 
and  verified  to  insure  that  subsequent  work  would  be  based  on  sound 
programs . 

5.  Computer  Program  for  the  Hiley  Dynamic  Formula 
The  Hiley  type  formula  is  thought  by  many  to  be  one  of  the 
better  dynamic  formulas  in  general  use  today.  Although  not  accepted 
by  all,  it  is  becoming  more  popular  in  the  United  States  and  has  been 
included  in  some  building  codes.  In  order  to  facilitate  comparison 
of  results  obtained  by  the  wave  equation  and  those  obtained  by  using 
the  Hiley  formula,  a  computer  program  was  developed  for  it  also. 
This  program  requires  that  input  data  be  provided  to  give  the  hammer 
data;  pile  material;  the  pattern,  either  rectangular  or  triangular, 
of  side  resistance  distribution;  coefficient  of  restitution;  weight 
of  the  pile;  area  of  the  pile  at  the  top,  mid-length,  and  bottom; 
the  ground  quake;  the  modulus  of  elasticity  of  the  pile  material; 
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Table  B 

COMPARISON  OF  COMPUTER  SOLUTION  WITH  MANUAL  SOLUTION 
FOR  ILLUSTRATIVE  EXAMPLE 


Displacement  of  Block  1 


Time  Interval 

Computer  Results 

Manual  Solution 

1 

0.030 

0.030 

2 

0.060 

0.060 

3 

0.089 

0.089 

k 

0.117 

0.118 

5 

o.ite 

0.1U5 

6 

0.171 

0.172 

7 

0.197 

0.197 

8 

0,221 

0.221 

9 

0.2U5 

0.2U5 

10 

0.267 

0.267 

11 

0.288 

0.288 

12 

0.309 

0.309 

13 

0.328 

0.328 

1U 

0.3V7 

0.3U7 

15 

0.36U 

0.365 

16 

0.381 

0.381 

17 

0.397 

0.397 

18 

0.U12 

0.1*12 

19 

O.U26 

O.U26 

20 

o*hho 

O.hhO 

21 

O.U53 

O.U53 

22 

O.U65 

0.1*65 

23 

O.U76 

0.1*76 

2k 

O.U85 

0.1*85 

2$ 

0.h9h 

O.U95 

26 

0.502 

0.502 

27 

0.508 

0.508 

28 

0.513 

0.51U 

29 

0.517 

0.518 

30 

0.520 

0.520 

31 

0.522 

0.522 

32 

0.523 

0.523 

33 

0.523 

0.521* 

3U 

0.522 

0.523 

35 

0.522 

0.521 

36 

0.517 

0.518 

37 

0.513 

o.5iU 

38 

0.508 

0.509 

39 

0.502 

0.50U 

END  OF  CALCULATIONS 
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and  the  pile  length  in  feet.  ._ The . program  will  automatically  range 
over  ultimate  ground  resistances  from  20  to  380  tons  in  increments 
of  20  tons,  writing  out  the  set  and  blows  per  inch  corresponding  to 
each  ground  resistance  as  well  as  repeating  the  input  data  for  each 
case.  It  computes  for  each  value  of  the  ground  resistance  the  set 
associated  with  all  the  load  carried  by  the  point,  one  half  of  the 
load  carried  by  the  point,  and  none  of  the  load  carried  by  the  point. 
Programming  the  Hiley  formula  may  not  be  justified  under  usual  cir- 
cumstances, but  it  was  done  in  this  case  as  a  convenience.  The 
program  and  instructions  for  its  use  are  included  as  Appendix  D. 
It  should  be  noted  that  the  Hiley  formula  is  discriminatory  to  long, 
heavy  piles  and  is  subject  to  the  same  general  limitations  as  any 
dynamic  formula. 
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V.  SOIL  ENGINEERING  CONSIDERATIONS 

As  has  been  mentioned  in  the  development  of  this  method  of 
pile  driving  analysis,  values  must  be  assigned  to  the  dynamic  soil 
properties  of  ground  quake  and  damping.  The  distribution  of  forces 
acting  on  the  pile  tip  and  sides  must  also  be  evaluated  to  complete 
the  information  required  for  a  solution.  Furthermore,  in  order  to 
approach  correlation  between  dynamic  driving  resistance  and  static 
resistance,  the  manner  in  which  frictional  resistance  along  the 
sides  of  the  pile  acts  in  both  the  static  and  the  dynamic  cases  must 
be  explored.  The  subject  of  load  testing  must  also  be  investigated 
insofar  as  the  way  in  which  it  is  done  and  the  interpretation  of 
results  affect  correlations. 

A.  Side  Friction  Under  Static  Loading 

1.  Granular  Materials. 

Some  recent  work  has  been  done  on  load  transfer  in  the  case 
of  end  bearing  piles  in  which  two  instrumented  H  piles  were  tested 
(ref.  12).  llhe  piles  were  1UBPQ9   and  li^P117  and  were  driven  ap- 
proximately UO  feet  through  granular  soil  varying  from  sandy  silt 
to  sand  and  gravel.  The  piles  were  end  bearing  on  weathered  silty 
shale.  In  spite  of  the  fact  that  the  piles  were  short  and  rela- 
tively stiff,  they  transferred  as  much  as  one  third  of  their  load 
to  the  surrounding  soil  by  friction.  Figure  17  shows  the  soil  con- 
ditions and  driving  records. 

Upward  shear  along  the  pile  develops  as  the  pile  shortens 
elastic ally  under  the  applied  static  load  at  the  top  of  the  pile. 
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This  movement  of  the  pile  causes  distortion  of  the  surrounding  soil 
resulting  in  development  of  a  frictional  force  resisting  the  dis- 
placement of  the  pile.  Generally,  it  would  be  assumed  that  the 
greatest  side  resistance  would  be  developed  near  the  top  of  the  pile 
and  would  diminish  with  depth.  This  action  may  be  visualized  by 
thinking  of  the  pile  as  a  spring  embedded  in  the  ground  with  a  load 
applied  at  the  top  of  the  spring.  One  might  also  expect  a  nearly 
linear  relationship  between  load  applied  and  load  transferred  as  long 
as  no  relative  movement  occurs  between  the  soil  and  the  pile.  How- 
ever, a  loading  is  reached  at  which  the  soil  begins  to  slip  along  the 
pile.  This  zone  of  shear  failure  begins  at  the  top  of  the  pile  and 
progresses  downward.  If  the  pile  tip  never  moved,  the  entire  applied 
load  would  be  carried  by  friction  along  the  pile  sides,  but  even  in 
the  case  of  progressive  downward  failure,  the  failure  zone  would 
never  extend  down  to  the  tip.  In  the  actual  case  the  pile  tip  does 
move  downward  as  the  applied  load  is  increased  so  that  it  accepts  a 
portion  of  it. 

The  test  results  reported  in  reference  1  were  presented  from 
a  viewpoint  of  over-all  transfer  of  load  to  the  soil  by  friction. 
Taking  the  data  presented,  and  dividing  each  of  the  test  piles  into 
sections  about  $   feet  long,  Table  C  has  been  prepared  in  an  effort 
to  show  the  variation  in  friction  force  on  each  section  of  the  pile 
as  the  test  load  was  increased.  In  each  of  the  test  piles,  section 
1  is  the  top  section.  A  study  of  this  table  will  show  the  validity 
of  the    description  of  the  interaction  between  pile  and  soil,  i.e., 
the  frictional  forces  have  a  tendency  to  develop  more  at  the  top  of 
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Increment  of  Total 
Applied  Load 


Section  of  U4BP89 
Pile  and  Soil 
Description 

1— Sand  &  Gravel 
2— Sand  &  Gravel 
3— Fine  to  Med. Sand 
b~Fine  to  Med. Sand 
5— Sand  &  Gravel 
6— Sand  &  Gravel 
7— Sand  &  Gravel 
8— Sandy  Silt,  tip 

in  weathered 

shale 


FRICTION  FORCE  ON  SECTION 
(IN  TONS) 


0-75  Tons         75-150  Tons 


Change  Total     Change  Total 
(aT"7£)  (atl£6) 


8 
7 
6 
0 
5 
b 
3 


150-225  Tons     225-300  Tons 
Change  Total     Change  Total 
(at"22T)  (at  300) 


8 
7 
6 
0 
5 
b 
3 


2 

2 
1 
3 
2 
0 
-1 


10 
9 
7 
3 
7 

b 

2 


5 
b 

-2 

k 

-3 

1 
5 


15 
13 
5 
7 
b 
5 
7 


b 

2 

3 
1 
b 
8 
5 


19 

15 

8 

8 

8 

13 
12 


10 


TOTAL  FRICTION  FORCE 

3k 

U8 

62 

93 

TOTAL  FORCE  ON  TIP 

Ul 

102 

163 

207 

Section  of  1UBP117 

Pile  and  Soil 

Description 

1— Sand  &  Gravel 

U 

h 

-1 

3 

-3 

0 

0 

0 

2— Sand  &  Gravel 

h 

h 

1 

5 

-3 

2 

0 

2 

3— Fine  to  Med. Sand 

0 

0 

b 

h 

-b 

0 

2 

2 

b--Fine  to  Med. Sand 

2 

2 

3 

5 

-5 

0 

1 

1 

5— Sand  &  Gravel 

3 

3 

1 

b 

13 

17 

11.5 

28.5 

6— Sand  &  Gravel 

3 

3 

1 

b 

b 

8 

-1.5 

6.5 

7— Sand  &  Gravel 

6.5 

6.5 

3.5 

10 

-1 

9 

5 

lb 

8— Sandy  Silt,  tip 

in  weathered 

shale 

0.5 

0.5 

1.5 

2 

b 

6 

-b 

2 

TOTAL  FRICTION  FORCE 

23 

37 

k2 

5b 

TOTAL  FORCE  ON  TIP 

52 

113 

183 

2b6 
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the  pile.  It  also  shows  that  local  Mlure  extends  over  a  wide  range 
of  total  applied  load.  For  example,  as  the  applied  load  is  increased 
from  75  to  150  tons,  the  friction  forces  on  section  7  of  the  Ik  BP  89  is 
decreased  by  one  ton;  from  150  to  225  tons'  applied  load, local  failure 
occurs  at  sections  3  and  5.  After  this  loading  the  granular  structure 
of  the  soil  seems  to  have  been  readjusted  by  the  shearing  forces  and 
displacements  so  as  to  have  an  increased  load-carrying  capacity. 

In  the  case  of  the  ll;  BP  117  pile,  as  the  applied  load  is  in- 
creased from  75  to  150  tons,  local  failure  has  taken  place  on  section 
1  only.  As  the  applied  load  is  increased  to  225  tons,  failure  continues 
on  section  1  and  is  expanded  to  include  sections  2,  3,  h,   and  7.  Upon 
further  increase  of  the  applied  load  to  300  tons,  it  is  evident  that 
failure  has  generally  progressed  down  to  and  including  section  k   even 
though  there  is  a  very  slight  increase  in  friction  force  on  the  top 
half  of  the  pile.  Sections  6  and  8  are  experiencing  failure  while 
sections  5  and  7  are  attempting  to  assume  the  load.  Even  at  this 
depth  it  is  interesting  to  note  that  section  5,  the  uppermost  of  these 
two  sections,  has  assumed  an  appreciably  higher  part  of  the  load. 

Both  of  these  piles  were  designed  for  150- ton  static  loads. 
When  test-loaded  to  150  tons,  it  is  noted  that  the  general  shape  of 
side  frictional  resistance  was  triangular  with  the  base  of  the  triangle 
at  the  top  of  the  pile  and  the  apex  toward  the  lower  sections.  In  the 
case  of  the  lU  BP  117  at  a  test  load  of  150  tons,  however,  the  general 
shape  of  the  frictional  resistance  acting  along  the  length  of  the  pile 
was  generally  of  a  rectangular  shape  down  almost  to  the  pile  tip.  Thee 
difference  may  be  explained  by  considering  the  loads  required  for  a 
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unit  axial  shortening  (spring  constant)  of  each  of  the  pile  shapes. 

AF 
From  the  expression  K  -  -—  ,  and  the  fact  that  the  length  and  modulus 

of  elasticity  of  each  pile  is  the  same,  their  spring  constants'  ratio 
will  vary  as  their  area  ratio  or  3U.UU/26.19.  It  can  be  seen  that  the 
1U  BP  117  is  1.32  times  as  "stiff"  as  the  lU  BP  89.  Since  the  over-all 
dimensions  of  the  piles,  the  applied  load,  and  the  surrounding  soil 
are  essentially  the  same,  this  relative  "stiffness"  must  account  for 
the  difference  found  in  the  pattern  of  frictional  resistances  developed. 
As  a  consequence  of  this,  it  is  interesting  to  note  that  the  lli  BP  117 
was  nearer  total  failure  at  a  test  load  of  300  tons  than  was  the 
lighter  lU  BP  89  pile. 

In  summary,  the  following  factors  may  be  important  in  con- 
sidering the  static  behavior  of  piles  in  granular  material: 

a.  The  time  which  has  elapsed  since  driving.  In  some  cases  a 
quick  condition  may  have  been  caused  during  driving  which  will  dis- 
appear as  pore  pressures  dissipate.  Relaxation  may  also  occur. 

b.  The  general  manner  in  which  load  is  transferred  to  the 

soil  may  vary  considerably  and  is  generally  dependent  upon  the  "stiffness" 
of  the  pile  in  resisting  axial  deformations  as  related  to  the  strength 
characteristics  of  the  soil.  The  less  "stiff"  piles  tend  to  give  up 
their  loads  to  the  soil  at  the  upper  strata  while  the  very  "stiff" 
piles  tend  to  carry  their  load  by  end  bearing, 

c.  Failure  of  piles  in  granular  materials  is  often  not  too  well 
defined  in  the  case  of  friction  piles.  Small  local  failures  along  the 
pile  may  result  in  a  decrease  in  load  carrying  capacity,  but  then  the 
soil  particles  in  undergoing  shear  may  assume  different  packings  and 
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become  denser  with  the  result  that  their  strength  has  been  restored 
or  even  increased  over  its  original  value.  The  pile  load  may  be 
increased  while  these  readjustments  are  taking  place  until  a  load  is 
finally  applied  which  is  sufficient  to  push  the  pile  through  the  soil 
continuously.  With  each  of  these  readjustments  the  soil  takes  on 
different  properties  causing  the  pile  to  undergo  variations  in  its 
bearing  capacities. 

d.  Unusual  conditions  which  may  be  pertinent  to  a  particular 
problem  such  as  a  quick  condition  or  hydrostatic  uplift. 

2.  Cohesive  Materials. 

As  in  the  case  of  piles  in  granular  material,  some  research 

work  has  also  been  done  with  instrumented  piles  in  cohesive  materials. 

H.  B.  Seed  and  L.  C.  Reese  report  in  reference  (2)  that  several  6-inch- 

diameter  pipe  piles  from  20  to  22  feet  long  were  driven  approximately 

15  feet  into  a  stratum  of  soft,  saturated  clay.  Electric  strain  gages 

were  installed  on  one  of  these  piles.  Test  loadings  to  failure  were 

made  at  various  intervals  after  driving  as  follows: 

Interval  After  Driving  Failure  Load 

(pounds) 

3  hours  1070 

1  day  2800 

3  days  3200 

7  days  5U00 

1U  days  5800 

23  days  6100 

33  days  6200 

This  data  show  that  the  ultimate  bearing  capacity  increased  very 

rapidly  with  time  for  the  first  week,  but  that  thereafter  the  increase 

was  very  gradual. 
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Most  of  the  load  was  removed  from  the  pile  by  side  friction 
with  only  about  10-15$  of  the  applied  load  reaching  the  pile  tip. 
The  increase  in  supporting  capacity  was  reflected  by  an  increase  in 
the  load  transfer  to  soil  at  all  depths  as  shown  by  Table  D  which  was 
derived  from  published  data  (ref.  26).  The  pile  was  divided  into 
three-foot  sections  with  section  1  being  at  the  top  of  the  pile.  From 
this  table  it  can  be  seen  that  the  resistance  is  generally  increasing 
with  depth  for  the  first  9   to  12  feet  of  embedment,  and  then  it  tends 
to  become  constant  with  depth  for  the  final  3  to  6  feet  of  embedment. 
It  is  interesting  to  contrast  this  behavior  with  that  already  described 
for  granular  materials. 

In  this  series  of  tests  relationships  were  established  for 
undisturbed  soil,  and  fully  remolded  soil.  Tests  of  soil  samples  taken 
next  to  the  pipe  wall  showed  that  the  strength  loss  caused  by  dis- 
turbance of  driving  the  pile  was  about  60$ to  70$  of  that  which  the  soil 
would  have  lost  due  to  complete  remolding.  The  clay  next  to  the  pile 
wall  subsequently  gained  strength,  and  at  the  end  of  the  test  period 
it  had  a  strength  of  60$  higher  than  the  strength  of  the  undisturbed 
soil.  This  increase  in  soil  strength  appeared  to  parallel  the  decrease 
in  excess  hydrostatic  pressure  indicating  a  relationship  between  the 
two.  For  all  values  of  time  at  all  depths  the  failure  of  the  friction 
pile  did  not  occur  in  the  soil,  but  at  the  interface  between  pile  and 
soil. 

In  another  series  of  tests  made  by  Vey  (ref.  35),  with  H  piles, 
he  found  that  in  some  cases  the  friction  force  reversed  itself  when  the 
applied  load  was  removed  and  held  the  stress  in  the  pile.  He  also 
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Table  D 

FRICTION  FORCE  ON  SECTION 
(IN  POUNDS) 

Test  7—33  Days  After  Driving 


Increment  of  Total 

0-2000 

lbs. 

2000-UOOO  lbs. 

I4OOO-6OOO  lbs. 

Applied  Load 

Change 

Total 
(for 

Change 

Total 
(for 

Change 

Total 

(for 

Section  of  6"  Pipe 

2000) 

Uooo) 

6000) 

Pile  and  Soil 

Description 

1 — Soft  sat.   clay 

110 

110 

0 

110 

50 

160 

2 — Soft  sat.   clay 

290 

290 

270 

560 

150 

710 

3~Soft  sat.   clay 

U80 

UffO 

Uoo 

880 

U20 

1300 

k — Soft  sat.   clay 

U80 

U80 

520 

1000 

500 

1500 

5 — Soft  sat.  clay 

1*75 

U75 

525 

1000 

560 

1560 

TOTAL  FRICTION  FORCE 

1835 

3550 

5230 

TOTAL  FORCE  ON  TIP 

165 

U50 

770 

First  Loading  to 

0-1070  lbs. 

Failure   (3  hrs. 

Change 

Total 

after  Driving) 

(for 
1070) 

1— Same  soil 

0 

0 

2~ 

200 

200 

3— 

0 

0 

U~ 

330 

330 

5- 

280 

280 

TOTAL  FRICTION  FORCE 
TOTAL  FORCE  ON  TIP 


810 
260 
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found  that  adherence  would  increase  between  the  soil  and  the  pile  with 
time,  and  reconsolidation  of  the  disturbed  clay  would  increase  its 
shear  strength. 

The  findings  of  these  two  series  of  tests  regarding  the  effect 
of  pile  driving  on  clay  cannot,  of  course,  be  assumed  to  be  valid  under 
all  conditions.  Tschebotarioff  has  shown  in  reference  (3U)  that  the 
remolding  effect  caused  by  driving  pipe  piles  closed  end  in  100  feet 
of  varved  clay  can  be  great.  Remolded  clay  layers  can  lead  to  the 
development  of  negative  friction  with  the  result  that  end  bearing 
piles  may  become  overloaded.  It  is  generally  understood  that  displace- 
ment piles  should  not  be  driven  in  highly  sensitive  clays;  therefore 
this  part  of  the  problem  will  not  be  discussed  further. 

In  summary,  it  may  be  said  that  the  static  behavior  of  piles 
driven  in  cohesive  material  may  depend  upon: 

a.  The  time  which  has  elapsed  since  driving.  In  some  cases 
the  increase  in  bearing  capacity  for  even  short  periods  of  time  may 
be  great;  in  other  cases  a  loss  of  bearing  capacity  may  result. 

b.  The  friction  force  acting  on  the  pile  tends  to  increase 
with  depth.  It  is  believed  that  the  load  transfer  would  decrease  with 
depth  if  the  shear  strength  of  the  clay  were  sufficient  to  accept  the 
loads  at  the  upper  strata  and  if  there  were  a  motivating  force  avail- 
able to  cause  these  upper  strata  to  come  into  intimate  contact  with 
the  pile.  Since  neither  of  these  conditions  usually  exists,  the  load 
transfer  is  found  to  generally  increase  with  depth. 

c.  Failure  of  piles  in  clay  may  be  more  dramatic  than  those  in 
granular  materials.  In  the  case  of  a  pile  not  resting  on  an  extremely 
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hard  strata  such  as  rock,  once  the  shear  strength  or  adhesion  has  been 
exceeded,  the  pile  may  begin  to  punch  its  way  through  the  supporting 
strata  and  continue  downward  until  the  load  is  reduced  or  the  pile 
encounters  a  firmer  strata. 

B.  Dynamic  Factors  Involved  in  the  Pile-Soil  Interaction 

The  main  obstacle  in  attempting  to  relate  the  penetration 
record  obtained  during  driving  of  a  pile  to  the  static  load  carrying 
capacity  is  that  of  obtaining  a  reasonably  accurate  estimate  of  the 
dynamic  resistance  to  penetration  between  the  pile  and  soil.  The  prob- 
lem is  the  same  whether  using  the  wave  equation  or  a  dynamic  pile 
driving  formula. 

The  dynamic  resistance  to  penetration  is  related  to  one  or 
more  of  the  following  factors: 

a.  The  static  resistance. 

b.  The  increase  in  strength  observed  to  occur  in  the  soil 
due  to  rapid  loading. 

c.  Frictional  resistance  between  the  pile  and  soil. 

d.  Viscous  resistance  developed  as  the  pile  penetrates  the 
soil. 

The  basic  s  trength  of  the  soil,  whether  due  to  cohesion  or 
internal  friction,  constitutes  the  static  load  capacity  of  the  pile 
and  soil.  At  least,  this  represents  the  limiting  value  since  the  pile 
may  not  develop  the  full  shear  strength  of  the  soil  when  acting  as  a 
friction  pile.  The  value  of  the  static  resistance  can  be  thought  of 
as  the  base  value  of  the  dynamic  resistance  in  the  idealized  case 
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where  the  strength  properties  of  the  soil  are  not  changed  by  the 
process  of  pile  driving.  This  is  rarely  the  case,  and  usually  the 
driving  of  the  pile  acts  to  decrease  or  increase  the  natural  strength 
of  the  soil. 

Several  investigators  have  studied  the  effect  of  rapid  loading 
times  on  the  compressive  strength  of  the  soil  (ref.  U,  11,  2$,   37), 
and  some  discussion  of  the  results  of  such  investigations  appears 
below.  It  seems  logical  to  assume  that  in  the  early  time  intervals 
of  pile  penetration  most  of  the  dynamic  resistance  is  a  result  of 
this  type  of  phenomenon,  the  magnitude  of  which  would  be  directly 
related  to  the  static  strength  of  the  soil. 

Frictional  resistance  in  nearly  all  cases  develops  between  the 
pile  sides  and  soil  as  a  result  of  the  penetration  of  the  pile.  The 
nature  and  magnitude  of  this  frictional  resistance  depends  upon  the 
type  of  soil,  moisture  content,  type  of  pile  and  pile  material.  If 
the  adhesion  or  friction  acting  between  pile  and  soil  is  greater  than 
the  shear  strength  of  the  soil,  the  frictional  resistance  may  be  de- 
veloped within  the  soil  itself.  The  dynamic  frictional  resistance 
would  be  of  about  the  same  value  as  the  static  resistance  since  the 
dynamic  coefficient  of  friction  between  the  pile  and  soil  materials 
is  about  equal  to  the  static  coefficient  of  friction  (ref.  2U).  Fric- 
tion as  a  dynamic  factor  will  be  neglected  since  there  is  some  doubt 
as  to  the  extent  that  actual  relative  movement  takes  place  between  the 
pile  and  soil  at  the  pile-soil  interface. 

No  published  work  has  been  found  that  deals  with  viscous 
forces  as  developed  by  penetration  into  soils.  One  investigator 
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attributes  part  of  the  influence  of  high  rates  of  strain  on  the  increase 
in  compressive  strength  of  soil  to  viscous  action  (ref.  37).  Since 
viscosity  is  defined  as  a  property  of  internal  friction  (ref.  2),  it 
will  be  considered  as  included  as  one  of  the  factors  contributing  to 
the  increase  of  compressive  strength  experienced  by  a  soil  subjected 
to  high  rates  of  strain. 

During  the  driving  of  a  pile,  it  is  assumed  that  the  resistance 
developed  at  the  tip  and  along  the  sides  can  be  expressed  in  terms  of 
the  resistance  developed  under  static  loading.  This  is  to  say  that 
the  dynamic  resistance  to  driving  is  composed  of  the  static  resistance 
plus  an  increment  of  resistance  that  develops  under  dynamic  loading 
and  is  expressed  as  a  percentage  of  the  static  value.  The  problem  is 
to  determine  what  form  the  relationship  between  static  and  dynamic 
resistance  takes  and  how  to  quantitize  it. 

Another  problem,  which  will  be  taken  up  later,  is  how  much  of 
the  static  resistance  is  present  during  driving.  It  cannot  be  generally 
assumed  that  all  of  the  static  resistance  is  present  during  driving, 
and  Chellis  points  out  ways  in  which  this  may  vary  (ref.  6).  On  the 
other  hand,  dynamic  resistance  may  develop  in  soil  strata  that  offer 
only  very  small  or  negligible  static  resistance. 

In  the  computer  solution,  the  total  resistance  to  driving  at 
the  pile  point  is  represented  by  the  expression: 

R  =  (Dm  -  Dm)  K«  (1  +  Jv)  (15) 

where   (Dl,  -  D^)  is  the  "elastic"  strain  of  the  soil  with 

a  limiting  value  of  Q, 

K1  is  the  spring  constant  of  the  soil, 


' 


' 


' 


52 


v  is  the  pile  tip  velocity  at  any  instant,  and 

J  is  the  damping  constant. 
This  relationship  is  shown  by  Smith  (ref.  28). 

Two  questions  arise  with  respect  to  this  mathematical  representa- 
tion of  the  dynamic  resistance  to  driving: 

(a)  Does  this  equation  approximate  the  actual  form  of  the  pile 
resistance;  i.e.,  does  the  resistance-time  relationship  approximate 
that  of  the  actual  case? 

(b)  What  values  of  J  should  be  assigned  for  different  soil  and 
driving  conditions? 

Smith  proposed  the  use  of  a  value  of  0.1$  for  J;   however,  since  the 
purpose  of  his  published  Xirork  did  not  include  the  correlation  of  re- 
sults with  load  tests,  he  gave  no  basis  for  the  use  of  this  value.  Some 
preliminary  attempts  to  correlate  wave  equation  results  with  the  results 
of  pile  load  tests  and  driving  records  using  values  of  0.1!?  and  0.05 
for  J  and  J'  indicated  that  this  value  for  J  was  not  appropriate  for 
clay  soils. 

In  an  attempt  to  answer  the  questions  posed  above,  a  search  of 
the  literature  was  made  for  any  experimental  work  that  might  lead  to 
a  confirmation  of  the  nature  of  the  dynamic  resistance  as  well  as  to  a 
determination  of  more  appropriate  values  for  the  coefficients  J  and  J*. 

Several  research  programs  over  the  past  15  years  have  investi- 
gated the  influence  of  time  of  loading,  or  strain  rate,  on  the  compres- 
sive strength  of  soils.  All  such  projects  have  used  rapid  loading  in 
triaxial  and  unconfined  compression  tests,  or  dynamically  loaded  footing 
models  as  the  method  of  investigation.  The  application  of  results  from 
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such  tests  to  pile  driving  must  be  approached  with  caution  and  should 
be  considered  only  as  indicative  of  the  general  nature  of  the  dynamic 
resistance  to  pile  driving  -  and  of  the  general  order  of  magnitude  for 
values  of  the  coefficients  of  damping. 

Before  describing  the  tests  made  by  individual  investigators, 
several  terms  will  be  defined: 

(1)  Time  of  loading:  The  time  from  the  beginning  of  a  test 
to  that  at  which  the  maximum  force  is  achieved.  Figure  (21)  shows 
this  quantity  graphically. 

(2)  Strain  rate:  The  rate,  expressed  in  per  cent  per  second, 
of  strain  imposed  on  the  sample  at  the  time  the  compressive  strength 
is  reached. 

(3)  Transient  load:  A  dynamically  applied  load  in  an  uncon- 
fined  compression  or  triaxial  test  that  reaches  a  maximum  in  times  of 
less  than  30  seconds  and  as  short  as  0.005  second. 

(U)  Static  test:  For  comparative  purposes,  a  test  that  is  made 
with  a  time  of  loading  of  ten  minutes  or  more,  unless  otherwise  indicated. 

(5)  Strength  ratio:  The  ratio  of  the  compressive  strength  for 
a  given  time  of  loading  to  the  strength  for  a  static  test. 

(6)  Strain-rate  effect:  Same  as  strength  ratio. 

2.  Work  at  Harvard 

In  a  test  program  sponsored  by  the  Panama  Canal  in  1QU8, 
Casagrande  and  Shannon  made  a  series  of  static  and  transient  unconfined 
and  triaxial  tests  on  soft  to  medium  clays  and  a  soft  rock  (ref.  k) . 
Triaxial  tests  of  the  vacuum  type  were  also  made  on  a  dry  sand.  Times 
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of  loading  for  the  transient  tests  ranged  from  0.01  seconds  to  30 
seconds. 

In  all  the  tests  on  clays,  the  transient  compressive  strength 
was  greater  than  the  compressive  strength  as  determined  by  static  tests. 
The  transient  strengths,  taken  at  a  time  of  loading  of  0.02  second, 
were  from  l.U  to  2.6  times  the  static  compressive  strength.  The  tests 
also  showed  that  the  increase  of  the  fast  transient  compressive 
strength  over  the  static  strength  is  greatest  for  specimens  at  the 
highest  water  content  and  least  for  the  specimens  at  lowest  water 
content.   (Figure  (17-A)  is  a  summary  of  the  average  results  from  the 
clay  tests. 

Vacuum-type  triaxial  tests  were  performed  on  a  clean,  medium 
sand  with  a  void  ratio  of  0.62  and  at  a  minor  principal  stress  of 
0.3  kg/cm  .  The  test  results  indicated  that  the  compressive  strength 
of  a  dry  sand  in  fast  transient  tests  is  about  15$  greater  than  the 
static  compressive  strength. 
3.  Work  of  Whitman 

Whitman  measured  the  increase  in  compressive  strengths  of 
soils  using  a  triaxial  machine,  and  testing  at  high  rates  of  strain 
(ref.  37).  Table  E  shows  his  results  for  several  types  of  cohesive 
soils.  The  last  column  of  this  table  (S-R)  is  the  ratio  of  the 
failure  load  at  a  strain  rate  of  1000  per  cent  per  second  to  the 
failure  load  at  a  strain  rate  of  0.03  per  cent  per  second.  It  is 
observed  that  the  strain-rate  effect  varies  from  1.3  to  2.0  for  these 
soils.  In  his  tests,  Whitman  achieved  times  of  loading  as  short  as 
.005  second.  Confining  pressures  up  to  100  p.s.i.  were  used. 
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Table  E 

INCREASE  OF  COMPRESSIVE  STRENGTH 
IN  COHESIVE  SOILS  DUE  TO  HIGH  STRAIN  RATES 

(from  ref .  37) 

DESCRIPTION                 p.l.  P.L.    w    qu    S-R 

1.  Normally  consolidated,  sensitive 

ocean  sediment,  undisturbed         63  k9           92       0.3   2.0 

2.  Very  plastic  clay,  remolded        2?  38    U8    7   1.6 

3.  Plastic  clay,  remolded             17  H     16    10   1.7 

it.  Medium  soft,  slightly  sensitive 

clay,  undisturbed                 2I4.  26     27    10   2.0 

5.  Slightly  plastic  clay  loam, 

remolded                        23  22     21    13  l.k 

6.  Plastic  clay,  remolded             27  38    hh        1$      1.7 

7.  Moderately  sensitive  silty  clay, 

undisturbed                      21  22     35    22   1.6 

8.  Impervious  compacted  fill           17  11     12    25   1.8 

9.  Tough  compacted  fill               Ul  21     26    35   2.0 
10.  Stiff,  dry  clay,  undisturbed        23  30     20   250   1.3 


In  comparing  results  of  tests  on  clay  with  and  without  confining 
pressure,  Whitman  observed  two  different  time  effects.   The  first  effect 
can  be  likened  to  that  encountered  in  an  extremely  viscous  fluid,  and 
is  associated  with  a  continuous,  plastic  deformation  of  the  soil.  The 
second  effect  is  associated  with  the  formation  of  discontinuities  such 
as  shear  planes  and  splits,  aid  occurs  when  the  soil  sample  is  tested  to 
failure  without  (or  small  values  of)  confining  pressure.  From  these 
results,  he  concluded  that  when  cohesive  samples  were  subjected  to  a 
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lateral  pressure,  the  strain-rate  effect  was  the  result  of  viscous 
phenomena,  and  was  the  same  regardless  of  the  magnitude  of  strain. 
The  same  effect  occurred  in  very  plastic  soils  regardless  of  confining 
pressure. 

In  dry  or  moist  sands  an  increase  in  compressive  strength  of 
10-15  per  cent  was  observed.  Only  limited  tests  were  conducted  with 
saturated  sands,  and  figure  (18)  shows  results  for  saturated  Ottawa 
sand  and  a  well  graded  sand.  Whitman  observed  that  the  strain-rate 
effect  in  saturated  sand  must  result  from  differences  between  the  pore 
water  pressures  in  slow  and  rapid  tests.  When  sand  is  tested  at 
constant  volume,  pore  water  migrates  to  the  central  part  of  the  sample. 
As  speeds  of  deformation  increase,  more  energy  must  be  expended  to 
overcome  resistance  to  the  flow  of  water. 
h.     Work  of  Seed  and  Lundgren 

Seed  and  Lundgren  investigated  the  effect  of  transient  loadings 
on  the  strength  and  deformation  characteristics  of  saturated  sands 
(ref.  25).  Their  tests  were  also  made  with  the  triaxial  apparatus,  and 
three  different  categories  of  tests  were  conducted: 

a.  Static  tests  in  which  the  loading  times  were  about 
10  minutes; 

b.  Slow  transient  tests  with  loading  times  of  about 
U  seconds; 

c.  Rapid  transient  tests  with  loading  times  of  about 
.02  second. 

In  the  rapid  transient  tests  a  constant  rate  of  deformation  of  UO  inches 

per  second  was  used.  All  tests  were  conducted  wi th  a  confining  pressure 

of  2  kg/sq.cm.,  and  the  samples  were  tested  at  various  void  ratios. 
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Results  of  the  tests  on  dense,  fine  sand  indicated  that  the 
compressive  strength  of  a  sample  subjected  to  rapid  transient  loading 
was  about  UO  per  cent  greater  than  that  obtained  by  static  drained 
tests.  Due  to  the  high  rate  of  loading  the  rapid  tests  were  actually 
undrained  tests,  since  drainage  could  not  occur  within  the  very  short 
time  of  the  test.  About  half  of  the  increase  in  strength  is  attributed 
to  a  dilatancy  effect  and  half  to  the  high  rate  of  loading  alone. 

It  is  a  well  known  fact  that  when  a  dense  sand  is  sheared  it 
increases  in  volume.  If  shear  occurs  under  conditions  that  prevent 
the  raovement  of  pore  water,  such  as  in  the  undrained  triaxial  test  or 
in  very  rapid  loading,  the  increase  in  volume  will  result  in  a  decrease 
in  pore  pressure  and  an  increase  in  strength,  as  shown  by  the  equation 

s  ■  (z   -  u)tan0 
in  which  the  value  of  u  (pore  pressure)  is  decreased.  These  results 
are  summarized  in  figure  (19)  and  apply  to  fine  and  coarse  sands  for 
very  short  loading  times. 

If  the  results  of  the  tests  made  by  Seed  and  Lundgren  are 
applied  in  equation  (1$)   the  following  results  are  obtained: 

Dynamic  Resistance  (RD)  ■  1.1+0  x  Static  Resistance  (Rg) 
or,     RD  -  Rg{l  +  JV) 
where  J  is  the  damping  factor  and  V  the  velocity  of  deformation. 

1.U0  Rs  ■  Rs(l  +  JV) 

T  .  1.U0-1.00 

J ^ 

From  figure  (20),  the  load  curve  from  one  test,  the  velocity  of 
deformation  is  estimated  to  be  about  3$   inches  per  second,  and 
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substituting  this  value  in  the  last  equation  gives 

J  -  .133 
for  the  case  of  rapid  loading  in  saturated  sand.  This  value  agrees 
closely  with  the  value  of  0.1$  for  J  suggested  by  Smith  (ref.  28). 

5.   Work  of  Cunny  and  Sloan 

Some  preliminary  tests  were  conducted  by  the  Waterways  Experi- 
ment Station,  Corps  of  Engineers,  Vicksburg,  Mississippi  in  which  small 
footings  were  dynamically  loaded  in  a  test  tank  filled  with  clay  or 
dry  sand  (ref.  11).  In  these  tests  a  loading  machine  capable  of 
applying  a  load  in  a  time  of  from  3  to  150  milliseconds  was  programmed 
to  apply  a  preset  load  daring  a  predetermined  loading  time.  As  a 
control,  static  tests  were  also  performed  on  the  same  soil  under 
identical  conditions  except  for  loading  time. 

Dynamic  and  static  tests  were  conducted  on  a  sand  compacted  to 
a  relative  density  of  96  per  cent.  Figure  (21)  shows  the  relationship 
between  load  and  deformation  for  the  static  and  one  dynamic  test  on  a 
9-inch-square  footing.  The  load  time  curve  for  the  dynamic  test  is 
shown  in  figure  (22).  From  this  curve  it  is  seen  that  the  loading 
time  was  11  milliseconds  and  the  maximum  load  reached  was  3500  pounds. 
Since  the  maximum  value  of  the  load  was  limited  by  the  setting  of  the 
machine,  it  is  not  necessarily  the  full  strength  of  the  soil  under 
dynamic  loading. 

By  relating  the  values  of  static  and  dynamic  loading  at  equal 
deformations,  it  is  possible  to  arrive  at  an  indication  of  the  value 
of  the  damping  factor  under  these  conditions  of  loading.  At  a  deformation 
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of  0.10  inch  the  static  load  is  found  to  be  1600  pounds  and  the  dynamic 
load  26^0  pounds  from  the  load-deformation  curve  of  figure   (21).     The 
velocity  of  deformation  at  this  point  is  estimated  to  be  about    UO  inches 
per  second,    or  3.33  feet  per  second.     The  corresponding  value  of  J 
using  equation  (if?)  as  before  is  found  to  be  0.197. 

6.     Comparison  of  Dynamic  Tests  with  Computer  Solutions 

In  order  to  determine  whether  the  dynamic  tests  described  above 
are  qualitatively  comparable  with  the  equations  used  for  computation 
of  the  total  resistance,   an  analysis  has  been  made  of  several  sets  of 
computer  results  in  the   same  terms  as  the  data  obtained  from  the 
laboratory  results.     This  analysis  could  only  be  carried  to  the  point 
of  maximum  deformation,   since  the  program  was  established  to  stop  at 
this  point. 

Figures   (23)  and  (2k)  are  plots  of  the  computer  results   for 
pile  tip  resistance-time,   deformation- time,   and  load-deformation 
relationships  for  a  53-lb.  H  pile  I4O  feet  long  and  a  concrete  pile  36 
inches  in  diameter  and  176  feet  long.     The  H  pile  is  driven  with  a 
Vulcan  #1  hammer  and  the  concrete  pile  is  driven  with  a  Raymond  k/0 
hammer.     The  tip  resistance-time  curve  for  the  H  pile,   figure   (2k), 
compares  closely  in  shape  with  same  curve  for  the  dynamic  laboratory 
test  of  the  9-inch-square  footing  shown  in  figure   (22).     This  would 
tend  to  indicate  that  the  load-time  relationship  resulting  from 
equation  (1$)  is  a  fair  representation  of  the  actual  dynamic  load 
mechanism,   at  least  with  respect  to  a  dynamic  test  on  a  small  footing. 

The  pile  tip  velocity- time  curves  in  figures  (23)   and  (2k) 
show  a  velocity  of  about  2  feet  per  second  at  the  time  of  maximum 
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resistance  for  the  U 0-foot  H  pile  and  a  velocity  of  about  one  foot 
per  second  for  the  176-foot  concrete  pile.  If  it  is  accepted  that 
most  of  the  increase  of  dynamic  resistance  over  the  static  resistance 
is  due  to  the  increase  in  compressive  strength  of  the  soil  due  to 
rapid  loading,  and  if  this  increase  amounts  to  a  factor  of  1.0  for 
clays  and  0.U  for  sands,  then  the  damping  coefficient  must  have  a 
value  such  that  when  multiplied  by  the  velocity  of  deformation  at 
the  time  of  maximum  load  results  in  values  of  1.0  for  clays  and 
O.U  for  sand.  Figures  (23)  and  (21*)  indicate  that  tip  velocities  at 
the  time  of  maximum  resistance  are  between  1.0  ard  2.0  feet  per  second, 
and  based  on  these  values,  the  damping  factor  (J)  would  have  a  value 
of  0.5  to  1.0  for  clay  and  0.1  to  O.a  for  sand. 

If  load  vs.  deformation  is  plotted  from  the  information 
available  in  the  curves  of  load  and  deformation  vs.  time  for  the 
12  BP  53  pile  shown  in  figure  (23),  a  curve  shown  in  figure  (2$)  is 
obtained.  This  curve  is  very  similar  in  shape  to  the  curve  of  dynamic 
load  vs.  deformation  from  rapid  load  tests  performed  on  sarid  in  the 
laboratory  as  shown  in  figure  (21).  The  fact  that  these  curves  are 
similar  in  shape  tends  to  indicate  that  the  variation  of  load  with 
deformation  in  the  wave  equation  solution  is  of  the  same  form  as 
that  occurring  in  laboratory  load  tests.  From  this  it  is  concluded, 
tentatively,  that  the  computation  of  total  driving  resistance  at  the 
tip  by  the  equation 

R  -  (Dp  -  D^  K«  (1  +  Jv) 
is  a  valid  representation  of  the  actual  physical  process  of  rapid 
penetration  of  the  soil  by  the  pile. 
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7.  Ground  quake,  which  has  been  defined  as  the  limit  of  elastic 
deformation  of  the  soil,  is  normally  thought  of  as  a  static  concept, 
but  it  is  an  essential  element  in  wave  equation  computations.  The 
elastic  ground  compression  occurring  during  static  loading  is  not 
necessarily  the  same  that  occurs  during  pile  driving,  and  no  attempt 
has  been  made  to  relate  the  two.  Therefore,  quake  as  used  in  this 
investigation  is  a  dynamic  property. 

A  method  of  measuring  quake  during  pile  driving  is  described 
by  Chellis  (ref .  6).  Essentially,  this  method  measures  the  total 
elastic  compression  of  the  pile  and  soil.  Pile  compression  alone  is 
computed  from  the  formula  C_  ■  R  L/aE,  where  P^  is  the  ultimate 
bearing  load  on  the  pile,  L  is  the  length  of  the  pile,  A  is  the  cross- 
sectional  area,  and  E  ±s   the  modulus  of  elasticity.  This  value  is 
subtracted  from  the  measured  compression  of  pile  and  soil  to  obtain 
the  value  of  quake.  This  method  of  measuring  quake  assumes  that  the 
pile  is  in  compression  simultaneously  along  its  entire  length,  which 
is  somewhat  erroneous  in  view  of  the  longitudinal  stress  wave  concept. 

In  extensive  studies  made  of  pile  driving,  Hiley  made  measure- 
ments of  quake  and  reported  that  quake  varies  considerably  for  different 
conditions  of  ground  and  may  average  0.1  inch  in  firm  gravel,  up  to 
0.2  inch  in  firm  clay  and  may  have  values  as  high  as  0.5  inch  in 
peaty  material  (ref.  17).  With  respect  to  driving  conditions,  Hiley 
summarized  his  observations  as  follows: 

Driving  Conditions    Quake  Stress  in  Pile 

Easy  driving  0.0  -  0.10  500  psi 

Medium  driving  0.10  -  0.20  1000  psi 

Hard  driving  0.10  -  0.30  1$00  psi 

Very  hard  driving  0.Q5  -  0.20  2000  psi 


■ 
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It  was  also  stated  that  quake  may  be  as  much  as  double  these  values 
if  soft  ground  is  immediately  below  the  pile  tip.  These  values 
for  quake  are  used  as  a  limiting  range  for  investigations  of  load 
tests. 
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Chapter  VI 
CORRELATION  OF  WAVE  EQUATION  RESULTS  WITH  PILE  LOAD  TESTS 

A.  Methods  of  Approach  and  Effects  of  the  Various  Parameters 

Although  there  are  many  load  tests  published,  there  are  com- 
paratively few  which  have  not  left  out  some  vital  part  of  the  in- 
formation required  for  a  thorough  evaluation  of  the  total  pile  driving 
and  loading  process.  Partly  due  to  this  limitation,  two  general 
methods  of  approaching  the  problem  of  correlating  wave  equation 
computations  with  load  test  data  were  tried. 

The  first  method  was  particularly  applicable  to  those  tests 
for  which  very  complete  information  was  published.  In  these  cases 
it  was  thought  that  most  of  the  variables  could  be  fairly  accurately 
evaluated  so  that  the  investigation  could  be  limited  to  a  trial  and 
error  approach  in  which  the  principal  unknowns  could  be  varied  over  a 
wide  range  of  values  and  tried  in  various  combinations.  Generally 
this  means  that  ground  quake  and  damping  constants  and  side  friction 
were  varied  with  all  other  factors  kept  constant  for  the  case  in 
question. 

The  second  method  of  approach  was  adopted  to  try  and  derive 
as  much  information  as  possible  from  published  load  test  data  which 
were  not  reported  as  completely  as  those  investigated  by  the  first 
method  of  approach.  It  also  considered  desirable  to  run  a  series  of 
data  to  study  more  fully  the  effects  of  the  various  parameters  on  the 
resistance-set  relationship  for  particular  piles  and  hammers.  It 
was  thought  that  both  of  these  objectives  could  be  met  by  selecting 
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a  particular  pile  type  for  which  a  number  of  reasonably  good  test 
data  were  available.  Furthermore,  it  was  desired  that  the  pile  meet 
the  following  requirements: 

(a)  Be  constructed  of  a  material  of  which  definite  information 
was  available  as  to  its  properties. 

(b)  Be  constructed  of  a  material  which  of  itself  would  least 
affect  the  soil  around  it  and  vice  versa. 

(c)  Have  uniform  properties  throughout  its  length. 

(d)  The  pile  as  driven  would  be  required  to  sustain  the  static 
loads  to  be  imposed  upon  it. 

The  reason  for  these  requirements  was  to  eliminate  as  many  of 
the  unknowns  from  the  problem  as  possible.  The  steel  H  pile  and  steel 
pipe  pile  each  satisfied  these  requirements  better  than  other  types. 
Because  there  were  more  usable  load  tests  readily  available  for  steel 
H  piles,  particularly  the  12-inch  H  bearing  pile  weighing  53  pounds 
per  foot,  this  shape  was  selected.  It  had  the  added  advantage  that 
most  of  the  records  available  indicated  that  one  type  of  hammer  had 
been  used,  the  Vulcan  Number  One;  however,  tests  were  available  with 
other  size  hammers  so  that  a  further  check  of  the  method  could  be  made. 

The  main  effort  has  been  directed  toward  establishing  correla- 
tion for  each  load  test  case  within  a  limiting  range  of  values  of  each 
of  the  variables  involved.  Correlation  of  wave  equation  computations 
with  the  actual  pile  load  test  involves  the  evaluation  of  five  different 
variable  factors.  These  are: 
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(a)  Ground  quake; 

(b)  Point  damping  factor; 

(c)  Side  damping  factor; 

(d)  Distribution  of  resistance  between  point  and  sides  of 
the  pile; 

(e)  Per  cent  of  static  resistance  acting  during  driving. 
These  factors  have  been  discussed  in  above  sections  with  the  exception 
of  the  last  named  variable,  the  per  cent  of  static  resistance  acting 
during  driving,  which  will  be  discussed  below. 

In  correlating  theoretical  calculations  with  load  tests  and 
the  corresponding  driving  records,  all  variables  except  one  must  be 
held  constant  while  varying  one  factor  at  a  time.  From  practical 
considerations  it  was  not  possible  to  treat  all  five  factors  alike  as 
variables.  Since  little  is  known  of  the  nature  of  the  side  damping 
factor  and  since  the  effect  of  side  damping  in  most  cases  is  relatively- 
small  compared  to  point  damping,  the  side  damping  factor  was  eliminated 
as  a  separate  influence  by  assuming  it  equal  to  one  third  of  the  point 
damping  factor  for  most  cases »  Smith  indicated  his  belief  that  side 
damping  would  be  less  than  point  damping  for  a  given  soil,  and  he 
recommended  a  value  which  was  one  third  of  the  point  damping  coeffi- 
cient (ref.  28). 

Although  up  to  five  different  distributions  of  resistance  were 
investigated  for  each  correlation  case,  this  factor  was  assumed  known 
when  it  came  to  the  actual  correlation  using  the  first  method.  For  the 
other  method,  distribution  was  treated  within  a  range  of  values.  The 
various  distributions  of  resistance  used  for  investigation  are: 
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(a)  End  bearing 

(b)  75  per  cent  end  bearing  and  25  per  cent  side  friction 

(c)  50  per  cent  end  bearing  and  50  per  cent  side  friction 

(d)  25  per  cent  end  bearing  and  75  per  cent  side  friction 

(e)  All  side  friction 

Since  the  pile  is  generally  divided  into  sections  5  to  10  feet  long 
for  purposes  of  calculation,  the  side  friction  acting  on  the  lower 
section  of  the  pile  is  indistinguishable  from  the  point  bearing  resist- 
ance, except  for  damping,  and  adds  to  it  for  all  practical  purposes. 

Three  methods  of  distributing  side  friction  along  the  embedded 
length  of  the  pile  have  been  used.  The  method  used  most,  and  the 
simplest,  assigns  a  uniform  value  of  resistance  for  a  unit  length  of 
pile  regardless  of  depth.  This  method  is  termed  rectangular  distribu- 
tion. In  another  method  the  side  friction  is  distributed  in  a  manner 
so  that  the  unit  resistance  is  directly  proportional  to  the  depth,  and 
is  called  triangular  distribution.  Both  of  these  methods  are  included 
in  the  computer  program,  and  selection  of  the  applicable  method  is 
made  by  use  of  a  control  statement  in  the  program.  Except  for  very 
long  piles,  the  difference  in  effect  on  the  resistance-set  relationship 
of  the  two  methods  is  small,  so  the  first  method  was  generally  used. 
The  third  method  of  distributing  frictional  resistance  is  by  estimating 
the  actual  resistance  developed  in  each  of  the  soil  strata  that  the 
pile  penetrates. 

The  per  cent  of  static  resistance  acting  during  driving  is  a 
factor  that  is  influenced  by  time  effects.  In  a  sensitive  clay,  pile 
driving  remolds  the  soil  resulting  in  a  reduction  in  strength  during 
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and  immediately  after  driving  of  the  pile.  If  the  load  test  is  con- 
ducted after  the  soil  has  had  time  to  regain  its  strength,  the  value 
of  the  static  resistance  should  be  reduced  when  making  the  pile-driving 
analysis.  Another  type  of  time  effect  may  result  when  a  pile  is  driven 
into  a  sand  layer  through  overlying  layers  of  soft  or  medium  clay. 
When  the  pile  is  test  loaded,  resistance  initially  developed  in  the 
clay  may  be  transferred  by  plastic  flow  to  the  tip  and  lost  from  the 
standpoint  of  static  load-carrying  capacity.  During  driving,  however, 
resistance  in  the  upper  clay  layers  will  be  present  and  will  appear 
to  result  in  heavy  damping.  Another  effect  is  the  liquefaction  of 
saturated  silts  due  to  a  quick  condition  being  created  by  the  driving 
action  and  resulting  in  greatly  reduced  resistance  to  penetration  of 
the  pile.  The  consideration  of  the  effects  of  such  factors  in  the 
calculations  is  very  approximate  and  without  complete  soil  data  cannot 
be  accounted  for. 

In  correlating  wave  equation  computations  with  load  tests  in 
which  sufficient  data  was  available  on  which  to  base  a  reasonable 
estimate  of  the  distribution  of  load  between  point  bearing  and  side 
friction,  the  approach  was  to  determine  the  various  values  of  quake 
and  damping  that  produced  a  resistance-set  curve  that  included  the 
point  representing  the  ultimate  load  and  final  penetration  of  the  test 
pile.  Hypothetically,  there  would  be  an  infinite  number  of  such 
curves,  but  a  limiting  range  of  values  was  used  for  these  variables 
that  was  determined  on  the  basis  of  Smith's  work  and  the  experimental 
work  on  the  dynamic  properties  of  soil  discussed  above.  The  range  of 
values  for  quake  and  damping  investigated  were: 


68 


(a)  Quake 

(1)  sard  0.05  to  0.20 

(2)  clay  0.0$  to  0.30 

(b)  Damping 

(1)  sand  0.10  to  0.20 

(2)  clay  0.1*0  to  1.00 

These  ranges  of  values  were  not  rigidly  adhered  to,  and  in 
some  instances  outside  values  were  investigated.  The  load  test  in- 
formation, input  data  used  for  the  computer  program  for  each  load 
test,  and  results  of  the  correlation  attempts  are  included  in  tables 
P  to  BB.  The  accompanying  figures  are  the  ultimate  resistance  versus 
set  (in  blows  per  inch)  curves  drawn  from  the  computer  results  and 
used  as  the  basis  for  correlation. 

Similar  curves  used  as  the  basis  of  correlation  by  the  second 
method  are  included  in  appendix  E.  The  length  of  piles  used  for  the 
calculations  was  chosen  to  generally  fit  the  load  test  data  available, 
but  since  the  effect  of  length  was  not  clearly  established  the  choice 
of  U0-,  80-  aid  120-foot  piles  was  somewhat  arbitrary.  Calculations 
were  made  for  12  BP  53  piles  in  these  lengths  using  a  Vulcan  Number  One 
hammer  with  a  hardwood  capblock,  a  pile  cap  weighing  700  pounds,  for 
all  permutations  of  the  ranges  of  ground  quake,  point  and  side  damping 
factors  and  side  friction  distributions  of  a  rectangular  shape.  For 
each  combination  of  the  variables,  set  in  blows  per  inch  was  computed 
for  each  value  of  ground  resistance  ranging  from  UO  to  and  including 
280  tons  in  20-ton  increments.  The  results  were  plotted  so  that  on 
any  graph  there  are  five  curves,  one  for  each  percentage  of  side 


69 


resistance  assumed.  For  any  one  curve  of  ultimate  resistance  vs.  set, 
all  other  variables  are  constant. 

These  results  show  the  following  generalized  effects  of 
changing  the  variables  indicated  below: 

(a)  Ground  quake.  Increasing  ground  quake  increases  the  set 
in  blows  per  inch,  i.e.  harder  driving,  for  a  given  ultimate  ground 
resistance  for  all  values  of  side  resistance  including  the  point 
bearing  case. 

(b)  Increasing  both  the  point  and  side  damping  increases  the 
set  in  blows  per  inch  for  a  given  ultimate  ground  resistance  for  all 
values  of  side  friction,  including  zero. 

(c)  Side  distribution— increases  in  per  cent  of  side  friction 
cause  dramatic  decreases  in  set  in  blows  per  inch  for  a  given  ultimate 
ground  resistance. 

(d)  Increased  pile  length — here  the  results  are  mixed.  In  the 
zero  friction  or  end  bearing  case,  the  ultimate  resistance  tends  to 
become  the  same  for  I4O-,  80-,  and  120-foot  cases  after  about  2h   blows 
per  inch.  Before  that  point  is  reached,  the  80-  and  120-foot  piles 
have  essentially  the  same  curve  of  ultimate  resistance  vs.  set,  and 
the  UO-foot  pile  for  a  given  resistance  within  this  range  tends  to  have 
an  appreciably  decreased  set  in  blows  per  inch.  For  about  $0%   friction, 
the  results  may  be  practically  identical  for  all  three  lengths.  For 
higher  values  of  friction  the  longer  piles  may  have  decreased  set  for 
given  resistance. 

(e)  Increased  pile  weight  may  result  in  decreased  set  in  blows 
per  inch  for  a  given  ultimate  resistance;  however,  in  the  lower  ranges 
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of  resistance  the  difference  may  be  slight.  For  exanple,  for  a  side 
resistance  of  2$  per  cent  and  all  factors  the  same  for  two  U5-foot 
H  piles  driven  with  a  Vulcan  0  hammer,  the  difference  in  set  for  a 
1U  BP  89  and  a  lU  BP  117  is  less  than  2  blows  per  inch.  The  values 
are  8  and  6.6  respectively, 

(f )  Increased  hammer  weight,  but  with  the  same  stroke,  results 
in  decreased  set  in  blows  per  inch  for  a  given  resistance  and  all 
other  factors  constant. 

(g)  A  softer  capblock  results  generally  in  increased  blows 

per  inch  for  a  given  resistance,  but  the  effect  is  really  quite  variable. 
Comparisons  were  made  for  an  ultimate  ground  resistance  of  100  tons  for 
the  following  piles  driven  with  a  Vulcan  Number  One  hammer: 

1)  Concrete  Piles — 18  inches  square  with  lengths 
of  UO,  80,  and  120  feet.  It  was  assumed  that  head  packing 
would  be  used,  but  this  same  soft  material  was  assumed  for 
the  capblock  also  and  the  results  compared  with  capblocks 
of  hardwood  and  micarta.  There  was  no  appreciable  differ- 
ence in  sets  for  any  of  these  conditions  for  any  of  the  con- 
crete piles. 

2)  Wood  Piles—for  the  80-  and  120-foot  lengths,  the 
difference  was  slight  for  the  three  capblocks.  For  the  U0- 
foot  pile  the  difference  between  a  micarta  and  a  hardwood 
capblock  were  slight,  but  the  very  soft  capblock  made  about 
10  blows  per  inch  difference  in  the  results. 

3)  Steel  piles— the  capblock  had  essentially  no  effect 
on  the  80-  and  120-foot  length  piles,  but  the  very  soft  cap- 
block  made  about  2  blows  per  inch  difference  with  the  UO-f oot 
pile. 
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The  effects  produced  by  varying  some  of  these  parameters  are 
not  surprising,  but  in  some  instances  they  may  be,  especially  if  one's 
judgment  is  based  on  results  produced  by  popular  dynamic  formulas. 
For  example,  the  results  produced  by  the  Hiley  formula  may  indicate 
very  slight  differences  between  the  cases  when  side  friction  is 
present  and  when  it  is  totally  absent.  On  the  other  hand,  it  is 
interesting  to  note  that  when  side  friction  is  a  large  percentage  of 
the  ultimate  resistance  to  driving,  the  set  vs.  ultimate  resistance 
curves  may  tend  toward  the  curve  produced  by  the  Engineering  News 
formula.  The  foregoing  results  produced  by  the  wave  equation  method 
of  analysis  would  seem  to  bear  out  its  general  suitability  under  widely 
varying  conditions  of  driving  equipment,  pile  type,  and  ground  conditions. 

As  can  be  seen  by  examining  the  following  case  studies,  some 
of  the  data  upon  which  to  base  correlation  is  rather  complete,  and 
some  is  not.  The  cases  studied  include  various  size  piles,  assorted 
pile  driving  equipment,  in  soils  varying  from  gravel  to  clay,  piles 
constructed  of  concrete  and  steel,  different  ratios  of  embedded  length 
to  driven  length,  and  for  varying  combinations  of  soil  strata.  It  is 
apparent  that  for  many  of  the  cases  correlation  can  be  obtained  in 
several  ways  by  various  combinations  of  values  for  ground  quake,  point 
and  side  damping,  and  side  friction.  Because  these  factors  are  not  now 
known  with  a  high  degree  of  certainty,  the  problem  becomes  one  of 
searching  for  a  pattern  within  the  likely  ranges  of  the  variables  for 
the  various  types  of  soils  encountered. 

In  some  of  the  following  cases  it  was  necessary  to  estimate 
the  values  of  Q  and  J  which  produce  correlation,  since  the  range  of 
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values  used  in  the  computations  did  not  extend  far  enough  to  include 
the  load  test  point  within  computed  load-set  curves.  The  values  for 
Q  and  J  so  estimated  are  based  on  the  curves  nearest  to  the  load  test 
and  the  trend  of  variation  between  curves.  It  is  believed  that  the 
estimated  values  are  reasonably  accurate  considering  the  nature  and 
accuracy  of  the  basic  data. 

Before  examining  the  case  studies  in  detail,  the  manner  of 
identifying  the  set  versus  resistance  curves  is  noted: 

Curve  A  —  0%   friction,  i.e.,  an  end  bearing  pile  with  no 
resistance  to  driving  offered  by  the  soil  along 
the  pile  sides. 
Curve  B  —  25%   friction,  i.e.,  75%   of  resistance  to  driving 
offered  by  the  soil  at  the  pile  point  and  25% 
along  the  pile  sides. 
Curve  C  —  50%   friction,  i.e.,  50%   of  resistance  to  driving 
offered  by  the  soil  at  the  pile  point  and  50% 
along  the  pile  sides. 
Curve  D  —  "]5%   friction,  i.e.,  25%   of  resistance  to  driving 
offered  by  the  soil  at  the  pile  point  and  75% 
along  the  pile  sides. 
Curve  E  —  10C$  friction.  All  resistance  to  driving 
distributed  along  the  pile  sides. 
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TABLE  F 


Case  Number  2001 

Reference:  (18) 

Pile:  Prestressed  concrete,  36  in.  diameter,  176  ft.  long 

Hammer:  Raymond  3/0,  12500  lb.  ram,   2U  in.   stroke 

Embedded  length:     U0  ft. 

Soil:  Soft  to  medium  clay  with  point  in  stiff  clay. 

Final  penetration:     20  blows  per  inch 

Load  test:     Tested  to  failure  at  190  tons  including  the  weight  of 
pile,   tested  one  week  after  driving 

Input  data:       W(l)  12,500    Ram 

W(2)  8,800     (includes  follower  and  driving  cap) 

W(3)-W(20)  U,955     (based  on  9.775  ft.   segments) 

V  10.15     (ram  velocity  at  impact  for  2 -ft. 

drop) 
K(l)  10,600,000     (capblock  spring  constant) 

K(2)  2,180,000     (cushion  and  pile) 

K(3)  -  -  (Pile) 

K(19)  20,750,000 

K'(m)  Soil  spring  constant  depends  on  value  of 

ultimate  resistance  which  is  varied. 

Distribution  of  resistance:     50$  end  bearing  and  50$  friction,  based 

on  computation  of  point  capacity  by  the 
formula:  P  *»  9AC.  Friction  considered 
to  have  a  rectangular  distribution. 

Results:         See  figures   (27)  and  (28) 

Correlation  is  obtained  for  the  following  values 
of  Q  and  J: 


JL 

J 

Jl 

0.1 

0.9 

0.3 

0.2 

0.3 

0.1 
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TABLE  G 

Case  Number  2002 

Reference:  (l8) 

Pile:      Prestressed  concrete,  36  in.  diameter,  1?6  ft.  long. 

Embedded  length:  UO  ft. 

Hammer:    Raymond  U/0,  15,000  lb.  ram,  3U  in.  stroke. 

Soil:      Soft  to  medium  clay  with  point  in  stiff  clay. 

Final  penetration:  Refusal,  last  3A  in.  at  99   blows. 

Load  test:  Tested  to  failure  at  230  tons  including  the  weight  of  pile. 

Input  data:  W(l)     15,000 

W(2)     8,800 

V       12.08 

(Rest  of  input  data  same  as  for  case  2001) 

Distribution  of  resistance:  Estimated  between  7$%   and  50$  point  bearing. 

Results:    See  figures  (29)  and  (30). 

Correlation  is  obtained  for  the  following  values 
of  Q  and  J; 

Q  J         J' 

0.2         0.8       0.27 
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TABLE  H 

Case  Number  2003 

Reference:  <l8) 

Pile:      Prestressed  concrete,  36  in.  diameter,  176  ft.  long. 

Embedded  length:  UO  ft. 

Hammer:     Raymond  3/0,  12,500  lb.  ram,  39  in.  stroke 

Soil:      Soft  to  medium  clay  with  point  in  stiff  clay 

Final  penetration:  10  blows  per  inch 

Load  test:  Tested  to  failure  at  165  tons  including  weight  of  pile. 

Input  data:  W(l)     12,500 

V         12.9 

(Rest  of  input  data  same  as  for  case  2001) 

Distribution  of  resistance:  Estimate  50$  end  bearing  and  50$  friction 

Results:    See  figures  (31)  through  (3k). 

Correlation  obtained  for  the  following  values 
of  Q  and  J; 


Q 

J 

j« 

0.1 

1.0 

0.33 

0.2 

0.5 

0.17 
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TABLE  I 


Case  Number  200ii 

Reference:  (18) 

Pile:      Prestressed  concrete,  36  in.  diameter,  80  ft.  long. 

Embedded  length:  1*3.5  ft. 

Hammer:     Raymond  3/0,  12, $00  lb.  ram,  39  in.  stroke. 

Soil:      Medium  and  stiff  clay  with  point  in  very  stiff  clay. 

Final  penetration:  20  blows  per  inch. 

Load  test:  Tested  to  failure  at  365  tons. 

Input  data: 


W(l) 

W(2) 

W(3)  -  W(10) 

V 

K(l) 

K(2) 

K(3)  -  K(9) 

el 

e« 


12,500 

8,800 

5,070 

12.9 

10,600,000 

k, 020, 000 

20,300,000 

0.8 

0.5 


Distribution  of  resistance:  From  computations  based  on  load  test  and 
soil  data,  it  is  estimated  that  87$  of  the  load  was 
carried  by  the  pile  point  at  failure.  Correlation  should 
be  made  for  distributions  between  100  per  cent  and  75  per 
cent  of  the  load  carried  by  the  point. 

Results:    See  figures  (35)  through  (37) 

Correlation  is  obtained  for  the  following  values 
of  Q  and  J: 


Remarks : 


Q 


0.1 
.05 


o.U5 
0.6 


J' 

0.15 
0.2 


In  this  case,  a  hole  2k   ft.  in  depth  and  k2   in.  diameter 
and  below  that  for  a  depth  of  27  ft.,  a  hole  35  in. 
diameter  was  pre-excavated  for  the  pile.   The  effect  of 
the  35  in.  pre-excavation  is  uncertain,  but  appears  to 
have  considerably  reduced  the  dynamic  resistance  to  driving, 
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TABLE  J 

Case  Number  2005 

Reference:  (18),  Test  7-1 

Pile:      Prestressed  concrete,  36  in.  diameter,  68  ft.  long. 

Hammer:     Raymond  3/0,  12, $00  lb.  ram,  39  in.  stroke 

Embedded  length:  31  ft. 

Soil:      Stiff  clay. 

Final  penetration:  kO   blows  per  inch. 

Load  test:  Tested  to  failure  at  260  tons. 

Input  data:  Same  as  case  200U,  except: 

W(3)  -  W(9)  U, 930 

K(3)  -  K(8)       20,900,000 

Distribution  of  resistance: 

From  calculations  based  on  load  test  and  soil  data  it  is 
estimated  that  k3%   of  the  load  was  carried  by  the  point 
of  the  pile  failure.  For  correlation,  $Qffo   of  the  load  is 
considered  to  be  carried  by  the  point.  Friction  is  con- 
sidered to  be  distributed  over  the  lower  31  ft.  of  the  pile, 

Results:    See  figure  (38). 

Correlation  is  obtained  for  the  following  values 
of  Q  and  J: 

Q         J         J' 
0.10       1.0        0.33 

0.20       0.5        0.17 
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TABLE  K 

Case  Number  6006 

Reference:  (l) 

Pile:      Prestressed  concrete,  5U  in.  diameter,  96  ft.  long. 

Hammer:     Raymond  3/0,  12,500  lb.  ram,  39  in.  stroke,  wood  cap block. 

Embedded  length:  60  ft. 

Soil:      Soft,  medium  and  stiff  clay  and  silty  clay  with  some  sand 
with  pile  of  point  in  dense,  fine,  gray  sand.  See  figure 
(39)  for  soil  profile. 

Final  penetration:  U5  blows  last  3A  in«  (60  blows  per  inch). 

Load  Test:  Tested  to  Li.20  tons,  failure  had  not  occurred.  Bayliss  cal- 
culated a  minimum  ultimate  load  capacity  of  U98  tons  based 
on  load  test  data  and  soil  characteristics  using  static 
formulas  (ref.  1). 


Input  data:                W(l) 

12,500 

W(2) 

6,000 

W(3)  -  W(12) 

6,200 

V 

12.9 

K(l) 

2,715,000 

K(2) 

2,025,000 

K(3)  -  K(ll) 

27,200,000 

el 

0.50 

e2 

0.50 

T 

0.00033 

Distribution  of  resistance: 

Based  on  Bayliss'  computations  (ref.  1)  it  is  estimated 
that  50  per  cent  of  the  load  was  carried  as  end  bearing 
with  friction  distributed  over  the  lower  60  ft.  of  the 
pile. 

Results:    See  figures  (39)  through  (1*2). 

Correlation  is  obtained  for  the  following  values 
of*  Q  and  J: 

Q  J  J' 


0.1  0.15  0.05 

0.15  0.10  0.05 

0.05  .2  .05       extrapolated 
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TABLE  L 

Case  Number  6007 

Reference:  (19) 

Pile:      Pipe  pile,  18  in.  diameter,  3/8  in.  wall  thickness,  80fb. 
long. 

Hammer:     Vulcan  No.  1,  5,000  lb.  ram,  36  in.  stroke. 

Embedded  length: 

Soil:  Predominantly  plastic  clays  to  a  depth  of  80  ft.  underlain 
by  uniformly  graded  fine  sand  having  a  relative  density  of 
approximately  80  per  cent.  Pile  tip  stopped  5  ft.  above  sand. 

Final  penetration:  2  blows  per  inch. 

Load  test:  Tested  to  failure  at  81  tons. 

Input  data: 


W(l) 

5,000 

W(2) 

1,000 

W(3)  - 

-  W(12 ) 

562 

V 

12.  k 

K(l) 

2,715,000 

K(2)  - 

-  K(ll) 

6,520,000 

el 

0.50 

e2 

1.0 

T 

0.00025 

Distribution  of  resistance:  Load  carried  as  friction. 

Results:  See  figure  (Ii3). 

Correlation  is  obtained  for  the  following  values 
of  Q  and  J: 

Q  J  J1 

0.1        0.15       0.05 
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TABLE  M 

Case  Number  6008 

Reference:  (19) 

Pile:      Steel  pipe  pile,  18  in.  diameter,  3/8  in  wall  thickness, 
96  ft.  long. 

Hammer:     Vulcan  No.  1,  5,000  lb.  ram,  36  in.  stroke. 

Embedded  length:  75  ft. 

Soil:      Plastic  clays  to  a  depth  of  80  ft.  underlain  by  a  uniformly- 
graded  fine  sand  having  a  relative  density  of  approximately 
80  per  cent.  The  pile  tip  was  driven  about  5  ft.  into  the 
sand  layer. 

Final  penetration:  226  blows  for  the  last  6  inches  (37.7  blows/in.) 

Load  test:  Tested  to  failure  at  2\\k   tons. 

Input  data:       W(l)  5,000 

W(2)  1,000 

W(3)  -  W(12)  676 

V  12.  h 

K(l)  2,715,000 

K(2)  -  K(ll)  5,U20,000 

ex  0.50 

e2  1.0 

T  0.00025 

Distribution  of  resistance: 

Based  on  two  different  load  tests,  one  with  the  pile  tip 
in  sand  and  the  other  (case  6007)  with  the  tip  5  feet 
above  the  sard  layer,  it  was  determined  that  159  tons 
or  6S   per  cent  of  the  load  was  carried  by  the  point  (the 
lower  5  ft.).  It  is  assumed  the  frictional  load  is 
distributed  over  the  embedded  length. 

Results:    See  figure  (kh). 

Correlation  is  obtained  with  the  following  values 
of  Q  and  J: 

Q         J         J» 


0.05     0.10     0.03 

0.10       0.05       0.05  (estimated) 


. 


■ 


. 
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TABLE  N 

Case  Number  6011 

Reference:  (20) 

Pile:      Hi-inch  H  pile  weighing  73  lb. /ft.  and  81  ft.  long. 

Hammer:     Vulcan  OR,  9,300  lb,  ram,  39  in.  stroke. 

Embedded  length:  80  feet. 

Soil:      Alternating  strata  of  silts,  sandy  silts,  and  silty  sands 

with  some  interspersed  clay  strata  for  a  total  thickness  of 
about  50  to  60  feet.  Clean  sands  from  UO  to  60  feet  thick 
lie  beneath  the  silts.  The  pile  penetrated  32  feet  into 
the  sand  layer. 

Final  penetration:  28  blows  per  foot  (2.3  blows/inch) 

Load  test:  The  pile  buckled  under  a  test  load  of  3^U  tons;  at  this 

point  the  pile  had  reached  a  gross  settlement  of  1.2  inches. 
From  the  load- settlement  diagram  it  appears  that  the  pile 
capacity  had  been  reached,  based  on  the  usual  criteria. 

Input  data:       W(l)  9,300 

.  W(2)  700 

W(3)  -  W120  6^8 

V  12.9 

K(l)  2,715,000 

K(2)  -  K(ll)  7,270,000 

e1  0.50 

e2  1.0 

T  0.00025 

Distribution  of  resistance: 

This  pile  was  equipped  with  strain  rods  in  order  that  the 
distribution  of  load  on  the  pile  could  be  determined  during 
the  load  test.  The  high  value  of  set  (low  value  of  blows 
per  inch)  for  this  pile  tends  to  indicate  the  possibility 
of  liquefaction  in  the  silt  layers.  Assuming  this  to  be 
the  case,  correlation  is  made  on  the  basis  of  the  load 
carried  to  the  sand  layer  which  was  170  tons. 

Results:    See  Figures  (U5)  and  (U6). 

Correlation  is  obtained  for  the  following  values 
of  Q  and  J: 

_Q_         J         _A1_ 
.1        .15        .05 


■ 


82 


TABLE  0 


Case: 

31 

Source ; 

Page  31,  reference  (16). 

Pile: 

6$*   12BP53 

Embedment:  55.  U  feet. 

Hammer:     Vulcan  #1 

Soil:      10'  silty  sand 

1U '  loam 

26'  fine  sand,  silty  sand,  and  clay 

10'  fine  to  medium  sand 

10'  fine  to  medium  sand  and  gravel 

then  fine  sand 

Final  penetration:  1.63  blows  per  inch. 

Load  test  results:  U5  tons  17  days  after  driving. 

Results  of  calculations  by  the-  wave  equation  method: 

a.  If  Q  =  0.10 

J  =  0.10 
J«=  0.03 
correlation  is  obtained  with  a  side  friction  of  2%. 

b.  If  Q  =  0.10 

J  =  0.20 
J'=  0.067 
correlation  is  obtained  with  a  side  friction  of  17$. 

c.  If  Q  =  0.20 

J  =  0.10 
J'=  0.03 
correlation  is  obtained  with  a  side  friction  of  2%. 

d.  If  Q  =  0.20 

J  =  0.20 
J'=  0.067 
correlation  is  obtained  with  a  side  friction  of  32$. 

Comments:   Hiley  Formula  gives  ultimate  resistance  of  52  tons. 
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TABLE  DD 


Case: 

3k 

Source : 

Page  3U,  reference  (16). 

Pile : 

12BP53 

Embedment : 

kl   feet. 

Hammer: 

Vulcan  #1 

Soil: 

35'  silty  clay- 
then  sand 

Final  penetration:  U.57  blows  per  inch. 

Load  test  results:  100  tons,  but  ultimate  resistance  not  reached. 

Results  of  calculations  by  the  wave  equation  method: 

a.  If  Q  =  0.10 

J  -  0.10 

J'  =  0.03 
correlation  cannot  be  obtained  for  these  values  for  a 
failure  load  of  100  tons;  however^  these  are  the  results: 

122  tons  for  no  side  friction 
lU6  tons  for  25%  side  friction 
185  tons  for  50%  side  friction 

b.  If  Q  =  0.10 

J     =     0.20 
J«   »     0.067 
correlation  is  obtained  with  no  side  friction. 

c.  If     Q     =     0.20 

J     =     0.10 
J'   =     0.03 
correlation  is  obtained  with  no  side  friction. 

d.  If  Q  =  0.20 

J  =  0.20 
J»  =  0.067 
correlation  is  obtained  with  a  side  friction  of  25%. 

Comments:  It  is  not  known  just  what  the  ultimate  resistance  would  have 
been  if  the  load  test  had  been  carried  all  the  way  to  failure, 

The  driven  length  of  the  pile  was  not  given;  57'  was  assumed. 

Hiley  Formula  gives  ultimate  resistance  of  96  tons. 


. 
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TABLE  P 

Case :      35 

Source:     Page  35,  reference  (16). 

Pile :      12BP53 

Embedment:  5U.U  feet. 

Hammer:     Vulcan  #1 

Soil:      15'  sandy,  silty,  clay  and  loam 

10'  sand 

10 *  sandy  loam 

2U1  sand 

then  sandy  loam 

Final  penetration:  U.57  blows  per  inch. 

Load  test  results:  100  tons,  but  ultimate  resistance  not  reached. 

Results  of  calculations  by  the  wave  equation  method: 

a.  If  Q  ■  0.10 

J  =  0.10 

J»  =  0.03 
correlation  cannot  be  obtained  for  these  values  for  a  failure 
load  of  100  tons;  however,  these  are  the  results: 

113  tons  for  no  side  friction 
lUl  tons  for  2$%  side  friction 
182  tons  for  50$  side  friction 

b.  If    Q     ■     0.10 

J     -     0.20 

J«     =     0.067 
correlation  is  obtained  with  a  side  friction  of  8%  for 
a  100-ton  ultimate  resistance.     Other  results  are: 

111  I  tons  for  25$  side  friction 
151  tons  for  50$6  side  friction 

c.  If  Q  =  0.20 

J  =  0.10 

J'  =  0.03 
correlation  is  obtained  with  a  side  friction  of  2%   for 
a  100- ton  ultimate  resistance.  Other  results  are: 

122  tons  for  a  2$%   side  friction 
163  tons  for  a  50^  side  friction 
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Case  35  (continued): 

d.  If  Q  =  0.20 

J  -  0.20 

J«  =  0.067 
correlation  is  obtained  with  a  side  friction  of  26$  for  a 
100-ton  ultimate  resistance.  Other  results  are: 

103  tons  at  2$%   side  friction 
129  tons  at  $C$   side  friction 

Comments:  It  is  not  known  just  what  the  ultimate  resistance  would 
have  been  if  the  load  test  had  been  carried  all  the  way 
to  failure. 

The  driven  length  of  the  pile  was  not  given;  65'  was 
assumed. 

Hiley  Formula  gives  ultimate  resistance  of  92  tons. 
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TABLE  Q 


Case: 

36 

Source : 

Page  36,  reference  (16) 

Pile: 

12BP53 

Embedment : 

U°.8  feet. 

Hanmer : 

Vulcan  #1 

Soil: 

15 '  sandy,  silty,  clay 

then  sand 
Final  penetration:  U.15  blows  per  inch. 

Load  test  results:  100  tons,  but  ultimate  resistance  not  reached. 
Results  of  calculations  by  the  wave  equation  method: 

a.  If  Q  =  0.10 

J  «  0.10 

J«  =  0.03 
correlation  cannot  be  obtained  for  these  values  for  a 
failure  load  of  100  tons.  Other  results  are: 

108  tons  for  no  side  friction 
135  tons  for  25$  side  friction 
178  tons  for  50/o   side  friction 


b.  If     Q     =     0.10 

J     =     0.20 

J'   =     0.067 
correlation  is  obtained  with  a  side  friction  of  1%  for 
a  100-ton  ultimate  resistance.     Other  results  are: 

115  tons  for  2%  side  friction 
1U3  tons  for  £0$  side  friction 

c.  If  Q  -  0.20 

J  =  0.10 

J»  =  0.03 
correlation  is  obtained  with  a  side  friction  of  1%   for 
a  100- ton  ultimate  resistance.  Other  values  are: 

115  tons  for  2$%   side  friction 
1^2  tons  for  $0$   side  friction 
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Case  36  (continued): 

d.  If  Q  =  0.20 
J  «  0.20 
J*  =  0.067 
correlation  is  obtained  with  a  side  friction  of  32$  fer 
a  100-ton  ultimate  resistance.  Other  values  are: 

91  tons  at  2$%   side  friction 
121  tons  at  5>0$  side  friction 

Comments:   It  is  not  known  just  what  the  ultimate  resistance  would 

have  been  if  the  load  test  had  been  carried  all  the  way  to 
failure. 

The  driven  length  of  the  pile  was  not  given;  60'  was 
assumed. 

Hiley  Formula  gives  ultimate  resistance  of  92  tons. 
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TABLE  R 

Case: 

67 

Source : 

Page  67,  reference  (16) 

Pile: 

llli1  10HU2 

Embedment:  95  feet 

Hammer:     Steam,  double  acting,  95>0  pounds  with  stroke  of  16  inches, 

Soil:      Silt  and  fine  sand,  pile  hit  hard  strata 

Final  penetration:  32  blows  per  inch 

Load  test  results:  60  tons 

Results  of  Calculations  by  the  wave  equation: 

a.  If  Q  =  0.10 

J  =  0.15 

J1  =  0.05 
correlation  cannot  be  obtained;  however,  100  tons  was  the 
calculated  value  for  no  side  friction. 

b.  If    Q     =     0.10 

J     =     0.20 

J'   =     0.067 
correlation  cannot  be  obtained;  however,    88  tons  was  the 
calculated  value  for  no  side  friction. 

c.  If     Q     =     0.20 

J«    =     0.05 
J     =     0.15 
correlation  is  obtained  with  a  side  friction  of  1$%. 

d.  If     Q     =     0.20 

J     =     0.20 
J>    =     0.067 
correlation  is  obtained  with  a  side  friction  of  1*0$. 

Comments:       Time  between  driving  and  load  testing  is  unknown. 

Hiley  Formula  gives  an  ultimate  resistance  of  25  tons. 
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TABLE  S 


Case: 

71 

Source: 

Page  71,  reference  (16) 

Pile: 

20 •  10HU2 

Embedment:  19  feet 

Hammer:     Gravity,  3,150-pound  ram  and  5-foot  drop 

Soil:      Sand  and  Gravel 

Final  penetration:  2  blows  per  inch 

Load  test  results:  6k   tons 

Results  of  Calculations  by  the  Wave  Equation: 

a.  If     Q     =     0.10 

J     =     0.15 

J»   =     0.05 
correlation  is  obtained  at  a  side  friction  of  both  75$ 
and  100$. 

b.  If     Q     =     0.10 

J     =     0.20 
J»   =     0.067 
correlation  is  obtained  for  a  side  friction  of  100$. 

c.  If  Q  =  0.20 

J  =  0.15 
J«  =  0.05 
correlation  is  obtained  for  a  side  friction  of  100$. 

d.  If  Q  =  0.20 

J  =  0.20 

J'  =  0.067 
correlation  cannot  be  obtained;  however,  the  ultimate 
resistance  is  calculated  to  be  60  tons  for  100$  friction, 
or  UO  tons  for  the  end  bearing  case. 

Comments:   Time  between  driving  and  load  testing  is  unknown. 

Hiley  equation  gives  an  ultimate  resistance  of  80  tons. 
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TABLE  T 


Case :      75 

Source:     Page  75,  reference  (16) 

Pile:      30'  10Hl*2 

Embedment:  23  feet 

Hammer:     Drop — 3000-pound  ram,  5-foot  drop 

Soil:      Fine  Sand 

Final  Penetration:  5  blows  per  2  inches 

Load  Test  Results:  90  tons 

Results  of  Calculations  by  the  Wave  Equation: 

a.  If  Q .  -  0.10 

J  =  0.15 
J'  =  0.05 
correlation  is  obtained  for  100$  side  friction. 

b.  If  Q  :  =     0.10 

J     =     0.20 

J'   =     0.067 
correlation  cannot  be  obtained  for  a  failure  load  of 
90  tons j  however,   these  are  the  results: 

$$  tons  for  no  side  friction 
6U  tons  for  50$  side  friction 
70  tons  for  100$  side  friction 

c.  If  Q  =  0.20 

J  =  0.15 

J'  =  0.05 
correlation  cannot  be  obtained  for  a  failure  load  of 
90  tons;  however,  these  are  the  results: 

5U  tons  for  no  side  friction 
60  tons  for  50$  side  friction 
60  tons  for  100$  side  friction 

d.  If  Q  =  0.20 

J  =  0.20 

J»  =  0.067 
correlation  cannot  be  obtained  for  a  failure  load  of 
90  tons,*  however,  these  are  the  results: 


. 
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TABLE  T  (cont'd.) 

U8  tons  for  no  side  friction 
55  tons  for  5Q&  side  friction 
60  tons  for  100$  side  friction 

Comments:   Time  between  driving  and  load  testing  unknown, 

Hiley  formula  gives  Ru  of  80  tons. 
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TABLE  U 

Case: 

76 

Source : 

Page  76,  reference  (16) 

Pile: 

UO»  10HU2 

Embedment :     3h  feet 

Hammer:  Drop — 3000-pound  ram,  5-foot  drop 

Soil:  Fine  Sand 

Final  Penetration:     2  blows  per  inch 

Load  Test  Results:     7f?  tons  on  second  loading,   £0  tons  first  loading 

Results  of  Calculations  by  the  Wave  Equation: 

a.  If  Q  =  0.10 

J  =  0.15 
J«  =  0.05 
correlation  is  obtained  with  side  friction  of  2$%. 

b.  If     Q     =     0.10 

J     =     0.20 

J«   -     0.067 
correlation  is  obtained  for  a  side  friction  of  35$  (the 
difference  in  this  case  between  2$%  and  $0$  friction 
is  only  1  ton). 

c.  If  Q  =  0.20 

J  «  0.15 
J'  =  0.05 
correlation  is  obtained  with  a  side  friction  of  1$%. 

d.  If     Q     =     0.20 

J     -     0.20 
J'    =     0.067 
correlation  is  obtained  with  a  side  friction  of  90$. 

Comments:       Time  between  driving  and  load  testing  unknown. 

Hiley  formula  gives  an  ultimate  resistance  of  60  tons. 
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TABLE  V 

Case :      89 

Source:     Reference  (12) 

Pile:      h$x   H1BP89 

Embedment:  Uk* 

Hammer:     Vulcan  0 

Soil:      Same  as  case  #  (l) 

Final  Penetration:    30  blows  per  inch 

Results  of  Load  Test:  300  tons 

Results  of  Calculations  by  the  Wave  Equation  Method: 

Correlation  is  obtained  with  2$%   friction  and 

Q  -  0.20 
J  =  0.20 
J»  -  0.067 


Comment:    This  pile  was  tested  $   days  after  driving.  The  300-ton 
load  appears  to  be  close  to  the  ultimate  resistance  of 
the  ground.  Load  increments  were  75  tons,  15>0  tons, 
225  tons,  225  tons  (repeated),  and  300  tons. 

Hiley  formula  gives  ultimate  resistance  of  230  tons. 
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TABLE  W 

Case :      113 

Source:     Page  113,  reference  (16) 

Pile:      kOx   precast  c oncrete  weighing  2.7  tons 

Embedment:  31.5' 

Hammer:     Vulcan  #1 

Driving  Cap:  0.525  tons 

Soil:      18'  coarse  sand 

26'  fine,  yellow  sand 

Final  Penetration:  10  blows  per  2.05  inches 

Load  Test  Results:  100  tons 

Results  of  Calculations  by  the  Wave  Equation: 

a.  If     Q     =     0.10 

J     =     0.15 
J'   «     0.05 
correlation  is  obtained  for  50$  friction. 

b.  If    Q     -     0.10 

J     =     0.20 
J"   =     0.067 
correlation  is  obtained  with  a  side  friction  of  85$. 

c.  If  Q  =  0.20 

J  =  0.15 

j»  =  0.05 

correlation  is  obtained  with  a  side  friction  of  75$. 

d.  If  Q  =  0.20 

J  =  0.20 
J«  =  0.067 
correlation  is  obtained  with  a  side  friction  of  95$. 

Comments:   Time  between  driving  and  load  testing  is  unknown. 

Results  using  the  Hiley  equation  are  122  tons. 
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TABLE  X 

Case: 

117 

Source: 

Reference  (12) 

Pile : 

\6X  1I4BP117 

Embedment : 

kh   feet 

Hammer: 

Vulcan  0 

Soil:       9'  sand  &  gravel 

12  *  fine  to  medium  sand 
15 *  sand  and  gravel 

61  sandy  silt 

2X  weathered,   silty  shale 
then  hard,   silty  shale 

Final  Penetration:     22  blows  per  inch 

Results  of  Load  Test:     300  tons 

Results  of  Calculations  by  the  Wave  Equation  Method: 

Correlation  is  obtained  with  2$%  friction  and 

Q  =  0.20 
J  =  0.20 
J»   =     0.067 


Comments:   This  pile  was  test  loaded  2  days  after  driving.  The 
300-ton  load  appears  to  be  close  to  the  ultimate 
resistance  of  the  ground.  Load  increments  were  75  tons, 
150  tons,  225  tons,  225  tons  (repeated),  aid  300  tons. 

Hiley  formula  gives  ultimate  resistance  of  lUO  tons. 
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TABLE  Y 

Case: 

lhh 

Source : 

Page  2hh,   reference  (16) 

Pile: 

120 •  12BP53 

Embedment : 

81.5* 

Hammer: 

Vulcan  #1 

Soil: 

51  clay  loam 

11'  sandy,  silty,  clay 

12*  medium  sand 

51  coarse  sand 

39'  medium  clay 

10'  stiff  clay 

10'  very  stiff  clay 

10'  medium  clay 

11'  hard  clay  with  pebbles 

then  medium  clay 
Final  Penetration:  8.0  blows  per  inch 
Load  Test  Results:  117  tons  ultimate  resistance 
Results  of  Calculations  by  the  Wave  Equation  Method: 

a.  If  Q  =  0.10 

J  =  0.10 

J»  =  0.03 
correlation  cannot  be  obtained  for  a  failure  load  of  117 
tons.  At  full  point  bearing  (0$  side  friction)  the 
ultimate  resistance  is  calculated  to  be  132  tons. 

b.  If     Q     =     0.10 

J     =     0.20 
J'   =     0.067 
correlation  is  obtained  with  a  side  friction  of  $%. 

c.  If  Q  -  0.20 

J  =  0.10 

J'  ■  0.03 

correlation  is  obtained  with  a  side  friction  of  0$. 

d.  If    Q     =     0.20 

J     -     0.20 
J»   =     0.067 
correlation  is  obtained  with  a  side  friction  of  20$. 

Comments:     Hiley  Formula  gives  ultimate  resistance  of  90  tons. 


. 
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TABLE  Z 


Case: 

11*5 

Source : 

Page  11*5,  reference  (16) 

Pile: 

106 «  12BP53 

Embedment :  1Q5 f 

Hammer:    Vulcan  #1 

Soil:      Same  as  case  #  l)|)| 

Final  Penetration:  10*3  blows  per  inch 

Load  Test  Results:  182  tons  ultimate  resistance 

Results  of  Calculations  by  the  Wave  Equation  Method: 

a.  If  Q  -  0.10 

J     =     0.10 
J'   =     0.03 
correlation  is  obtain?  d  with  a  side  friction  of  19$. 

b.  If     Q     -     0.10 

J     =     0.20 
J'   -     0.067 
correlation  is  obtained  with  a  side  friction  of  35$. 

c.  If    Q     =     0.20 

J     =     0.10 
J'   -     0.03 
correlation  is  obtained  with  a  side  friction  of  36$. 

d.  If    Q     ■     0.20 

J     =     0.20 
J'   ■     0.067 
correlation  is  obtained  with  a  side  friction  of  \\7%* 

Comment:     Hiley  Formula  gives  ultimate  resistance  of  86  tons. 


,  ■ 


' 
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TABLE  AA 

Case: 

1U6 

Source : 

Page  lk69   reference  (16) 

Pile: 

100'  12BP53 

Embedment: 

97.1* 

Hammer:  Vulcan  #1 

Soil:  Same  as  Case  #lMh 

Final  Penetration:     U.75>  blows  per  inch 

Load  Test  Results:         91  tons 

Results  of  Calculations  by  the  Wave  Equation  Method: 

a.  If     Q     =     0.10 

J     -     0.10 

J'  =  0.03 
correlation  cannot  be  obtained  for  a  failure  load  of 
91  tons.  At  full  point  bearing  (0$  side  friction)  the 
ultimate  resistance  is  calculated  to  be  110  tons. 

b.  If  Q  -  0.10 

J  =  0.20 
J'  =  0.067 
correlation  is  obtained  with  a  side  friction  of  2%. 

c.  If  Q  «  0.20 

J  =  0.10 
J»  -  0.03 
correlation  is  obtained  at  2%   side  friction. 

d.  If     Q     -     0.20 

J     =     0.20 
J'   =     0.067 
correlation  is  obtained  with  a  side  friction  of  17$. 

Comment:     Hiley  Formula  gives  ultimate  resistance  of  72  tons. 


' 
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TABLE  BB 

Case: 

161 

Source : 

Page  161,  reference  (16). 

Pile: 

50«  12BPS3 

Embedment:     32  feet 
Hammer:  Vulcan  #1 


Soil: 


3 

25 
18 
15 
9 
20 
22 


fill 

coarse  sand 
medium  sand 
fine   sand 
medium  clay- 
very  stiff  clay- 
medium  clay- 


then  very  stiff  clay 
Final  Penetration:     it.83  blows  per  inch. 
Load  Test  Results:     167  tons 
Results  of  Calculations  by  the  Wave  Equation  Method: 

a.  If    Q     =     0.10 

J  =  0.10 
J'  =  0.03 
correlation  is  obtained  with  a  side  friction  of  "£>%. 

b.  If    Q     =     0.10 

J     =     0.20 
J'   -     0.067 
correlation  is  obtained  with  a  side  friction  of  $5%. 

c.  If  Q  -  0.20 

J  •  0.10 
J«  =  0.03 
correlation  is  obtained  with  a  side  friction  of  $3%. 

d.  If  Q  ■  0.20 

J  =  0.20 
J>  =  0.067 
correlation  is  obtained  with  a  side  friction  of  70$. 


Comment:    Hiley  Formula  gives  ultimate  resistance  of  96  tons. 
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C.  Discussion  of  Results, 

Each  of  the  cases  correlated  by  the  first  method  are  dis- 
cussed below: 

1.     Case  number  2001,   table   (F).     This  and  the  next  two  cases 
are  36-inch  diameter  prestressed  concrete  piles  driven  at  Lake  Maricaibo, 
Venezuela.     The  load  carrying  capacity  of  the  pile  tip  was  calculated 
using  the  equation,   P  ■  9AC,   and  a  value  of  cohesion  (C)  of  2560  p.s.f ., 
giving  a  capacity  of  163,000  pounds.     This  value  is  U2  per  cent  of  the 
total  load  carried  by  the  pile  at  failure,   and  the  assumption  that  the 
load  is  carried  one  half  in  friction  and  one  half  by  point  bearing 
would  be  reasonable.     If  the  load  carried  as  friction  is  considered  to 
be  evenly  distributed  on  the  embedded  length  of  the  pile,   each  ten-foot 
section  would  develop  a  resistance  equal  to  k7f500  pounds.     If  the 
distribution  of  the  side  friction  resistance  were  assumed  to  vary 
directly  as  the  depth,    the  distribution  would  be  as  shown  in  figure 
(26),  which  also  compared  both  distributions  with  that  estimated  from 
the  shear  strength  of  the  soil.     In  the  calculated  values  the  full 
value  of  the  shear  strength  of  the  soft  clay  layers  was  assumed  to  be 
effective  and  a  reduced  value  of  800  p.s.f.  was  used  for  the  stiff 
clay  layer  between  elevations  -1!?3  and  -lii3  ft.     This  value  is  con- 
sistent with  Nordlund's  estimate   (ref.    23)  and  the  findings  of  others 
that  have  investigated  the  friction  values  of  piles  in  stiff  clay 
(ref.  33,   38).     Comparing  the  load  distribution  in  figure   (26),   it 
appears  that  the  rectangular  distribution  is  as  good  a  representation 
of  the  estimated  actual  load  distribution  as  the  triangular  method 
would  be.     From  figure   (27)  it  is  seen  that  values  of  Q  of  0.1  and 
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J  of  0.9,  and  Q  of  0.2  and  J  of  0.3  produce  nearly  identical  curves 
that  give  a  resistance  value  slightly  higher  (200  tons  vs.  190  tons) 
for  the  set  of  20  blows  per  inch. 

Case  2002.  This  is  the  same  pile  as  case  2001  redriven  to 
refusal  with  a  heavier  hammer.  The  load  capacity  was  not  increased 
a  great  deal  by  the  redriving,  going  only  from  190  tons  to  230  tons. 
Correlation  was  obtained  only  with  a  quake  of  0.20  using  a  value  of 
J  of  0.80  as  shown  in  figure  (29).  This  is  quite  a  departure  from  the 
value  of  J  obtained  in  the  previous  case  and  wouli  tend  to  indicate 
that  for  harder  driving,  quake  increased.  Correlation  could  not  be 
obtained  for  a  quake  of  0.1  with  values  of  damping  as  high  as  1.2. 

Case  200U.  The  tip  of  the  pile  in  this  case  penetrated  into 
clay  having  a  shear  strength  of  10,000  p.s.f.  as  determined  from 
unconfined  compression  tests.  Results  of  correlation  tend  to  indicate 
a  smaller  value  of  quake  than  experienced  in  the  previous  cases. 
Figure  (38)  shows  the  results  from  wave  equation  computations. 

Case  2003.  In  this  case  correlation  was  obtained  with  values 
of  Q  and  J  of  the  same  approximate  values  as  in  case  2001  as  indicated 
in  figures  (3U)  and  {36),     Figures  (31)  to  (33)  are  included  as  a  matter 
of  interest,  and  show  that  correlation  could  be  obtained  with  different 
distributions  of  resistance. 

Case  2005.  This  pile  is  of  the  same  type  as  the  previous  one 
but  was  driven  into  a  somewhat  softer  clay  having  a  shear  strength  of 
about  3000  p.s.f.  Figure  (39)  is  the  basis  of  correlation,  and  indicates 
values  of  J  in  the  range  more-or-less  expected. 
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Case  6006.  The  pile  in  this  case  was  not  loaded  to  failure, 
but,  from  analysis  of  the  load  test  and  driving  records,  the  engineer 
supervising  this  test  estimated  the  value  of  side  friction  and  point 
bearing  resistance  developed  under  static  loading.  From  this  the 
ultimate  load  carrying  capacity  of  the  pile  was  computed  to  be  a  mini- 
mum of  U98  tons.  This  value  was  used  in  correlating  the  load  test 
with  wave  equation  computations.  The  correlation  for  a  value  of  Q  of 
0.20  was  obtained  by  extrapolation,  but  it  is  believed  sufficiently 
accurate  considering  the  general  order  of  accuracy  of  the  other  vari- 
ables and  data.  Figure  (1*0)  is  the  load  settlement  curve  obtained 
from  the  load  test,  and  figures  (Ul)  and  (U2)  show  the  basis  for 
correlation. 

Cases  6007  and  6008.  These  are  two  pipe  piles  driven  through 
80  feet  of  plastic  clay  underlain  by  a  thick  sand  stratum.  One  pile 
was  driven  about  5  feet  into  the  sand  and  the  other  was  stopped  $   feet 
above  the  sand  layer.  The  object  of  the  load  tests  was  to  determine 
the  load  carried  by  the  point  in  sand.  Correlation  for  the  pile  that 
extends  into  the  sand  (case  6008)  was  obtained  in  the  lower  end  of  the 
expected  range  of  values  for  the  point  damping  factor  and  quake,  with 
values  of  0,0$  and  0.10  respectively.  No  attempt  was  made  to  account 
for  the  fact  that  the  pile  itself  is  largely  embedded  in  clay.  For 
the  pile  stopped  above  the  sand  layer,  correlation  was  obtained  for 
values  of  Q  of  0.10,  J  of  0,1$  and  J»  of  0.0$.  Since  the  load  is 
carried  entirely  by  friction,  all  of  the  damping  is  due  to  J',  the 
side  damping  factor.  These  values  are  not  consistent  with  the  expected 
range  of  values  for  clay,  and  indicate  that  the  value  of  quake  is  smaller 
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than  expected,  which  would  be  consistent  with  the  easy  driving  condi- 
tions, that  the  relationship  between  J  and  J'  is  in  error,  or  that 
side  damping  does  not  develop  when  a  steel  pipe  pile  is  driven  through 
plastic  clay.  With  respect  to  case  6008,  pile  point  in  sand,  it  is 
possible  that  the  clay  was  somewhat  sensitive  and  regained  strength 
during  the  period  before  the  load  test  was  made.  If  the  sensitivity 
of  the  clay  were  such  that  one  half  of  its  strength  were  lost  due  to 
remolding  and  then  regained  after  a  period  of  time,  the  ultimate  re- 
sistance during  driving  would  be  reduced  by  about  UO  tons.  This  would 
produce  correlation  at  a  value  of  Q  of  0.10  and  J  of  0.l£. 

Case  6011.  The  high  load  capacity  developed  by  this  pile  and 
the  low  resistance  to  penetration  indicate  that  a  quick  condition  in 
the  overlying  silt  layers  occurred.  If  this  were  the  case,  the  resistance 
to  driving  would  be  developed  on  the  portion  of  the  pile  embedded  in 
sand  with  no  resistance  occurring  in  the  silt  layers.  Correlation  is 
obtained  on  this  basis  as  shown  in  figures  (U5)  and  (U6). 

The  results  of  correlation  attempts  for  these  cases  are  sum- 
marized in  figure  (67).  The  piles  driven  with  the  tip  in  sand  are 
closely  grouped  and  fall  within  a  range  of  values  for  J  of  0.0S>  to  0.2, 
and  for  values  of  Q  of  0.05  to  0.10.  This  result  tends  to  confirm 
Smith  »s  recommended  value  for  J  of  0,1$,   for  piles  driven  in  sand. 

Results  of  correlation  for  piles  driven  in  clay  are  spread 
over  a  range  of  values  for  J  of  0.1!?  to  about  l.U,  and  for  values  of 
Q  of  0.1  and  0.2.  Case  6007,  which  has  a  J  value  for  correlation  of 
0.1$,  is  a  100$  friction  pile,  and  the  point  damping  factor  (J)  has  no 
influence.  This  possibly  indicates  that  the  one-third  relationship 


Case 


Test 
Load 


2001 
2002 
2003 
2004 
2005 
6006 
6007 
6003 
6011 


190 
230 
I65 
365 
260 
498 
81 
244 
170 


Blows 

per 
in«_ 


20 
132 

10 

20 

4o 
60 

2 
38 

2 


$  load 
at  point 


50 

50  -  75 

50 

50  -  25 

100  -  75 

50 
0 
75  -  50 
100  -  75 


Value 8  of  J  for 
Q=  0.05    Q=  0.1 


0.2 


0.1 


0.9 

1.0 

0.45 

1.0 

0.15 

0.15 

0.05 

0.15 


Q=  0.2 


Soil 


0.3 

0.8 

0.5 

0.15 

0.5 


Cloy 
Clay 
Clay 
Clay 
Clay 
Sand 
Sand 
Sand 
Sand 


0.2 


0.4         0.6 
Point  Lamping  (J) 


0.8 


Summary  of  Correlation  Results  -  First  Method 
Figure  (67) 
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between  J  and  J'  is  too  gross  an'  assumption  in  cases  where  g"H  or  most 
of  the  load  is  carried  by  friction.  In  general,  the  results  in  clay- 
are  consistent  with  values  of  J  expected,  based  on  dynamic  tests  made 
on  soils  in  the  laboratory. 

D.  Summary  of  Correlation  Results  for  Second  Method  of  Approach 

For  the  cases  studied,  the  following  values  of  ground  quake, 
point  damping,  and  side  damping  seem  to  be  appropriate  for  the  soil 
types  indicated.  The  percentages  of  side  friction  given  are  for  steel 
H  piles  unless  otherwise  indicated.  Furthermore,  the  principal  case 
studies  used  in  arriving  at  the  values  are  shown: 

1.  Coarse  Sand:  (Case  Numbers  113  and  l6l) 

Q  =  0.10 
J  =  0.15 
Ji  =  0.05 
Side  Friction  is  3>$%. 

2.  Sand  and  Gravel  Mixed:  (Case  Number  71) 

Q  =  0.10 
J  =  0.15 
J»  =  0.05 
Side  Firction  is  between  75$  and  100$. 

3.  Fine  Sand:  (Case  numbers  75,  76,  and  113) 

Q  =  0.15 
J  =  0.15 
J«   =     0.05 

Side  Friction  is  100$  (use  5$  for  straight-sided 

concrete  pile). 

U.  Sand  layers  combined  with  layers  or  strata  of  clay  or 

loam  or  both,  but  sand  layers  predominate:  (Case  numbers 
31,  3ks   35,  36,  HUi,  1U5,  11*6,  and  l6l  apply) 

Q     =     0.20 

J     =     0.20 

J»  =    0.067 

Side  Friction  is  2$% 
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5.  Silt  and  Fine  Sand  underlain  by  Hard  Strata:  (Case 
number  67) 

Q  =  0.20 

J  =  0.20 

J'  -  0.067 

Side  Friction  is  1|C$. 

6.  Sand  and  Gravel  Underlain  by  Firm  Strata  such  as  Shale: 
(Case  numbers  89  and  117) 

Q  ■  0.10  to  0.20 
J  ■  0.15  to  0.20 
J'  =  0.05  to  0.067 
Side  friction  of  2%. 

Using  the  average  values  of  Q,  J,  and  J'  and  side  friction 
obtained  for  the  various  types  of  soil  conditions,  the  variations  of 
calculated  ultimate  resistances  from  load  test  results  for  the  specific 
cases  will  be  investigated. 

Hiley  equation  results  are  also  shown  with  this  comparison  in 
Table  CC. 

It  should  be  noted  that  the  values  of  ultimate  resistance  shown 
to  be  calculated  by  the  Hiley  formula  were  obtained  by  using  values  of 
C-^,  C2,  and  Co  recommended  by  Chellis  (ref.  5)  for  use  in  the  absence 
of  specific  information.  No  attempt  was  made  to  determine  what  changes 
in  these  coefficients  would  be  necessary  to  achieve  closer  agreement 
with  load  test  results.  Comparison  of  Hiley  formula  ultimate  re- 
sistance calculations  has  therefore  been  restricted  to  those  case 
studies  where  complete  data  were  not  provided.  Whether  someone  having 
experience  with  the  Hiley  formula  could  have  justified  use  of  different 
values  for  these  coefficients  and  thereby  achieved  closer  agreement 
with  load  tests  is  not  known. 
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TABLE  CC 

COMPARISON  OF  WAVE  EQUATION  SOLUTION 
WITH  HILEY  FORMULA  AND  LOAD  TESTS 


Ultimate  Resistance 

Per  Cent 
Wave 
Equation 
% 

Error 

Case 
Number 

Computed  by- 
Wave  Equation 
(tons) 

Computed  by 

Hiley  Formula 

( tons ) 

Load 

Test 

(tons) 

Hiley 
Formula 
1  % 

31 

Uo 

52 

h$ 

-11.1 

-15.5 

3U 

100 

96 

100 

0.0 

-  k.o 

35 

103 

92 

100 

3.0 

-  8.0 

36 

91 

92 

100 

-  9.0 

-  8.0 

67 

60 

25 

60 

0.0 

-58.U 

71 

6k 

80 

6U 

0.0 

25.0 

75 

90 

80 

90 

0.0 

-ll.l 

76 

50 

60 

5o 

0.0 

20.0 

89 

300 

230 

300 

0.0 

-23.3 

117 

300 

mo 

300 

0.0 

-53.3 

113 

100 

122 

100 

0.0 

22.0 

1UU 

122 

90 

117 

U.3 

-23.0 

Ui5 

133 

86 

182 

-26.9 

-52.7 

1U6 

97 

72 

91 

6.6 

-20.8 

161 

103 

96 

167 

-38.3 

-U2.5 

. 
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It  must  be  remembered  that  the  number  of  cases  investigated 
here  is  small,  and  the  values  of  point  and  side  damping,  ground  quake, 
and  side  friction  reported  as  giving  good  correlation  may  be  replaced 
by  more  accurate  values  as  more  driving  records  and  load  tests  are 
analyzed  and  as  more  basic  information  about  the  dynamic  properties 
of  soils  becomes  known.  Furthermore,  side  resistance  was  taken  as 
rectangular,  or  uniform  distribution,  along  the  side  of  the  piles 
with  about  the  top  10  feet  of  embedment  discounted  in  each  case. 
Different  side  distribution  patterns  would,  of  course,  alter  the  values 
found  above  somewhat. 

On  the  other  hand,  the  values  reported  are  within  the  ranges 
that  one  might  anticipate,  and  they  do  give  good  correlation  as  well 
as  providing  results  generally  more  accurate  than  predictions  from 
the  Hiley  formula. 
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Chapter  VII 
RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

No  attempt  wxll^  be  made  here  to  make  an  exhaustive  and  detailed 
list  of  further  research  possibilities,  but  several  will  be  discussed 
briefly.  Probably  as  many  have  already  occurred  to  the  reader,  and 
certainly  many  more  will  present  themselves  to  anyone  beginning  work 
in  this  general  area.  The  main  factors  reqi  iring  further  evaluation 
are  ground  quake,  point  and  side  damping,  and  distribution  of  side 
frictional  resistance.  Specifically,  the  recommendations  are  as  follows: 

A.  Construct  an  apparatus  in  which  model  piles  may  be  loaded 
both  statically  and  dynamically.  If  economically  feasible,  the  model 
pile  may  be  instrumented  and  recording  pressure  cells  embedded  in  the 
soil.  Also  it  would  be  desirable  to  be  able  to  enclose  the  entire 
soil-pile  system  so  that  a  variable  confining  pressure  could  be 
applied.  With  this  set-up  the  model  pile  could  be  encased  to  avoid 
side  friction  effects  in  the  initial  experiments.  The  model  pile 
could  then  be  struck  by  a  small  drop  hammer  and  the  set  determined. 
This  procedure  could  be  repeated  for  different  sized  hammers  and 
strokes,  and  the  results  compared  with  computer  results  for  the  test 
set-ups.  The  ultimate  resistance  could  be  determined  from  test  load- 
ing statically  an  identical  pile  in  the  same  controlled  soil.  In  this 
way  values  of  quake  and  damping  might  be  bracketed  more  closely  for 
various  soil-pile  systems. 

B.  Obtaining  a  picture  of  the  interaction  between  pile  and 
soil  may  prove  useful.  It  is  thought  that  this  may  be  done  for  both 
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the  static  and  dynamic  cases.  A  photosensitive  material  could  be 
used  for  a  model  pile  which  would  be  embedded  in  another  photosensitive 
material  representing  the  ground.  The  stress  waves  propagated  by- 
striking  the  model  pile  at  its  head  by  a  small  drop  hammer  could  be 
captured  in  successive  time  intervals  by  using  a  high-speed  camera. 
Even  though  such  a  study  may  give  only  qualitative  results,  it  may 
assist  in  understanding  how  much  of  the  ground  is  acting  with  a  pile 
and  the  order  of  magnitude  of  stresses  within  this  zone.  Such  informa- 
tion should  prove  helpful  in  evaluating  both  ground  quake  and  the 
distribution  of  side  frictional  resistance. 

C.  Study  the  problem  from  both  an  analytical  and  statistical 
viewpoint.  Seek  the  cooperation  of  major  construction  corporations, 
government  bodies,  and  others  in  obtaining  as  many  good  load  tests 

as  possible.  For  example,  it  is  understood  that  the  University  of 
Michigan  or  the  Bureau  of  Public  Roads  is  on  the  verge  of  publishing 
a  number  of  very  carefully  controlled  tests.  Run  a  wide  range  of 
curves  for  various  pile  and  hammer  combinations  as  was  done  herein  for 
the  53-pound  bearing  pile,  and  see  if  dynamic  soil  properties  cainot 
be  deduced  from  them  which  will  give  good  correlation.  If  sufficient 
driving  records  and  load  tests  become  available,  it  may  be  advisable 
to  write  a  computer  program  which  will  assist  in  matching  up  wave 
equation  results  with  field  data. 

D.  A  wide  range  of  wave  equation  solutions  may  be  run  off 
for  the  common  soil  sampling  rigs  with  a  view  toward  correlating  the 
information  with  the  soil  data  and  with  pile  driving  and  testing 
records.  In  this  connection  it  should  be  noted  that  the  pipe  section 
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of  the  sampler  should  be  broken  down  into  one-foot  increments  and 
the  time  intervals  of  the  calculations  proportionately  reduced  in 
order  to  obtain  stable  calculations.  It  should  also  be  kept  in  mind 
that  since  lateral  support  of  the  long  pipe  section  is  lacking  in 
this  case,  some  energy  will  be  lost  in  transverse  vibrations.  How 
seriously  this  will  affect  the  results  is  not  known. 

E.  The  computer  program  included  herein  may  be  modified 
so  that  it  can  also  handle  the  case  of  a  long  hammer,  which  should 
be  represented  as  several  weights  in  the  mathematical  model,  and  the 
case  where  the  ram  strikes  the  pile  other  than  at  its  top.  A.  E.  L. 
Smith  illustrates  how  the  mathematical  model  for  these  cases  may  be 
represented  (ref.  28). 

F.  More  detailed  study  should  be  made  of  pile  tip  velocities 
at  the  time  of  maximum  tip  force,  as  determined  from  the  computer  re- 
sults, for  a  wide  range  of  pile  types,  lengths,  and  types  of  hammers 
in  order  to  establish  a  more  valid  basis  for  relating  results  of  rapid 
loading  tests  made  in  the  laboratory  with  pile  driving  action. 

G.  With  respect  to  the  computer  program,  investigation  should 
be  made  into  the  factors  that  cause  instability  in  the  computations. 
In  a  few  computer  results  for  this  work,  instability  was  encountered 
unexpectedly.  These  instances  involved  light  piles,  driven  with 
relatively  heavy  hammers,  or  low  values  of  ultimate  resistance  relative 
to  the  normal  capacity  of  the  pile.  Investigations  should  be  made  to 
determine  additional  criteria  for  use  in  making  up  the  input  data  so 

as  to  avoid  instability. 

As  previously  indicated,  these  are  but  a  few  of  the  possibilities 
for  further  work;  however,  it  is  hoped  that  these  listed  will  be  suffi- 
cient to  stimulate  the  interest  of  potential  investigators. 
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CONCLUSION 

The  prediction  of  a  pile,s  ultimate  static  bearing  capacity 
from  its  dynamic  behavior  during  driving  has  been  an  elusive  goal 
of  civil  engineers  for  many  years.  Hundreds  of  ingenious  empirical 
dynamic  formulas  have  been  devised  which  may  give  good  results  under 
particular  conditions  of  driving  equipment,  pile,  and  soil  combina- 
tions; but  the  difficulty  has  been  in  selecting  the  proper  formula 
to  obtain  reliable  results.  Furthermore,  these  dynamic  formulas 
have  limited  ranges  of  applicability  since  each  ignores  important 
aspects  of  the  problem.  Because  of  these  limitations,  their  validity 
as  a  general  method  of  approaching  the  problem  at  hand  is  open  to 
question.  Civil  engineers  have  recognized  these  shortcomings  and 
have  tended  to  treat  the  dynamic  formulas  as  crude  guides  to  the 
ultimate  bearing  capacity  of  a  pile  in  granular  soils,  recognizing 
their  complete  inadequacy  for  piles  driven  in  cohesive  soils. 

The  application  of  the  wave  equation  to  the  pile  driving 
problem,  as  envisioned  here,  represents  an  attempt  to  obtain  a 
general  method  of  analysis  which  will  assist  in  predicting  the 
ultimate  bearing  capacity  of  a  pile  from  its  observed  behavior  under 
the  last  hammer  blow  during  driving  subject  to  the  limitations 
enumerated  in  paragraph  E  below.  The  wave  equation  method  of  analysis 
is  also  useful  for  studying  stresses  occurring  during  driving  and 
for  selecting  appropriate  pile  driving  equipment  to  meet  economically 
specific  field  conditions.   The  major  contribution  by  A.  E.  L.  Smith, 
a  pioneer  in  this  work,  was  in  devising  a  mathematical  representation 
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of  the  pile  driving  problem  suitable  for  digital  computer  solution. 
He  then  developed  such  a  computer  program,  and  his  published  work 
which  followed  emphasized  the  stress  determination  problem.  He  did, 
however,  suggest  that  the  method  could  be  used  in  relating  the  static 
bearing  capacity  of  a  pile  to  its  dynamic  behavior,  and  he  offered 
values  of  damping  and  ground  quake  which  he  believed  to  be  accurate 
enough  to  make  use  of  the  method  practical. 

To  examine  the  possibility  of  correlating  pile  driving 
records  with  load  tests  by  this  method  of  analysis  has  been  the 
principal  objective  of  this  thesis.  Before  making  detailed  case 
studies,  it  was  necessary  to  examine  the  validity  of  the  method 
itself,  to  review  the  soil  engineering  aspects  of  the  problem,  and 
to  develop  the  necessary  computer  programs.  These  findings,  together 
with  the  results  from  2i;  case  studies,  are  as  follows: 

A.  Validity  of  the  mathematical  model  used.  The  model  de- 
vised by  A.  E.  L.  Smith  has  been  accepted  as  correct  by  those  com- 
menting on  his  work,  but  it  had  not  been  proven  to  be  equivalent 

to  the  wave  equation.  The  formal  proof  of  this  equivalency  has  been 
provided  in  Chapter  III. 

B.  Resistance  to  driving.  It  is  shown  in  Chapter  V  that 

the  expression  used  for  the  resistance  of  a  pile  to  driving  corresponds 
approximately  to  published  experimental  findings  of  others  on  the  re- 
sistance developed  in  a  soil  due  to  high  rates  of  loading. 

C.  Computer  programs.  Although  A  E.  L.  Smith  developed 
computer  programs  for  the  wave  equation  method  of  analysis,  they  were 
not  available  as  he  had  not  published  them.  Computer  programs  are 
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published  herewith  which  were  developed  independently.  Their  validity 
is  shown  by  comparing  with  both  a  manual  solution  for  a  sample  problem 
which  utilizes  all  parts  of  the  program  and  with  results  published  by 
A.  E.  L.  Smith.  Complete  instructions  for  using  these  programs  are 
provided  in  Appendices  A,  B,  and  C  so  that  someone  unfamiliar  with 
computer  work  may  utilize  them. 

The  computer  program  developed  for  the  Hiley  formula  is  also 
subject  to  the  limitations  inherent  in  an  empirical  dynamic  pile 
driving  formula. 

D.  Ground  damping,  quake,  and  side  friction.  From  published 
experimental  work  by  others,  approximate  values  of  damping  to  be  used 
with  the  wave  equation  method  of  analysis  for  some  soils  have  been 
derived  in  Chapter  V.  Starting  with  a  range  of  quake  reported  by 
others,  a  range  of  values  of  damping  were  found  which  caused  the 
solution  of  the  wave  equation  method  to  correlate  with  load  tests. 
This  range  of  values  of  damping  was  compatible  with  the  values  derived 
from  the  published  experimental  work.  Using  the  values  of  quake  and 
damping  thus  established,  values  of  uniformly  distributed  side  fric- 
tional  resistance  were  determined  for  various  soil  and  pile  combinations 

E.  Usefulness  of  the  wave  equation  method  for  predicting 
static  bearing  capacity.  Attempts  to  correlate  2k   pile  driving  records 
with  load  tests  have  shown  surprising  results.  It  appears  that  the 
wave  equation  method  of  pile  driving  analysis  may  become  an  accurate 
method  which  may  be  used  together  with  other  factors  in  predicting  the 
ultimate  static  bearing  capacity  of  a  pile,  as  indicated  by  the  results 
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summarized  in  figure  (67)  and  table  (CC).  Among  these  other  factors 
which  must  be  evaluated  are  the  relaxation  or  set-up  of  the  soil  with 
time  after  driving,  the  group  effect  of  piles,  development  of  negative 
friction,  consolidation  of  a  clay  layer  beneath  a  pile  tip,  long-term 
changes  in  water  table,  corrosion  or  deterioration  of  the  pile  from 
other  causes,  and  any  other  special  problems  which  may  be  associated 
with  a  particular  construction  project. 

Furthermore,  the  problem  has  been  treated  as  one-dimensional j 
i.e.,  lateral  vibrations  have  not  been  included  in  the  analysis.  This 
should  not,  however,  be  a  serious  limitation  on  the  method  for  most 
pile  driving  situations  since  a  large  proportion  of  the  pile  will  be 
embedded  as  final  penetration  is  reached.  Hysteresis  losses  in  the 
pile  have  not  been  included  in  the  analysis,  but  their  neglect  is 
considered  justified  at  this  stage  of  development  of  the  wave  equation 
method.  Since  good  correlation  was  possible  using  values  derived 
from  the  published  dynamic  soil  tests,  it  appears  that  the  wave  equation 
method  is  well  suited  for  analysis  of  any  size  or  type  of  pile  being 
driven  by  any  hammer  in  not  only  granular  soils,  but  also  cohesive  soils. 

F.  The  dynamic  interaction  between  pile  and  soil  is  incom- 
pletely understood  and  offers  a  fruitful  area  for  further  research.  Even 
a  qualitative  picture  of  stress  transfer  should  prove  helpful.  Addi- 
tional experimental  work  to  isolate  the  parameters  affecting  ground 
quake,  point  damping,  side  damping,  and  distribution  of  point  and  side 
resistance  to  driving  is  needed  as  well  as  the  evaluation  of  these 
factors. 
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G.  Correlation  with  many  driving  records  and  load  tests 
should  be  made  before  taking  the  values  reported  herein  as  correct. 
Future  investigators  may  use  freely  the  computer  programs  and  plotted 
results  of  calculations  already  made,  and  presented  here,  to  facilitate 
such  work. 
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APPENDIX  A 

DETAILS  OF  THE-  "BASIC  COMPUTER  PROGRAM" 
WHICH  USES  THE  WAVE  EQUATION 


General t  This  program  is  written  in  Fortran  language  for  the 
IBM  7090  computer,  but  it  can  be  used  equally  well  with  the  IBM  709  or 
IBM  709U  computer.  It  has  as  its  objective  the  detailed  analysis  of 
stresses,  velocities,  accelerations,  displacements,  and  permanent  set 
of  the  pile  as  these  factors  vary  with  time.  For  details  of  the  theory 
and  the  mathematical  model,  refer  to  chapter  III. 

Termination:  Analysis  is  terminated  automatically  when  the 
velocities  of  all  the  blocks  of  the  mathematical  model  are  simul- 
taneously negative  or  equal  to  zero.  Negative  velocities  are  taken 
as  being  in  the  direction  from  the  bottom  of  the  pile  toward  its  top. 
Termination  is  also  scheduled  when  the  permanent  set  of  the  pile  be- 
comes constant.  Effects  beyond  either  of  these  times  are  considered 
to  be  of  secondary  importance  for  purposes  of  this  study.  Generally 
these  conditions  will  be  met  before  calculations  have  been  made  for 
300  time  intervals,  and  in  most  cases  much  sooner.  The  time  interval 
is  usually  taken  as  0.00025>  or  0.00033  seconds  because  of  the  high 
velocity  of  stress  wave  travel  in  materials  used  for  piles ;   therefore, 
this  program  investigates  what  happens  during  the  first  0.1  second 
or  less  after  the  hammer  hits  the  pile.  In  the  usual  case  the  pile 
will  achieve  its  permanent  set  well  within  this  time. 

The  program  will  also  stop  automatically  if  the  velocity  of 
either  the  pile  cap  or  pile  tip  exceeds  twice  the  velocity  the  ram 
had  at  the  moment  of  impact.  This  feature  is  designed  to  help  detect 
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instability  (see  Chapter  III,  C)  in  the  calculations,  but  should  not 
be  relied  upon  entirely.  The  surest  way  of  detecting  instability  is 
by  plotting  the  displacements  and  velocities  of  each  block  in  the 
model  against  time  and  then  examining  the  curves  so  obtained  for  very 
sharp  peaks  or  discontinuities.  Should  such  irregularities  be  present, 
the  time  interval  may  be  reduced  and  the  program  rerun  and  rechecked 
for  stability.  The  importance  of  this  ^heck  cannot  be  over-emphasized. 

Should  the  program  run  the  full  300  time  increments,  no 
indication  is  given  in  the  output  as  to  why  the  program  stopped.  In 
such  a  case  the  results  should  be  examined  especially  carefully,  for 
in  all  probability  either  instability  was  present  or  the  program 
stopped  before  the  final  permanent  set  was  reached.  In  such  a  case 
the  variable  DP  may  be  redimensioned  to  a  larger  value  (second  state- 
ment above  statement  number  207  in  the  beginning  of  the  program),  the 
statement  immediately  after  statement  number  $00  should  be  changed  to 
read  DO  $01     I  =  1,X  where  X  will  be  the  new  dimension  of  variable 
DP,  and  the  statement  immediately  after  statement  number  302  should 
be  similarly  changed  to  read  DO  101  N  =  3,X.  In  all  other  cases  of 
program  termination  the  reason  for  stopping  will  be  printed  out  auto- 
matically in  plain  language;  for  example,  "VELOCITY  OF  PILE  CAP  EXCEEDED 
TWICE  THE  RAM  VELOCITY." 

Preparation  of  Input  Data:  All  data  is  punched  on  standard 
data  cards  in  the  exact  order  indicated  below  and  then  placed  immedi- 
ately after  the  *  DATA  card  at  the  end  of  the  source  program. 
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CARD  ORDER        DESCRIPTION 

1  Case  number — any  number  up  to  10  places  may  be 
assigned,  and  it  should  be  punched  within  the 
first  10  places  on  the  card.  Blanks  will  be  read 
as  zeros;  use  no  decimal  points. 

2  If  there  is  another  case  with  more  data  following 
this  set  of  data,  punch  +1  in  spaces  9  and  10  on 
this  card;  if  this  is  the  last  case  for  this  run, 
punch  -1  in  places  9  and  10. 

3  In  the  first  U  places  of  this  card  punch  the  num- 
ber of  weights  in  the  mathematical  model.  Within 
spaces  5  through  8  punch  the  number  of  weights 

in  the  mathematical  model  minus  1.   Use  no  decimals, 

h  Punch  within  the  first  10  spaces  the  decimal  form 

of  the  time  interval  to  be  used  in  seconds. 

The  next  M  cards     M  is  the  number  of  weights  in  the  mathematical 

model.  One  card  is  to  be  used  for  each  weight. 
Punch  the  value  in  pounds  of  each  weight  in  the 
first  10  spaces  of  each  card.  The  first  card  is 
for  weight  1  (the  ram)  and  then  consecutively 
down  to  the  pile  tip.  Use  a  decimal  point  and 
the  number  can  be  located  anywhere  within  the 
first  10  spaces. 

The  next  (M-l)  cards  Punch  within  the  first  10  spaces  of  each  of  the 

cards  the  spring  constants  beginning  with  the  cap- 
block  and  extending  to  the  top  of  the  pile  tip 
section  of  the  model.  The  values  are  to  be  in 
pounds  per  inch.  Use  a  decimal  point  and  the  num- 
ber can  be  located  anywhere  within  the  first  10 
spaces. 

The  next  M  cards     Punch  within  the  first  10  spaces  of  each  card  the 

soil  spring  coefficients  beginning  with  block  1. 
At  least  the  first  2  or  3  will  be  zero  normally. 
Use  a  decimal  point  and  the  number  in  pounds  per 
inch  can  be  placed  anywhere  within  the  first  10 
spaces. 

The  next  card      Punch  with  a  decimal  point  anywhere  within  the 

first  10  spaces  on  this  card  the  value  of  the 
capblock  coefficient  of  restitution. 

The  next  card      Punch  with  a  decimal  point  anywhere  with  in  the 

first  10  spaces  on  this  card  the  value  of  the 
coefficient  of  restitution  of  the  head  packing 
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in  the  case  of  a  concrete  pile.  For  all  other 
cases  without  head  packing  punch  1.00. 

The  next  card      Punch  with  a  decimal  point  anywhere  within  the 

first  10  spaces  the  velocity  of  the  ram  at  the 
instant  of  impact  in  feet  per  second. 

The  next  card      Punch  with  a  decimal  point  anywhere  within  the 

first  10  spaces  the  value  of  ground  quake  in  inches. 

The  next  card      Punch  with  a  decimal  point  anywhere  within  the 

first  10  spaces  the  value  of  point  damping  to  be 
used. 

The  next  card      Punch  with  a  decimal  point  anywhere  within  the 

first  10  spaces  the  value  of  side  damping  to  be 
used. 

The  next  card      If  the  pile  cap  is  allowed  to  transmit  tension, 

enter  a  +1  in  spaces  9  and  10;  if  not,  enter  a 
-1  there.  Use  no  decimals. 

The  next  card      If  side  friction  is  present,  enter  a  +1  in  spaces 

9  and  10;.  if  not,  enter  a  -1  there.  Use  no 
decimals. 

Running  the  Program:  The  data  deck  of  cards  just  described 
is  placed  immediately  after  the  source  program  deck  and  is  ready  to  be 
turned  over  to  a  machine  operator  for  processing. 

Results:  The  computer  will  clearly  write  out  in  the  output 
the  input  data;  the  number  of  each  time  interval;  and  identify  the  dis- 
placements, velocities,  permanent  sets,  forces,  and  ground  resistances 
acting  on  each  block.  The  results  are  printed  in  decimal  form  with  a 
multiplier  notation.  For  example,  0.2207  E06  means  0.2207  x  10  or 
220,700. 

The  Program;  The  complete  program  is  presented  so  that  any 
who  wish  to  use  it  may  do  so  freely.  Each  line  of  the  program  is  punched 
on  a  separate  card  and  must  be  assembled  in  the  exact  order  shown.  Any- 
one preparing  data  cards  from  this  program  should  verify  painstakingly 
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that  no  mistakes  have  crept  in  during  punching,  for  the  slightest  dis- 
crepancy can  lead  to  totally  erroneous  results.  For  that  reason  those 
preparing  cards  from  the  program  given  here  are  cautioned  to  try  it 
out  with  test  cases  before  relying  on  the  results. 

The  notation  used  in  the  program  is  the  same  as  was  used  in  the 
derivation  of  the  basic  equations  except  as  dictated  by  Fortran  require- 
ments. The  more  important  variations  in  terminology  are: 
S     instead  of  K   for  spring  constants 
SP    instead  of  K*  for  soil  spring  coefficients 

RES1  instead  of  e-^  for  coefficient  of  restitution  for  the 

capblock 

RES2  instead  of  e?  for  coefficient  of  restitution  for 

head  packing 

Z    instead  of  J   for  point  damping  factor 

ZP    instead  of  J'  for  side  damping  factor 
Additionally,  the  distinction  between  the  present  time  interval,  the 
last  time  interval,  and  the  second  previous  time  interval  is  made  by 
use  of  the  letters  L,  M,  and  S  respectively.  For  example,  in  the  case 
of  velocities, 

VL    is  used  instead  of  V 

VM    is  used  instead  of  v 

VS    is  used  instead  of  V* 
Other  variables  introduced  in  the  program  are  defined  as  they  occur 
within  the  program  or  have  been  explained  under  preparation  of  data 
above. 
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*ID   REESE        C2165   CIVIL  ENGR   G  WAVE  EQ   5  $0 

*  XEQ 

#  LABEL 

C     WAVE  EQUATION 

DIMENSION  DS(IOO),  DM(IOO),  DL(100),VS(100),  VM(IOO),  VL(IOO) 

DIMENSION  DPSM(lOO),  DPSL(lOO) 

DIMENSION  S(IOO),  W(IOO),  DP(300),  R(lOO),  SP(100),CS(100) 

DIMENSION  CM(lOO),  CL(100),FS(100),  FM(IOO),  FL(100),DPSS(100) 
207   DO  500  I  =  1,100 

DS(I)  =  0.0 

DM(I)  =0.0 

DL(I)  =0.0 

VS(I)  =  0.0 

VM(I)  =  0.0 

VL(I)  =  0.0 

DPSM(I)  =0.0 

DPSL(I)  =0.0 

S(I)  =  0.0 

W(I)  =  0.0 

R(l)  =  0.0 

SP(I)  =  0.0 

CS(I)  =  0.0 

CM(I)  =0.0 

CL(I)  =0.0 

FS(I)  =0.0 

FM(I)  =  0.0 

FL(I)  =  0.0 

500  DPSS(I)  =0.0 
DO  501  I  =  1,300 

501  DP(I)  =  0.0 

READ  INPUT  TAPE  5,27,NCASE,ITURN 

27  FORMAT  (I10/T10) 

0  READ  INPUT  TAPE  5,    70,M,TM1,T,  (W(I),I=1,M),  (S(l),I=l,IMl) 

1  ,(SP(I),I=1,M),  RES1,  RES2,  V,Q,Z,ZP,CAP,RUB 
70       FORMAT  (2LU/(F10.U)) 

WRITE  OUTPUT  TAPE  6,30,NCASE 

30  FORMAT  (12H  CASE  NUMBER  15) 

WRITE  OUTPUT  TAPE  6,28,  (I,W(I),S(l),SP(l),I=l,M) 

28  FORMAT  (5H  DATA//58H         M  W  S 
1                        SP/(IU,3F20.2)) 

•       WRITE  OUTPUT  TAPE  6,29,RES1,RES2,ZP,Z,Q,V 

29  FORMAT  (36H  COEF.   OF  REST.   OF  CAPBLOCK   (RES1)   =  FU.2/ 

1  36H  COEF.   OF  REST.   OF  PILECAP  (RES2)  =     FU.2/ 

2  36H  SIDE  DAMPING  FACTOR  (ZP)  =  FU.2/ 

3  36H  POINT  DAMPING  FACTOR  (Z)  »  FU.2/ 
U  36H  GROUND  QUAKE  (Q)  =  FU.2/ 
5  36H  INITIAL  RAM  VELOCITY  (V)  =                    F6.2) 

WRITE  OUTPUT  TAPE  6, 31, T, CAP, RUB 

31  FORMAT  (36H  TIME  INTERVAL  =  F8.6/ 
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1  36H  CAP  = 

2  36H  RUB  = 

WRITE  OUTPUT  TAPE  6,1 
FORMAT  (8L|H 
1  FL 

V-*12.*-T 


F6.2/ 
F6.2) 


M 


R//) 


DL 


VL 


DP 


DS(1J 

VS(1)  =    V  +(-DS(l)-*S(l))*T*32.2/W(D 

VS(2)  =  (DS(l)*S(l))*T*32.2/W(2) 

DM(1)  =  DS(1)  +  VS(l)*12.-*T 

DM(2)  =  VS(2)*12.*-T 

CM(2)  =  DM(2) 

VM(l)  =  VS(l)  -((DM(l)~DM(2))ttS(l))*T*32.2/W(l) 

VM(2)  =  VS(2)+((DM(l)-DM(2))-x-S(l)-(DM(2)-DM(3))-"-S(2))-x-T-32.2M2) 

VM(3)  ■  VS(3)  +((DM(2)-DM(3))*S{2)-R3)*T*32.2/W(3) 

CM(1)  =  DM(1)  -  DM(2) 

CHECK  «■  2.-*  V 

WRITE  OUTPUT  TAPE  6,301,DS(l),  VS(l) 

301  FORMAT  (8H  N  =   1,6H     1  2E1$.U//) 

WRITE  OUTPUT  TAPE  6,302, DM(l),VM(l),DM(2),VM(2) 

302  FORMAT  (OH  N  =   2,6H     1  2El5.U/lljH  2  2E15.U//) 
DO  101  N  =  3,  300 

WRITE  OUTPUT  TAPE  6,303,N 

303  FORMAT  (£H  =  13) 
MLESS1  =  M-l 

DO  130  I  =  1  M 

130  DL(I)  =  DM(I)  +  VM(I)*12.*T 
DO  131  I  =  1,MLESS1 

131  CL(I)  =  DL(I)  -  DL(I+1) 
13   IF(DL(M))  12,12,lU 

12   DE  =  0.0 

GO  TO  2k 
1U   IF(DL(M)  -  Q)  16,16,18 
16   DE  -  0.0 

GO  TO  2k 
18   VALUE  =  DL(M)  -  Q 

IF  (VALUE-DE)  22,22,20 
20   DE  -  VALUE 

GO  TO  2k 
22   DE  =  DE 
2k       DP(N)  =  DE 

R(M)  =  (DL(M)-DP(N))*SP(M)*(1.+Z#VM(M)) 

IF  (DP(N))  25,25,26 
26   IF  (DP(N)  -  DP(N-l))  197,197,25 
2$       CONTINUE 
3   VAL  =  S(l)*CL(l) 

IF(CL(l)-CM(l))  5,5, li 
k       FL(1)  -  VAL 

GO  TO  33 
5   CL(1)  =  CM(1) 

FL(1)«-VAL/(RES1**2 )-(l,/(RESI**2 )-l.  )*3(l)aCL(l) 
33   VAL2  =  S(2)*CL(2) 
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IF(CL(2)-CM(2))  35,35,3U 

3U 

FL(2)   =  VAL2 

GO  TO  37 

35 

CL(2)   =  CM(2) 

FL(2)   =  VAL2/(RES2-:(-x-2)-(l./(RES2-::-x-2)-l.)-«-S(2KL(2) 

37 

IF (CAP)  38,36,36 

38 

IF(FL(2))  39,36,36 

39 

FL(2)  =  0.0 

36 

CONTINUE 

7 

IF(FL(1))  8,133,133 

8 

FL(1)   =  0.0 

133 

DO  132  I  =  3.MLESS1 

132 

FL(I)  =  CL(I)*S(I) 

U2 

IF(RUB)U9,U9,UU 

ill; 

DO  U8  I  =  3,MLESS1 

DPSL(I)  =  DPSM(I) 

CHANGE  =  DL(I)  -  Q 

SUM  +  DL(I)  +  Q 

IF(DPSL(I)-CHANGE)l45,ii6,U6 

l£ 

DPSL(I)  =  CHANGE 

U6 

IF(DPSL(I)-SUM)U8,l|8,ii7 

U7 

DPSL(I)  =  SUM 

U8 

R(I)  =  (DL(I)  -  DPSL(I))*SP(lMl.+ZP*VM(I)) 

U9 

DO  50  I  =  1,M 

VL(I)  =  VM(I)   +  (FL(I-l)  -  FL(I)  -  R(I))*T*32.2/W(l) 

50 

CONTINUE 

DO  100  I  =  1,M 

703 

WRITE  OUTPUT  TAPE  6,30U,I,DL(I), VL(l),DPSL(I),FL(I),R(l) 

301; 

FORMAT  (11H                       I3,5E15.U) 

100 

CONTINUE 

DO  75  K  =  1,M 

X  =  1. 

IF  (VL(K))           75,75,800 

75 

CONTINUE 

GO  TO  199 

800 

CONTINUE 

IF(VL(2)  -  CHECK)99,99,190 

99 

IF(VL(M)-CHECK)102,102,19U 

102 

DO  98  K  -  1,M 

DS(K)   =  DM(K) 

DM(K)   =  DL(K) 

DL(K)   =  0.0 

VS(K)  =  VM(K) 

VM(K)  =  VL(K) 

VL(K)   =  0.0 

CS(K)  =  CM(K) 

CM(K)   =  CL(K) 

CL(K)  =  0.0 

FS(K)   =  FM(K) 

FM(K)   =  FL(K) 

FL(K)  =  0.0 
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DPSS(K)  =  DPSM(K) 

DPSM(K)  =  DPSL(K) 
98   DPSL(K)  =0.0 

WRITE  OUTPUT  TAPE  6, 30$,  DP  (N) 
30^   FORMAT  (UjH  KlS.h) 

101   CONTINUE 

GO  TO  196 

190  WRITE  OUTPUT  TAPE  6,1?1,N 

191  0  FORMAT  (6I4H  VELOCITY  OF  PILE  CAP  EXCEEDED  TWICE  THE  RAM  VELOCITY  W 
1  HEN  N  WAS  13) 

GO  TO  196 
19U   WRITE  OUTPUT  TAPE  6,1#,N 

195  0  FORMAT  (6I4H  VELOCITY  OF  PILE  TIP  EXCEEDED  TWICE  THE  RAM  VELOCITY  W 
1  HEN  N  WAS  13) 

GO  TO  196 

197  WRITE  OUTPUT  TAPE  6,198,N 

198  FORMAT  (26H  DP  BECAME  CONSTANT  AT  N  =  13) 
GO  TO  196 

199  WRITE  OUTPUT  TAPE  6,200,N 

200  FORMAT   ($2H  ALL  VL  WERE  SIMULTANEOUSLY  NEGATIVE  OR  ZERO  AT  N  =  IU) 

196  IF(ITURN)205,205,207 
20$       CALL  EXIT 

END 
*  DATA 
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APPENDIX  B 
DETAILS  OF  THE  "VARY  Ru  PROGFJM "  WHICH  USES  THE  WAVE  EQUATION 

General ♦  This  program  is  identical  to  the  BASIC  PROGRAM  except 
that  it  provides  the  capability  of  varying  certain  input  data  in  order 
to  obtain  wave  equation  solutions  for  a  range  of  input  variables.  The 
program  provides  for  a  maximum  of  500  time  cycles,  instead  of  300,  to 
allow  for  the  use  of  smaller  time  increments  and  a  greater  number  of 
pile  segments.  Otherwise,  the  termination  features  are  the  same  as  for 
the  BASIC  PROGRAM. 

Program  Features.  The  program  is  designed  with  the  following 
features  that  are  different  from  or  in  addition  to  those  of  the  basic 
program. 

1.  The  ultimate  resistance  (Ru)  is  varied  over  a  preselected 
range  of  values  in  any  increment  desired.  The  set  is  computed  for 
each  value  of  Ru.  The  initial  value  of  Ru,  the  increment  of  increase 
of  Ru  and  the  number  of  times  the  increment  is  applied  are  provided  as 
input  data. 

2.  The  effective  length  of  embedment  of  the  pile  can  be  varied 
and  is  determined  by  input  data.  This  is  the  length  of  the  pile  over 
which  the  resistance  to  driving  acts. 

3.  Five  patterns  of  resistance  distribution  between  point  bear- 
ing and  side  friction  can  be  applied  for  computation.  These  patterns 
vary  from  full  point  bearing  to  full  side  friction  with  intermediate 
distributions  of  75,  50,  and  25  per  cent  point  bearing  with  the  balance 
distributed  as  side  friction. 
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U.  The  method  of  distributing  side  friction  can  be  rectangu- 
lar or  triangular  as  desired.  Input  data  determine  the  selection  of 
the  method. 

5.  The  capability  of  recycling  the  entire  program  up  to 
three  times  with  the  new  values  of  any  variable  is  provided.  This 
feature  was  used  to  provide  new  values  of  Q  and  J,  although  other 
variables  could  be  used.  This  capability  is  provided  with  a  c  ontrol 
statement,  and  the  new  values  of  the  variables  to  be  changed  must  be 
put  into  the  program  itself,  a  somewhat  awkward  arrangement. 

6.  Any  number  of  sets  of  data  can  be  used  in  a  single  run, 
as  provided  for  in  the  input  data.  This  is  also  a  feature  of  the 
basic  program. 

Input  Data.  The  input  data  is  the  same  as  for  the  basic 

program  except  for  the  additional  data  required. 

CARD  ORDER  DESCRIPTION 

2a  After  the  second  card  of  the  basic  program, 

areas  of  the  pile  head,  center  and  tip  in  square 
inches  are  punched  on  one  card  in  the  first 
three  ten-column  fields.  Values  are  punched  with 
a  decimal  point. 

After  the  last  card    In  the  first  four  columns  punch  1  to  indicate 
of  the  basic  program   triangular  distribution  of  resistance  on  the  side 

or  -1  to  indicate  rectangular  distribution. 

Next  card        In  the  first  10  columns  punch  the  initial  value 

of  Ru  in  pounds,  and  in  the  second  10-column 
field  the  increment  of  Ru  in  pounds.  A  decimal 
point  must  be  used.  These  are  variables  RU1 
and  ADD. 

Next  card        In  the  first  four  columns  punch  the  number  rep- 
resenting the  number  of  times  the  increment  of 
Ru  is  to  be  applied;  this  establishes  themaximum 
value  of  Ru.  In  the  second  four  columns  (5  -  8) 
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punch  the  number  of  blocks  in  the  mathematical 
model  to  which  ground  resistance  is  not  to  be 
applied;  this  establishes  the  effective  embedded 
length  of  the  pile.  These  variables  are  ICYCLE 
and  NOLOAD  and  are  punched  without  a  sign  or 
decimal  point. 

Next  card        A  number  of  1  through  $   is  punched  in  each  of  the 

first  three  four-column  fields  for  the  variables 
LA,  LB,  and  LC  respectively.  These  numbers  pro- 
vide the  values  for  the  indexing  parameters  of 
a  DO  statement  that  determines  the  distribution 
of  resistance  between  point  and  side. 

Next  and  last  card    A  number  with  a  value  between  1  through  k   is 

punched  in  the  first  four  columns.  This  number 
controls  the  number  of  times  the  program  is  to 
be  repeated  with  new  variables.  No  sign  or  deci- 
mal is  used.  This  variable  is  called  IRES.  If 
the  number  1  is  used  the  program  will  cycle  the 
full  three  times,  and  if  h   is  used  it  will  stop 
after  all  computations  using  the  initial  data 
are  made. 

Output  Data.  In  addition  to  listing  all  input  data,  the  follow- 
ing output  is  provided: 

1.  Value  of  Ru  in  pounds  for  the  particular  computation. 

2.  Set  in  blows  per  inch. 

3.  Compression  and  tension  stresses  at  the  head,  center  and 
tip  of  the  pile,  in  pounds  per  square  inch. 

U.  The  number  of  time  cycles  run. 

Use  of  the  Program.  This  program  is  designed  to  provide  the 
data  necessary  to  plot  Ultimate  Resistance  vs.  Set  curves  for  various 
values  of  Q,  J, and  J»,  and  for  one  or  more  of  the  five  different  dis- 
tributions of  resistance  between  point  and  side.  The  values  of  Q,  J, 
and  J'  are  initially  provided  in  the  input  data.  If  additional  values 
are  desired  for  computations,  they  are  provided  through  the  use  of  the 
"computed-go-to"  statement  following  statement  number  600.  The 
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particular  values  must  be  placed  in  the  program  following  statements 
700,  701,  and  702.  This  feature  could  easily  be.  expanded  if  desired. 
If  it  is  desired  to  use  all  five  distributions  of  resistance  between 
point  and  side,  the  values  of  LA,  LB,  and  LC  must  be  1,  £,  and  1  re- 
spectively. When  the  index  variable  of  the  DO  statement  (five  state- 
ments after  777)  takes  different  values,  the  distribution  of  resistance 
is  as  indicated  below: 

Value  Distribution 

1  All  point  bearing 

2  7$  per  cent  point  bearing,  2$   per  cent  friction 

3  50  per  cent  point  bearing,  $0  per  cent  friction 
h  2$  per  cent  point  bearing,  75  per  cent  friction 
5        All  side  friction 

By  assigning  proper  values  to  LA,  LB,  and  LC,  the  desired  distributions 
are  obtained  in  accordance  with  the  operation  of  a  DO  statement  as  ex- 
plained in  reference  (21). 
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*ID   REESE        C2l6$   CIVIL  ENGR   G  WAVE  EQ  1$     250  VARY  RU 

*  XEQ 

*  LABEL 

C  WAVE  EQUATION  VARY  RU  101 

DIMENSION  DS(lOO),   DM(lOO),   DL(100),VS(100),   VM(lOO),   VL(lOO) 

DIMENSION  DPSM(lOO),   DPSL(lOO) 

DIMENSION  S(100)     W(100),   DP(500)     R(100),   SP(100),CS(100) 

DIMENSION  CM(IOO),   CL(100),FS(lOO),  FM(lOO),   FL(lOO),DPSS(100) 
207       READ  INPUT  TAPE  5,27,NCASE,ITURN 

27  FORMAT  (IIO/IIO) 

READ  INPUT  TAPE  5,71,AREA1,AREA2,AREA3 
7.1   FORMAT  (3F10.ii) 

OREAD  INPUT  TAPE  $,  70,  M,IML,T,  (W(l),I=l,M),  (S(l),I=l,IMl) 
1  ,(SP(I),I=1 M),  RESl,  RES2,  V,Q,Z,ZP,CAP,RUB 
70   FORMAT  (2lU/(F10.U)) 

READ  INPUT  TAPE  5,60U,MUD 

6Dh      FORMAT(lU) 

READ  INPUT  TAPE  5,72,RU1,ADD,ICYCLE,N0L0AD 

72  FORMAT  (2P10.U/2IU) 

READ  INPUT  TAPE  5, 73, LA, LB, LC, IRES 

73  FORMAT  (3IU/IU) 

WRITE  OUTPUT  TAPE  6,30,NCASE 

30  FORMAT  (12H  CASE  NUMBER  1$) 

WRITE  OUTPUT  TAPE  6,28,  (I,W(l),S(I),SP(l),I=l,M) 

28  FORMAT  (5H  DATA//58H    M  W  S 
1         SP/(IU,3F20.2)) 

WRITE  OUTPUT  TAPE  6,29,RES1,RES2,ZP,Z,Q,V 

29  FORMAT  (36H  COEF.   OF  REST.   OF  CAPBLOCK   (RESl)  =  FU.2/ 

1  36H  COEF.   OF  REST.   OF  PILECAP   (RES2)  =     FU.2/ 

2  36H  SIDE  DAMPING  FACTOR   (ZP)  =  FU.2/ 

3  36H  POINT  DAMPING  FACTOR  (Z)  =  FU.2/ 
U  36H  GROUND  QUAKE  (Q)  =  FU.2/ 
5  36H  INITIAL  RAM  VELOCITY  (V)  =  F6.2) 

WRITE  OUTPUT  TAPE  6,60$, MUD 
60$      FORMAT  (6H  MUD  ■  16) 

WRITE  OUTPUT  TAPE  6, 31, T, CAP, RUB 

31  FORMAT   (36H  TIME  INTERVAL  =  F8.6/ 

1  36H  CAP  -  F6.2/ 

2  36H  RUB  ■  F6.2//) 
WRITE  OUTPUT  TAPE  6,608,RUl,ADD,ICYCLE,NOLOAD 

608  FORMAT  ($H  RU  =  E20.U,10H  ADD  =  E20.U/9H  ICYCLE  -  16, 
1  13H           NOLOAD  =  16) 

WRITE  OUTPUT  TAPE  6,609,LA,LB,LC,IRES 

609  FORMAT  ($H  LA  =  lU,10H  LB  =  IU,10H  LC   =  lU,10H  LRES 

1     IU) 

CHECK  ■  2.*  V 
777   CONTINUE 
RU  =  RU1 
DO  600  IAM  -  1, ICYCLE 


■ 


'■ 
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ITEM  =  M  -  NOLO AD 
AITEM  =  ITEM 
DO  602  J  =  LA,LB,LC 
FMAX1  =0.0 
FMAX2  =0.0 
FMAX3  =0.0 
FTEN1  =0.0 
FTEN2  =0.0 
FTEN3  =0.0 
DO  500  I  -  1,100 
DS(I)   =  0.0 
DM(I)  =  0.0 
DL(I)   =  0.0 
VS(I)  =  0.0 
VM(I)  =  0.0 
VL(l)   =  0.0 
DPSM(I)   =0.0 
DPSL(I)   =  0.0 
R(I)   =  0.0 
SP(I)  =  0.0 
CS(I)   =  0.0 
CM(I)   =  0.0 
CL(I)   =  0.0 
FS(I)   =  0.0 
FM(I)   =  0.0 
FL(I)  =  0.0 

500  DPSS(I)   =0.0 
DO  501  I  =1,500 

501  DP(I)   =  0.0 
DS(1)   =  V*12.*T 

VS(1)  =  V  +(-DS(l)*S(l))*T*32.2/W(l) 

VS(2)  =  (DS(1)*S(1))*T*32.2/W(2) 

DM(1)  =  DS(1)  +  VS(i)*12.*T 

DM(2)  =  VS(2)*12.*T 

CM(2)  =  DM(2) 

VM(1)  =  VS(1)  -((DMCl)-DM(2))*S(l))*T*32.2Ml) 

VM(2)  =  VS(2)+((DM(1)-DM(2))^S(1)-(DM(2)-DM(3))*S(2))J-T^32.2M2) 

VM(3)  =  VS(3)+((DM(2)-DM(3))-"-S(2)-R(3))J"-T-^32.2/W(3) 

CM(1)  =  DM(1)  -  DM(2) 

AJ   =  J 

PART  =   (5.-AJ)/U. 

SP(M)  =  PART#RU/Q 

SIDE  =  RU  -  PART*RU 

LOAD  =  NOLOAD  +  1 

IF(MUD)650,65l,652 
650       DO  601  N  =  LOAD,IML 
601       SP(N)  =  SIDE/(AITEM*Q) 

SP(M)   =  SP(M)  +  SIDE/(AITEM*Q) 

GO  TO  651 
652       DO  603  N  =  3,IM1 

NT  -  N-2 
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NT1  =  2-*NT  -  1 

IM2  =  M-2 

ANT1  =  NT1 

AIM2  =  IM2 
603   SP(N)  ■  SIDE*ANTl/((AIM2**2)*Q) 

SP(M)  =  SP(M)  +  SIDE*(2.*AIM2-l.)/((AIM2-^2)-MQ) 
651   CONTINUE 

DO  101  N  =  3,500 

MLESS1  =  M-l 

DO  130  I  =  1  M 

130  DL(I)  =  DM(I)  +  VM(I)*12.*T 
DO  131  I  ■  1.MLESS1 

131  CL(I)  =  DL(I)  -  DL(I+1) 
13   IF(DL(M))  12,12,1k 

12   DE  =  0.0 

GO  TO  2k 
Ik      IF(DL(M)  -  Q)  16,16,18 
16   DE  =  0.0 

GO  TO  2k 
18   VALUE  =  DL(M)  -  Q 

IF  (VALUE-DE)  22,22,20 
20   DE  *  VALUE 

GO  TO  2k 
22   DE  =  DE 
2k       DP(N)  =  DE 

R(M)  =  (DL(M)-DP(N))*SP(M)*(l.+Z*VM(M)) 

IF  (DP(N))  25,25,26 
126   IF  (DP(N)  -  DP(N-l))  197,197,25 
25   CONTINUE 
3   VAL  =  S(1)*CL(1) 

IF(CL(1)-CM(1))  5,5,U 
k       FL(1)  =  VAL 

GO  TO  33 
5   CL(1)  =  CM(1) 

FL(l)=VAL/(RESl*-*2  )-(l./(RESl**2)-l.  )*S(l)*CL(l) 
33   VAL2  =  S(2)*CL(2) 

IF(CL(2)-CM(2))  35,35,3U 
3il-   FL(2)  =  VAL2 

GO  TO  37 

35  CL(2)  =  CM(2) 

FL(2)  =  VAL2/(HES2**2)-(l./(HES2«*2)-l.)*S(2)iiCL(2) 

37  IF(CAP)  38,36,36 

38  IF(FL(2))  39,36,36 

39  FL(2)  =  0.0 

36  CONTINUE 

7  IF(FL(1))  8,133,133 

8  FL(1)  =  0.0 

133   DO  132  I  =  3.MLESS1 

132  FL(I)  =  CL(I)*S(I) 
LF(FL(2))213,211,210 

213   IF(FTEN1-FL(2))211,211,2LU 
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2lU   FTEN1  =  FL(2) 
GO  TO  211 

210  IF(FMAXl-FL(2))2l5,211,211 
205   FMAX1  =  FL(2) 

211  CONTINUE 

MC  =  M/2  +  1 
IF(FL(MC))2l6,217,2l8 

216  IF(FTEN2-FL(MC))217,217,219 

219  FTEN2  =  FL(MC) 
GO  TO  217 

218   IF(FMAX2-FL(MC))220,217,217 

220  FMAX2  -  FL(MC) 

217  CONTINUE 
IF(FL(IMl))221,222,22U 

221  IF(FTEN3«-FL(na))222,222,223 
223   FTEN3  =  FL(IMl) 

GO  TO  222 
22U   IF(FMAX3-FL(IML))225,222,222 
225   FMAX3  =  FL(IKL) 

222  CONTINUE 

k2       IF(RUB)U9,U9,UU 
hk       DO  U8  I  =  3,MLESS1 

DPSL(I)   =  DPSM(I) 

CHANGE  =  DL(I)  -  Q 

SUM  -  DL(I)   +  Q 

IF(DPSL(I)-CHANGE)li5,U6,U6 
b$       DPSL(I)   =  CHANGE 
1*6       IF(DPSL(I)-SUM)U8,U8,U7 
kl       DPSL(I)   =  SUM 

U8       R(I)  =  (DL(I)  -  DPSL(I))*SP(l)*(l.+ZP*VM(l)) 
U9       DO  50  I  =  1,M 

VL(I)  -  VM(I)  +  (FL(I-l)  -  FL(I)  -  R(I))*T*32.2MD 
50       CONTINUE 

DO  75  K  =  1,M 

X  =  1. 

IF  (VL(K))  75,75,800 

75       CONTINUE 

GO  TO  199 
800       CONTINUE 

IF(VL(2)  -  CHECK)99, 99, 190 
99       IF (VL(M)-CHECK)102, 102,191 
102       DO  98  K  =  1,M 


DS(K 
DM(K 
DL(K 
VS(K 
VM(K 
VL(K 
CS(K 
CM(K 
CL(K 
FS(K 


=  DM(K) 
-  DL(K) 
=  0.0 
=  VM(K) 
=  VL(K) 
=  0.0 
=  CM(K) 
=  CL(K) 
=  0.0 
=  FM(K) 
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VM(K)   =  FL(K) 
FL(K)  =  0.0 
DPSS(K)   =  IPSM(K) 
DPSM(K)   -  DPSL(K) 
98       DPSL(K)  =  0.0 
101       CONTINUE 
GO  TO  196 

190  WRITE  OUTPUT  TAPE  6,191,N 

191  0~  FORMAT  (6I4H  VELOCITY  OF  PILE  CAP  EXCEEDED  TWICE  THE  RAM  VELOCITY  W 
1  HEN  N  WAS  13) 

GO  TO  196 
19U       WRITE  OUTPUT  TAPE  6,195,N 

195  0  FORMAT  (6I4H  VELOCITY  OF  PILE  TIP  EXCEEDED  TWICE  THE  RAM  VELOCITY  W 
1  HEN  N  WAS  113) 

GO  TO  196 

197  WRITE  OUTPUT  TAPE  6,1?8,N 

198  FORMAT  (26H  DP  BECAME  CONSTANT  AT  N  =  13) 
GO  TO  196 

199  WRITE  OUTPUT  TAPE  6,200,N 

200  FORMAT   (52H  ALL  VL  WERE  SIMULTANEOUSLY  NEGATIVE  OR  ZERO  AT  N  =  Ik) 

196  CONTINUE 

WRITE  OUTPUT  TAPE  6,607,RU 
607       FORMAT   (27H  ULTIMATE  RESISTANCE   (RU)   =  E15.U) 

WRITE  OUTPUT  TAPE  6,330,Q,W(1),V 
330       FORMAT  (1*11  Q  =  E10.2,7H  W(l)  =  El5.U.i|H  V  =  E15.U) 

WRITE  OUTPUT  TAPE  6,606,  (SP(L),L=1,M) 
606       FORMAT  (17H  NEW  VALUES  OF  SP/(8El5.U)) 

WRITE  OUTPUT  TAPE  6,305,DP(N) 

305  FORMAT  (7H  SET  =  E10.U) 
BLOW  =  l./DP(N) 

WRITE  OUTPUT  TAPE  6,805,BLOW 
80$       FORMAT  (66H 

1  BLOWS/EN=  E10.3) 

FTEN1  =  FTENi/aREAI 

FTEN2  =  FTEN2/AREA2 

FTEN3  -  FTEN3/AREA3 

WRITE  OUTPUT  TAPE  6,307,FTEN1,FTEN2,FTEN3 
307  0  F0RMT(30H  MAXIMUM  TENSION  AT  HEAD  =         E15.U/ 

1  30H  MAXIMUM  TENSION  AT  CENTER  =     E15.U/ 

2  30H  MAXIMUM  TENSION  AT  TIP  =  E15.U) 
FMAX1  =  FMAXl/AREAl 

FMAX2  =  FMAX2/AREA2 
FMAX3  ■  FMAX3/AREA3 
VJRITE  OUTPUT  TAPE  6,306,FMAX1,FMAX2,FMAX3 

306  0  FORM/IT  (26H  MAXIMUM  FORCE  AT  HEAD  =       E15.U/ 

1  26H  MAXIMUM  FORCE  AT  CENTER  -  El5.lt/ 

2  26H  MAXIMUM  FORCE  AT  TIP  =         E15.U//) 
602       CONTINUE 

RU  =  RU  +  ADD 
600       CONTINUE 
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GO  TO  (700, 701, 702, 703), IRES 

700  CONTINUE 
Z  =  0.20 
ZP  =  0.067 
IRES  =  2 
GO  TO  777 

701  CONTINUE 
Q  =  0.20 
IRES  =  3 
GO  TO  777 

702  CONTINUE 
Z  =  0.15 
ZP  =  0.05 
IRES  =  k 
GO  TO  777 

703  CONTINUE 
IF(ITURN)205, 205,207 

205   CALL  EXIT 

END 
-x-     DATA 
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APPENDIX  C 

DETAT  LS  OF  THE  "RESEARCHER  COMPUTER  PROGRAM" 
WHICH  USES  THE  WAVE  EQUATION 


General:  This  program  is  written  in  Fortran  language  for  the  IBM 
7090  computer,  but  it  can  be  used  equally  well  with  the  709  or  709U 
computer.  It  has  as  its  objective  the  systematic  production  of  data 
from  which  set  vs.  ultimate  ground  resistance  curves  can  be  plotted 
for  any  case.  Its  chief  advantage  lies  in  the  fact  that  for  every  set 
of  data  run,  the  program  automatically  begins  with  a  value  of  ground 
resistance  of  20  tons  and  finds  the  set  associated  with  all  the  load 
carried  by  the  pile  tip,  1$%   of  load  carried  by  the  pile  tip,  $0%   of 
load  carried  by  the  pile  tip,  2$%   of  load  carried  by  the  pile  tip,  and 
finally  all  the  load  carried  by  friction  along  the  sides  of  the  pile. 
It  then  assumes  an  increase  of  ground  resistance  of  20  tons  and  repeats 
the  calculations.  It  automatically  assigns  new  values  of  ground  re- 
sistance in  this  manner  up  to  and  including  280  tons.  Essentially, 
then,  it  progresses  through  the  basic  program  five  times  at  each  value 
of  ultimate  ground  resistance,  or  70  times  for  each  set  of  data.  It 
also  writes  out  the  maximum  compression  and  tension  stresses  in  pounds 
per  square  inch  for  each  of  these  70  conditions  at  the  head,  mid-length, 
and  tip  of  the  pile.  The  side  friction  distribution  may  be  chosen  as 
rectangular  or  triangular  in  shape.  In  both  of  these  cases  the  top  10 
feet  of  the  pile  is  considered  free  of  frictional  resistance. 

Termination;  Termination  of  each  of  the  70  conditions  for  each 
set  of  data  is  programmed  exactly  like  the  basic  program.  The  reason 
for  each  termination  is  written  out,  and  the  computer  immediately  moves 
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to  the  next  set  of  conditions.  No  over-all  termination  between  sets 
of  data  should  occur. 

Stability:  Besides  the  checks  to  detect  instability  in  the 
calculations  as  described  under  the  basic  program,  a  further  check 
is  obtained  when  plotting  set  vs.  ultimate  ground  resistance.  This 
plot  should  also  produce  smooth  curves.  If  not,  a  set  of  data  for 
each  suspicious  point  can  be  prepared  and  run  with  the  basic  program 
to  obtain  a  detailed  analysis  so  that  evaluation  can  be  made  of  the 
validity  of  the  results.  Reduced  time  intervals  can  be  used  to  elim- 
inate the  instability,  or  smaller  pile  lengths  in  the  mathematical 
model  together  with  reduced  time  intervals  may  be  used. 

Preparation  of  Input  Data;  All  data  is  punched  on  standard 
cards  exactly  as  described  in  the  basic  program  except  that  two  addi- 
tional cards  are  required  for  this  program. 

The  first  extra  card  is  inserted  between  cards  2  and  3.  In 
the  first  10  spaces  of  this  new  card  are  punched,  with  a  decimal  point, 
the  cross  sectional  area  of  the  pile  at  its  top,  within  spaces  11 
through  20  is  punched  in  the  same  manner  the  pile  cross-sectional  area 
at  mid-length  of  the  pile,  and  similarly  the  area  of  the  pile  tip  in 
spaces  21  through  30. 

The  second  extra  card  is  included  as  the  very  last  card  of 
data.  If  side  distribution  is  rectangular,  punch  -1  in  spaces  3  and 
U;  if  side  distribution  is  to  be  triangular,  punch  +1  in  spaces  3  and  k. 

In  all  other  respects,  data  is  prepared  exactly  as  described 
under  instructions  for  the  "basic  computer  program." 


- 


■ 
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Running  the  program:  The  data  deck  of  cards  just  described 
is  placed  immediately  after  the  program  deck  (after  the  *  DATA  card) 
and  is  then  ready  to  be  processed  by  a  machine  operator. 

Results :  The  input  data  will  be  repeated  in  the  output 
information.  The  time  cycle  that  the  calculations  were  terminated 
for  the  first  condition  and  the  reason  therefor  will  first  appear. 
The  ultimate  ground  resistance;  the  soil  spring  coefficients;  set; 
and  the  maximum  tension  and  compression  at  the  pile  head,  mid-length, 
and  tip  will  be  written  out  for  each  condition.  The  results  are 
printed  in  decimal  form  with  a  multiplier  notation.  For  example, 
0.1^11  E01  means  0.1^11  x  101  =  1.511. 

The  Program:  The  complete  program  is  presented  so  that  any 
who  wish  to  use  it  may  do  so  freely.  Each  line  of  the  progran  is 
punched  on  a  separate  card  and  must  be  assembled  in  the  exact  order 
given  below.  Extreme  care  must  be  used  in  punching  to  prevent  errors 
from  creeping  in,  and  for  that  reason  it  is  recommended  that  test 
cases  be  run  to  check  the  program  before  relying  on  its  results. 

'The  notation  used  is  as  described  under  this  section  of 
instructions  for  the  "basic  computer  program." 
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x-ID   REESE        C2165   CIVIL  ENGR   G  WAVE  EQ  k   100     RESEARCHER 

*  XEQ 

*  NOBIN 

*  LABEL 

C     WAVE  EQUATION  VARY  RU 

DIMENSION  DS(lOO),  DM(IOO),  DL(100),VS(100),  VM(lOO),  VL(lOO) 

DIMENSION  DPSM(lOO),  DPSL(lOO) 

DIMENSION  S(100)  W(100),  DP(300),  R(100),  SP(100),CS(100) 

DIMENSION  CM(lOO),  CL(100),FS(lOO),  FM(IOO),  FL(lOO),DPSS(100) 
207   READ  INPUT  TAPE  5,27,NCASE,ITURN 

27  FORMAT  (110/110) 

READ  INPUT  TAPE  5,71,AREA1,AREA2,AREA3 
71       FORMAT (3F10.U) 

0  READ  INPUT  TAPE  5,    70,  M,IM1,T,  (W(I),I=1,M),  (S(l),I=l,IMl) 

1  ,(SP(I),I=1  M),  RES1,  RES2,   V,Q,Z,ZP,CAP,RUB 
70      FORMAT  (2IU/(F10.U)) 

READ  INPUT  TAPE  5,60U,MUD 
601;       FORMAT  (IU) 

WRITE  OUTPUT  TAPE  6,30,NCASE 

30  FORMAT  (12H  CASE  NUMBER  15) 

WRITE  OUTPUT  TAPE  6,28,  (I,W(l),S(I),SP(l),I=l,M) 

28  FORMAT  (5H  DATA//58H         M  W  S 
1                        SP/(Iii,3F20.2)) 

WRITE  OUTPUT  TAPE  6,29,RES1,RES2,ZP,Z,Q,V 

29  FORMAT  (36H  COEF.   OF  REST.   OF  CAPBLOCK  (RES1)  =  Flu 2/ 

1  36H  COEF.   OF  REST.   OF  PILECAP   (RES2)  *    Flu  2/ 

2  36H  SIDE  DAMPING  FACTOR  (ZP)  -  Flu  2/ 

3  36H  POINT  DAMPING  FACTOR  (Z)  -  Flu  2/ 
k  36H  GROUND  QUAKE  (Q)  -  FU.2/ 
5  36H  INITIAL  RAM  VELOCITY  (V)  =                    F6.2) 

WRITE  OUTPUT  TAPE  6  605,MUD 
605       FORMAT   (6H  MUD  =  16) 

WRITE  OUTPUT  TAPE  6, 31, T, CAP, RUB 

31  FORMAT  (36H  TIM  INTERVAL  «  F8.6/ 

1  36H  CAP  =  F6.2/ 

2  36H  RUB  =  F6.2//) 
CHECK  =  2.*  V 

RU  =  UOOOO. 

DO  600  I  =  1,1U 
ITEM  =  M-3 
AITEM  -  ITEM 
DO  602  J  =  1,5 
FMAX1  =  0.0 
FMAX2  =0.0 
FMAX3  =0.0 
FTEN1  =0.0 
FTEN2  =0.0 
FTEN3  =0.0 
DO  500  I  =  1,100 
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DS(I)  =  0.0 
DM(l)  =0.0 
DL(I)  =  0.0 
VS(I)  =  0.0 
VM(I)  -  0.0 
VL(I)  =  0.0 
DPSM(I)  =0.0 
DPSL(I)  =0.0 
R(I)  =  0.0 
SP(I)  =  0.0 
CS(I)  =  0.0 
CM(I)  =  0.0 
CL(I)  =  0.0 
FS(I)  =  0.0 
FM(I)  =  0.0 
FL(I)  =  0.0 

500  DPSS(I)  =0.0 

DO  ^01  I  -  1,300 

501  DP(I)  =0.0 
DS(1)  =  V*12.*T 

VS(1)  =  V  +(-DS(l)*S(l))*T*32.2/to(D 

VS(2)  =  (DS(1)*S(1))#T*32.2/W(2) 

DM(1)  -  DS(1)  +  VS(l)#12.*T 

DM(2)   =  VS(2)#12.*T 

CM(2)   =  DM(2) 

VM(1)   =  VS(1)  -((DM(l)-DM(2))*S(l))-*T*32.2/W(l) 

VM(2)  -  VS(2)+((DM(1)-DM(2)>S(1)-(DM(2)-DM(3))*S(2))«T*32.2M2) 

VM(3)  ■  VS(3)+((DM(2)-DM(3))*S(2)-R(3))*T*32.2/W(3) 

CM(1)   =  DM(1)   -  DM(2) 

AJ  =  J 

PART  =  (5.-AJ)/U. 

SP(M)   =  PART*RU/Q 

SIDE  =  RU  -  PART*RU 

IF(MUD)650,65l,6$2 

650  DO  601  N  =  U,IM1 

601       SP(N)  =  SIDE/(AITEM*Q) 

SP(M)   =  SP(M)  +  SIDE/(AITEMfcQ) 

GO  TO  651 
652      DO  603  N  =  3,IM1 

NT  =  N-2 

NT1  =  2#NT  -  1 

IM2  =  M-2 

ANT1  -  NT1 

AIM2  =  IM2 
603   SP(N)  =  SIDE*ANTl/((AIM2**2)*Q) 

SP(M)  =  SP(M)  +  SmE*(2.*AIM2-l.)/((AIM2**2)*Q) 

651  CONTINUE 

DO  101  N  =  3,  300 
MLESS1  =  M-l 
DO  130  I  =  1,M 
130   DL(I)  =  DM(I)  +  VM(I)*12.*T 
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DO  131  I  =  l.MLESSl 

131  CL(I)  =  DL(I)  -  DL(I+1) 
13   IF(DL(M))  12,12,1U 

12   DE  -  0.0 

GO  TO  2k 
Ik      IF(DL(M)  -  Q)  16,16,18 
16   DE  =  0.0 

GO  TO  2k 
18   VALUE  =  DL(M)  -  Q 

IF  (VALUE-DE)  22,22,20 
20   DE  =  VALUE 

GO  TO  2k 
22   DE  =  DE 
2k       DP(N)  =  DE 

R(M)  =  (DL(M)-DP(N))#SP(M)*(1.+Z*VM(M)) 

IF  (DP(N))  25,25,26 
26   IF  (DP(N)  -  DP(N-l))  197,197,25 
25   CONTINUE 
3   VAL  =  S(1)*CL(1) 

IF(CL(1)-CM(1))  5,5,U 
k       FL(1)  =  VAL 

GO  TO  33 
5   CL(1)  =  CM(1) 

FL(1)*VAL/(RES1**2  )-(l./(RESl*-*2  )-l.  )*S(l)*CL(l) 
33   VAL2  =  S(2)*CL(2) 

IF(CL(2)-CM(2))  35,35,3l| 
3U   FL(2)  =  VAL2 

GO  TO  37 

35  CL(2)  ■  CM(2) 

FL(2)  =  VAL2/(RES2^x-2)-(l./(H5S2^2)-l.)--S(2)-x-CL(2) 

37  IF(CAP)  38,36,36 

38  IF(FL(2))  39,36,36 

39  FL(2)  =  0.0 

36  CONTINUE 

7  IF(FL(1))  8,133,133 

8  FL(1)  =  0.0 

133   DO  132  I  =  3?MLESS1 

132  FL(I)  =  CL(I)-x-S(l) 
IF(FL(2))213,211,210 

213   IF(FTEN1-FL(2))211, 211, 2lU 
21U   FTEN1  =  FL(2) 
GO  TO  211 

210  IF(FMAX1-FL(2))215,211,211 

215  FFiAXl  =  FL(2) 

211  CONTINUE 

MC  =  M/2  +  1 
IF(FL(MC))216,217,218 

216  IF(FTEN2-FL(MC))217,217,219 

219  FTEN2  =  FL(MC) 
GO  TO  217 

218   IF(FMAX2-FL(MC ) )220, 217, 217 

220  FMAX2  =  FL(MC) 
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217 

CONTINUE 
IF(FL(IMl))221,222,22l* 

221 

IF(FTEN3-FL(IM1))222,222,223 

223 

FTEN3  =  FL(IMl) 
GO  TO  222 

221* 

IF(FMAX3-FL(IM1))225,222,222 

225 

FMAX3  =  FL(IMl) 

222 

CONTINUE 

U2 

IF(RUB)U9,1*9,1*U 

kh 

DO  U8  I  =  3,MLESS1 

DPSL(I)  =  DPSM(I) 

CHANGE  =  DL(I)  -  Q 

SUM  =  DL(I)  +  Q 

IF  (DPSL(  I )  -CHANGE  )U5 ,  1*6,  1*6 

h$ 

DPSL(I)  =  CHANGE 

1*6 

IF(DPSL(I)-SUM)U8,U8,i*7 

1*7 

DPSL(I)  =  SUM 

U8 

R(I)  =  (DL(I)  -  DPSL(I))*SP(I)*(1.+ZPWM(I)) 

k9 

DO  50  I  =  1,M 

VL(I)  =  VM(I)  +  (FL(I-l)  -  FL(I)  -  R(I))*T*32.2/W(l) 

50 

CONTINUE 

DO  75  K  =  1,M 

X  =  1 

IF  (VL(K))     75,75,800 

75 

CONTINUE 
GO  TO  199 

800 

CONTINUE 

IF(VL(2)  »-  CHECK) 99, 99, 190 

99 

IF(VL(M)-CHECK)l02, 102,191* 

102 

DO  98  K  =  1,M 
DS(K)  =  DM(K) 
DM(K)  =  DL(K) 
DL(K)  =  0,0 
V3(K)  =  VM(K) 
VM(K)  =  VL(K) 
VL(K)  =  0.0 
CS(K)  =  CM(K) 
CM(K)  =  CL(K) 
CL(K)  =  0.0 
FS(K)  =  FM(K) 
FM(K)  =  FL(K) 
FL(K)  =  0.0 
DPSS(K)  =  DPSM(K) 
DPSM(K)  =  DPSL(K) 

98 

DP3L(K)  =0.0 

101 

CONTINUE 
GO  TO  196 

190 

WRITE  OUTPUT  TAPE  6,191,N 

191 

0  FORMAT  (6IiH  VELOCITY  OF  PILE  CAP  EXCEEDED  TWICE  THE  RAM 

1  HEN  N  WAS  13) 
GO  TO  196 

VELOCITY  W 

1U; 


19h      WRITE  OUTPUT  TAPE  6,195,N 

195  0  FORMAT  (610.  VELXITI  OF  PILE  TIP  EXCEEDED  TWICE  THE  RAM  VELOCITY  W 
1  HEN  N  WAS  13) 

GO  TO  196 

197  WRITE  OUTPUT  TAPE  6,1?8,N 

198  FORMAT  (26H  DP  BECAME  CONSTANT  AT  N  =13) 
GO  TO  196 

199  WRITE  OUTPUT  TAPE  6,200,N 

200  FORMAT  (52H  ALL  VL  WERE  SIMULTANEOUSLY  NEGATIVE  OR  ZERO  AT  N  -  Ik) 

196  CONTINUE 

WRITE  OUTPUT  TAPE  6,607,RU 
607       FORMAT  (27H  ULTIMATE  RESISTANCE  (RU)  =  E15.U) 

WRITE  OUTPUT  TAPE  6,606, (SP(L).L=1,M) 
606       FORMAT  (17H  NEW  VALUES  OF  SP/(8El5.U)) 

WRITE  OUTPUT  TAPE  6,305,DP(N) 

305  FORMAT  (7H  SET  =  E10.U) 
BLOW  -  l./DP(N) 

WRITE  OUTPUT  TAPE  6, 805, BLOW 
805       FORMAT (66H 

1  BLOWS/LN=  E10.3) 

FTEN1  =  FTENl/AREAl 

FTEN2  =  FTEN2/AREA2 

FTEN3  =  FTEN3/AREA3 

WRITE  OUTPUT  TAPE  6,307,FTEN1,FTEN2,FTEN3 
307  0  FORMAT (30H  MAXIMUM  TENSION  AT  HEAD  =         E15.U/ 

1  30H  MAXIMUM  TENSION  AT  CENTER  -     E15.1*/ 

2  30H  MAXIMUM  TENSION  AT  TIP  ■  E15.U) 
FMAX1  =  FMAXl/AREAl 

FMAX2  ■  FMAX2/AREA2 
FMAX3  ■  FMAX3/AREA3 
WRITE  OUTPUT  TAPE  6,306,FMAX1,FMAX2,FMAX3 

306  0  FORMAT  (26H  MAXIMUM  FORCE  AT  HEAD  =       E15.U/ 

1  26H  MAXIMUM  FORCE  AT  CENTER  =  E15.U/ 

2  26H  MAXIMUM  FORCE  AT  TIP  =         E15.H//) 
602       CONTINUE 

RU  «  RU  +  UOOOO. 
600      CONTINUE 

IF(ITURN)205,205,207 
205       CALL  EXIT 

END 
*  DATA 
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APPENDIX  D 
DETAILS  OF  THE  "HILEY  FORMULA  COMPUTER  PROGRAM" 

General  t  This  program  is  written  in  Fortran  Language  for  the 
IBM  7090  computer,  but  may  be  used  equally  well  with  the  709  or  709H 
computer.  It  has  as  its  objective  the  systematic  production  of  data  of 
set  vs.  ultimate  ground  resistance  using  the  Hiley  type  formula.  The 
program  will  take  given  pile  and  hammer  information  and  beginning  with 
an  ultimate  ground  resistance  of  $   tons  will  range  in  increments  of 
5  tons  until  an  ultimate  ground  resistance  of  2l|0  tons  is  reached.  At 
each  value  of  ground  resistance  the  set  is  computed  and  written  out. 
The  program  is  arranged  so  that  any  type  of  hammer  or  pile  can  be  used. 
The  coefficients  for  use  with  the  Hiley  formula  are  as  recommended  by 
Chellis  (ref .  $)   and  are  built  into  the  program  except  as  indicated 
under  input  data  preparation  section  below. 

Termination :  No  conditions  for  termination  are  built  into 
this  program.  If  extensive  use  were  planned  for  this  program,  an  "if" 
statement  should  be  added  so  that  it  will  cause  the  calculations  to 
terminate  when  the  computed  set  becomes  negative. 

Input  Data  Preparation;  All  data  is  punched  on  standard  data 

cards  in  the  exact  order  indicated  below.  Upon  completion  of  punching 

the  data  deck,  it  is  placed  immediately  after  the  #  DATA  card  of  the 

source  program. 

CARD  ORDER        DESCRIPTION 

1  If  another  set  of  data  follows  this  set  for 

another  case,  punch  +1  in  spaces  9  and  10;  if 
not,  punch  -1.  Use  no  decimals. 
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2  If  the  hammer  is  double  acting,  punch  +1  in 
spaces  9  and  10;  if  not,  punch  -1.  Use  no 
decimals . 

3  If  the  pile  is  steel,  punch  +1  in  spaces  9  and 
10;  if  wood,  -1  in  spaces  9  and  10;  if  concrete, 
punch  0  in  space  10.  Use  no  decimals, 

h  If  triangular  side  distribution  is  desired,  punch 

+1  in  spaces  9  and  10;  if  rectangular  side  dis- 
tribution is  desired,  purch  -1  in  spaces  9  and 
10. 

5  In  the  first  10  spaces  punch  the  hammer  effi- 
ciently using  a  decimal  point. 

6  In  the  first  10  spaces  punch  the  weight  of  the 
ram  in  pounds  using  a  decimal  point. 

7  In  the  first  10  spaces  punch  the  coefficient 
of  restitution  using  a  decimal  point. 

8  In  the  first  10  spaces  punch  the  weight  of  the 
pile  including  the  shoe  and  driving  cap  for 
drop  hammers  and  single  acting  steam  hammers. 
Weight  of  pile  including  weight  of  anvil  in  case 
of  double  acting  or  differential  acting  steam 
hammers.  All  in  pounds  and  with  a  decimal  point. 

9  In  the  first  10  spaces  punch  the  cross  sectional 
area  in  square  inches  of  the  pile  at  its  head 

or  the  area  of  the  capblock  if  one  is  used.  Use 
a  decimal  point. 

10  Same  as  card  9  but  for  the  cross  sectional  pile 
area  at  mid-length  of  the  pile. 

11  Same  as  card  9  but  for  the  cross  sectional  pile 
area  at  the  pile  tip. 

12  Punch  in  the  first  10  spaces  the  ground  quake 
in  inches,  using  a  decimal  point. 

13  Punch  in  the  first  10  spaces  the  modulus  of 
elasticity  in  pounds  per  square  inch  using  a 
decimal  point. 

Ik  Punch  in  the  first  10  spaces  the  pile  length  in 

feet  using  a  decimal  point. 
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15  Punch  in  the  first  10  spaces  the  height  of  free 

fall  of  the  ram,  in  inches,  for  drop  hammers 

or 
the  normal  (shortest)  stroke  of  ram  in  inches 
for  single  acting  steam  hammers 

or 
for  double  acting,  differential  acting  steam  and 
diesel  hammers  use  the  rated  energy  per  blow  in 
foot  pounds  as  published  by  the  manufacturers. 

Running  the  Program;  The  data  deck  of  cards  just  described  is 
placed  immediately  after  the  source  program  deck  and  is  ready  for 
delivery  to  a  machine  operator  for  processing. 

Results :  In  the  output  will  be  repeated  the  original  input 
information,  values  of  the  ultimate  ground  resistance,  per  cent  of 
load  carried  by  the  point,  and  the  set.  The  results  are  printed  in 
decimal  form  with  a  multiplier  so  that  0.6192  E00,  for  example,  means 
0.6192  x  10°,  or  0.6192. 

The  Program:  The  complete  program  is  presented  so  that  any  who 
wish  to  use  it  may  do  so.  Each  line  of  the  program  should  be  punched 
on  a  separate  card  and  assembled  in  the  exact  order  indicated.  Anyone 
preparing  cards  from  this  print-out  is  cautioned  to  verify  the  program 
wi th  test  cases  before  relying  on  the  results  obtained. 
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*ID   REESE        C2165   CIVIL  ENGR   G  HILET     7  150 

*  XEQ 

*  NOBIN 

*  LABEL 

C  HILEY  FORMULA  PROGRAM 

206  READ  INPUT  TAPE  5,20,ITURN,IEQUIP,IMAT,MUD 

20  FORMAT (110) 

READ  INPUT  TAPE  5,21,E,WR,RES,WP,APH,APAV,APT,Q,ELAS,PL,H 

21  FORMAT  (F10.U) 
IF(IMAT)  1,3,5 

1  WRITE  OUTPUT  TAPE  6,2 

2  FORMAT  (12H  WOODEN  PILE) 
GO  TO  7 

3  WRITE  OUTPUT  TAPE  6,U 

h  FORMAT  (liiH  CONCRETE  PILE) 
GO  TO  7 

5  WRITE  OUTPUT  TAPE  6,6 

6  FORMAT (11H  STEEL  PILE) 

7  IF(IEQUIP)  8,12,10 

8  WRITE  OUTPUT  TAPE  6,9 

9  FORMAT (3I4H  SINGLE  ACTING  OR  DROP  HAMMER  USED) 
GO  TO  12 

10  WRITE  OUTPUT  TAPE  6,11 

11  FORMAT (26H  DOUBLE  ACTING  HAMMER  USED) 

12  IF(MUD)  13,205,15 

13  WRITE  OUTPUT  TAPE  6,lU 

1U  FORMAT(33H  SIDE  DISTRIBUTION  IS  RECTANGULAR) 
GO  TO  17 

15  WRITE  OUTPUT  TAPE  6,16 

16  FORMAT  (32H  SIDE  DISTRIBUTION  IS  TRIANGULAR) 

17  WRITE  OUTPUT  TAPE  6,l8,E 

18  FORMAT (13H  EFFICIENCY-  E10.U) 
WRITE  OUTPUT  TAPE  6,19,WR 

19  FORMAT  (16H  WEIGHT  OF  RAM-  E10.U) 
WRITE  OUTPUT  TAPE  6, 22, RES 

22  FORMAT  (29H  COEFFICIENT  OF  RESTITUTION-  E10.U) 
WRITE  OUTPUT  TAPE  6,23,WP 

23  FORMAT  (13H  WT  OF  PILE-  E10.U) 
WRITE  OUTPUT  TAPE  6,2ii,APH 

2U  FORMAT  (26H  AREA  PILE  HEAD  IN  SQ  IN-  E10.U) 
WRITE  OUTPUT  TAPE  6,25,APAV 

25  FORMAT  (20H  AVERAGE  PILE  AREA=  E10.U) 
WRITE  OUTPUT  TAPE  6, 26, APT 

26  FORMAT  (16H  AREA  PILE  TIP-  E10.U) 
1-JRITE  OUTPUT  TAPE  6,27.Q 

27  FORMAT  (8H  QUAKE-  E10.U) 
WRITE  OUTPUT  TAPE  6,28,ELAS 

28  FORMAT  (17H  YOUNGS  MODULUS-  E10.U) 
WRITE  OUTPUT  TAPE  6, 29, PL 

29  FQ..MAT  (20H  PILE  LENGTH  IN  FT-  E10.U) 
IF  (IEQUIP)  31,35,33 
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31  WRITE  OUTPUT  TAPE  6,32,H 

32  FORMAT  (2I4H  HAMMER  DROP  IN  INCHES=  EIO.U) 

GO  TO  35 

33  WRITE  OUTPUT  TAPE  6,  3U,  H 

3U  FORMAT  (UOH  RATED  ENERGY  OF  HAMMER  IN  FOOT  POUNDS =  EIO.U) 
35  RU  =  lOOOO. 

DO  20U  J  =  1,U8 

WRITE  OUTPUT  TAPE  6,30,RU 
30  FORMAT  (20H  THE  VALUE  OF  RU  IS  EIO.U//) 

PI  -  RU/APH 

P2  =  RU/APAV 

P3  =  RU/APT 

DO  203  I  =  1,3 

IF(I-l)  60,UO,U2 
UO  PLMOD  =  PL 

WPC  =  WP/2. 

WRITE  OUTPUT  TAPE  6,1a 
U.  FORMAT  (29H  ALL  OF  LOAD  CARRIED  BY  POINT) 

GO  TO  60 
U2  IF(I-2)  60,UU,52 
UU  WRITE  OUTPUT  TAPE  6,  U7 
U7  FORMAT (31H  ONE  HALF  LOAD  CARRIED  BY  POINT) 

IF (MUD)  U6,60,U8 
U6  PLMOD  =0.75  *PL 

WPC  =  WP 

GO  TO  60 
U8  PLMOD  =  U.  *  PL  /  6. 

WPC  -  WP 

GO  TO  60 
52  IF(I-3)  60,  5U,  60 
5U  WRITE  OUTPUT  TAPE  6,  55 

55  FORMAT(29H  ALL  LOAD  CARRIED  BY  FRICTION) 
IF (MUD)  56,60,^8 

56  PLMOD  =  PL/2 
WPC  =  WP 

GO  TO  60 
58  PLMOD  =  2.  -x-  PL  /  3. 

WPC  =  WP 
60  C3  -  Q 

IF(IMAT)  70,80,90 
70  CI-  0.00010  *  PI 

C2  =  0.00008*P2-*PLMOD*l5  00000.  /ELAS 

GO  TO  100 
80  CI-  0.00025  *  PI 

C2  -  O.OOOOOU*P2*PLMOD*3000000./ELAS 

GO  TO  100 
90  CI  =  0.00008  -x-  PI 

C2  =  O.OOOOOOU-x-P2-x-PLMOD-x-30000000./ELAS 
100  DP= ( ( (E-x-WR-x-H )-x-(WR+WPC-x-RES-x-x-2 ) )/( (WR+WPC  )-x-RU ) )- ( C1+C2+C3  )/2 . 

LF(IEQUIP)  10ii,203,l°2 
102  DP  =  DP  +  (Cl+C2+C3)/2. 
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DP  =  12,-x-DP 
DP  =  DP  -(Cl+C2+C3)/2. 
10U  BLOW  =  l./DP 

WRITE  OUTPUT  TAPE  6, 201, DP 

201  F0RMAT(8H  SET  IS  E20.U) 
WRITE  OUTPUT  TAPE  6, 202, BLOW 

202  FORMAT (2 OH  BLOWS  PER  INCH  ARE  E20.U/) 

203  CONTINUE 

20U  RU  =  RU  +  10000. 

IF  (ITURN)  20^,205,206 
205  CALL  EXIT 
END 
*     DATA 


. 
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APPENDIX  E 

This  appendix  includes  some  of  the  graphs  of  computer  solutions 
which  were  used  both  for  correlation  with  actual  pile  cas^es  and  for 
studying  the  effect  of  varying  the  size  of  hammer,  length  of  pile, 
ground  quake,  point  damping,  and  side  damping.  On  each  sheet  several 
curves  are  drawn  for  various  values  of  frictional  resistance  expressed 
as  a  percentage  of  the  total  ultimate  bearing  capacity  of  the  pile  and 
indicated  as  follows: 

Curve  A  —  End  bearing;  no  side  friction 

Curve  B  —  75  per  cent  end  bearing;  25  per  cent  side  friction. 

Curve  C  —  £0  per  cent  end  bearing;  50  per  cent  side  friction. 

Curve  D  —  25  per  cent  end  bearing;  75  per  cent  side  friction. 

Curve  E  —  10C$  friction. 
Rectangular  side  distribution  was  used.  The  description  of  the  driving 
equipment,  pile,  and  soil  properties  for  which  the  computations  were 
made  appear  on  each  sheet,  but  for  convenience  of  use  in  future  correla- 
tion work  an  index  is  included  which  summarizes  this  information  for 
each  sheet.  The  figure  numbers  and  page  numbers  are  identical  in  this 
appendix. 
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INDEX  OF  FIGURES  IN  APPENDIX  E 

1.  Piles  Driven  by  a  Vulcan  #1  Hammer 

a.  UO-foot  12BP53  Piles: 

Figure 
Point  Damping   Side  Damping   Ground  Quake        Number 

El 

E2 

E3 
EU 

E5 
E6 

E? 
E8 
E9 
E10 
Ell 
E12 
E13 
ELI* 


ElS 

E16 
E17 
E18 
E19 
E20 


0.0 

0.0 

0.10 

0.10 

0.033 

0.10 

0.15 

0.0 

0.10 

0.20 

0.067 

0.10 

0.30 

0.10 

0.10 

o.Uo 

0.133 

0.10 

0.50 

0.167 

0.10 

0.10 

0.033 

0.20 

0.20 

0.067 

0.20 

0.30 

0  10 

0.20 

o.Uo 

0.133 

0.20 

o.5o 

0.167 

0.20 

0.10 

0.033 

0.30 

0.20 

0.067 

0.30 

b.  80-foot  12BP53  Piles 

• 
• 

0.0 

0.0 

0.10 

0.10 

0.033 

0.10 

0.20 

0.067 

0.10 

0.30 

0.10 

0.10 

o.Uo 

0.133 

0.10 

o.5o 

0.167 

0.10 
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(80-foot  12BP53  Piles  cont'd. ) 


Point  Damping 

Side  Damping 

Ground  Quake 

Figure 
Number 

0.10 

0.033 

0.20 

E21 

0.20 

0.067 

0.20 

E22 

0.30 

0.10 

0.20 

E23 

o.Uo 

0.133 

0.20 

E2U 

0.50 

0.167 

0.20 

E25 

0.10 

0.033 

0.30 

E26 

0.20 

0.067 

0.30 

E27 

c.  120-foot  12BP53  Piles 

0.0 

0.0 

0.10 

E28 

0.10 

0.033 

0.10 

E29 

0.20 

0.067 

0.10 

E30 

0.30 

0.10 

0.10 

E31 

O.UO 

0.133 

0.10 

E32 

o.5o 

0.167 

0.10 

E33 

0.10 

0.033 

0.20 

E3U 

0.20 

0.067 

0.20 

E35 

0.30 

0.10 

0.20 

E36 

o.Uo 

0.133 

0.20 

E37 

0.50 

0.167 

0.20 

E38 

0.10 

0.033 

0.30 

E39 

0.20 

0.067 

0.30 

EUO 

d.  18-inch-square  concrete  pile  UO  feet  long: 

0.15         0.0  0.10  EUl 


2$h 


2.  Piles  Driven  with  a  Vulcan  0  Hammer 

a.  12BP53  Piles  80  feet  long: 

Point  Damping   Side  Damping   Ground  Quake 

0.20  0.067  0.10 

0.U0         0.133  0.10 

b.  12BP53  Piles  16   feet  long: 

0.10         0.033  0.10 

c.  1UBP89  Piles  U$  feet  long: 

0.15  0.0$  0.10 

0.20  0.067  0.10 

0.20  0.067  0.20 

d.  1UBP117  Piles  h$   feet  long: 

0.15  0.05  0.10 

0.20  0.067  0.20 


Figure 
Number 

EU2 

EU3 


EU6 
EU7 

EI48 
EU9 
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